2Â¥: Originally posted by Mr__Right at 2017-08-18 11:00:47
Ê×ÏÈÄãҪʶÊý£¬ÄÇÊÇÒ»¹²9¸ö£¬ÏàÁÚµÄÊÇ8¸ö¡£

È»ºó£¬ÕÒµ½ÖÐÐÄλÖÃÖ¸±ê£¬



È»ºó,¸ù¾ÝÖ¸±êÔö¼õµÄ2^3=8ÖÖÅÅÁÐ×éºÏ£¬È«²¿¸³ÖµÎª0

3Â¥: Originally posted by imyourkobe at 2017-08-18 21:50:25
ºÜÈÝÒ׵ġ£ÓöþάÊý×é»òÈÝÆ÷À´´æ´¢Êý¾Ý£¬Ð´Ò»¸ö¼òµ¥µÄËã·¨ÅжÏÁù¸öÁ¬ÐøµÄÊý¾ÝΪ1µÄ¸öÊýÀÛ¼Ó¡£Âú×ãÔòÐ޸ġ£²»Âú×ã¼ÌÐøÅжÏÏÂÒ»×é¡£

2Â¥: Originally posted by Mr__Right at 2017-08-18 11:00:47
Ê×ÏÈÄãҪʶÊý£¬ÄÇÊÇÒ»¹²9¸ö£¬ÏàÁÚµÄÊÇ8¸ö¡£

È»ºó£¬ÕÒµ½ÖÐÐÄλÖÃÖ¸±ê£¬



È»ºó,¸ù¾ÝÖ¸±êÔö¼õµÄ2^3=8ÖÖÅÅÁÐ×éºÏ£¬È«²¿¸³ÖµÎª0

#define _CRT_SECURE_NO_WARNINGS

#include <stdio.h>
#include <math.h>
#include <iostream>
#include <fstream>
#include <string>  
#include <iostream>  
#include <vector>
using namespace std;

struct Location {
        double fai;                       
        int *grain_number;
}num[1000000];

#define seed 40000 //Á£×Ó¸öÊý
bool main(){
        int grain_num[seed], correction_num[seed], grain_size[seed],
                divided_grain_count[seed], div_area,
                grain_count, grain_sum, divided_grain_sum, cut_off,
                grain_size_distribution[seed], max_grain_size, div_distribution;

        double boundary_line, boundary_line1, boundary_line2;
        int count, count2, meshx, meshy, meshX, meshY;
        char        fname[100];
        FILE        *datin0;
        //        ³õÆÚ»¯
        grain_count = 0, count2 = 0, grain_sum = 0, divided_grain_sum = 0, count = 0;
        max_grain_size = 0;
        for (int i = 0; i < seed; i++){
                grain_num[i] = i;
                correction_num[i] = i;
                grain_size[i] = -1;
                divided_grain_count[i] = -1;
                grain_size_distribution[i] = 0;
        }
        // Ö¸¶¨Îļþ¼ÐºÍ±£´æ
        char*  file_name("C:\\testp.loc");
        ofstream fout1("save1.xls");
        ofstream fout2("save2.xls");
        ofstream fout3("save3.xls");
        ofstream fout4("save4.xls");
        ofstream fout5("save5.xls");

        boundary_line = 0.4;        // ÅжÏÊýÖµ´óС
        cut_off = 0;        // cut_offÒÔϵIJ»¼ÆËã

        //        ±ß½çÉ趨
        meshx = 1200;
        meshy = 150;
        meshX = meshx + 2; meshY = meshy + 2;

        // ´ò¿ªÎļþ¼Ð
        std::ifstream fin;
        fin.open(file_name, std::ios::in);
        // Îļþ¼ÐÈ·ÈÏ
        if (!fin){
                std::cout << " file open error! " << "\n";
                return false;
        }
        if (!fout1){
                cerr << "open error!" << endl;
                exit(1);
        }
        if (!fout2){
                cerr << "open error!" << endl;
                exit(1);
        }
        if (!fout3){
                cerr << "open error!" << endl;
                exit(1);
        }
        if (!fout4){
                cerr << "open error!" << endl;
                exit(1);
        }
        if (!fout5){
                cerr << "open error!" << endl;
                exit(1);
        }
        //¶ÁÈëÁ£×Ó
        for (int i = 0; i<meshY; i++){
                for (int j = 0; j<meshX; j++){
                        if (i <= 0 || i >= meshY - 1 || j <= 0 || j >= meshX - 1){
                                num[count].fai = 0;
                                count++;
                                continue;

                        }
                        else{
                                fin >> num[count].fai;

                                count++;


                        }
                }
        }

        //        ÅжÏÊýÖµ´óС´óÓÚ0.4µÈÓÚ1СÓÚ0.4µÈÓÚ0
                count = 0;
        for (int i = 0; i<meshY; i++){
                for (int j = 0; j<meshX; j++){
                        if (num[count].fai < boundary_line)
                        {
                                num[count].fai = 0;
                        }
                        else{
                                num[count].fai = 1;
                        }
                        fout1 << num[count].fai << " ";
                        count++;
                }
                fout1 << endl;
        }
        // grain_number
        count = 0;
        for (int i = 0; i<meshY; i++){
                for (int j = 0; j<meshX; j++){
                        if (i == 0 || i == meshY - 1 || j == 0 || j == meshX - 1){
                                num[count].grain_number = &grain_num[0];
                                count++;
                                continue;
                        }
                        else{
                                if (j == 1){
                                        grain_count++;
                                }
                                num[count].grain_number = &grain_num[grain_count];

                                if (i > 1 && num[count - meshX].fai == 1. && num[count].fai == 1.){
                                        grain_num[grain_count] = *num[count - meshX].grain_number;
                                }
                                if (num[count].fai == 0){
                                        num[count].grain_number = &grain_num[0];
                                        if (num[count + 1].fai == 1 && j>2 && j<meshX - 2){
                                                grain_count++;
                                        }
                                }
                                count++;
                        }
                        fout2 << num[count].grain_number << " "; //ºÅÂ븳ֵ
                        //cout << grain_count << " ";
                }
                fout2 << endl;
        }
        count = 0;
        for (int i = 0; i<10; i++){
                for (int j = 0; j<meshX; j++){

                        if (i == 0 || i == meshY - 1 || j == 0 || j == meshX - 1){
                                count++;
                                continue;
                        }
                        else{
                                if (num[count - meshX].fai == 1 && num[count].fai == 1){            //ÉÏÏÂֵΪ1 ÓÐphase
                                        if (*num[count].grain_number < *num[count - meshX].grain_number){    //Èç¹ûÏÂÃæ·¬ºÅСÓÚÉÏÃæµÄ·¬ºÅ
                                                correction_num[*num[count - meshX].grain_number] = *num[count].grain_number;//°Ñ*num[count].grain_numberµÄÖµ¸³¸øcorrection_num[*num[count - meshX].grain_number]
                                        }
                                        else{
                                                correction_num[*num[count].grain_number] = *num[count - meshX].grain_number;//·ñÔò*num[count - meshX].grain_number¸³¸øcorrection_num[*num[count].grain_number]
                                        }
                                }
                                count++;
                        }
                        fout3 << *num[count].grain_number << " ";
                }
                fout3 << endl;
        }
        // ¼ÆËãgrain_sizeÁ£¾¶´óС
        count = 0;
        for (int i = 0; i<meshY; i++){
                for (int j = 0; j<meshX; j++){
                        if (i == 0 || i == meshY - 1 || j == 0 || j == meshX - 1){
                                count++;
                                continue;
                        }
                        else{
                                grain_size[*num[count].grain_number] += 1;
                                count++;
                        }
                        fout4 << grain_size[*num[count].grain_number] << " ";
                }
                fout4 << endl;
        }
        //ÕâÀï²åÈëÁ£¾¶³¤¶ÈСÓÚ9µÄСÓÚµÈÓÚ0 ´óÓÚ9µÄÖµ²»±ä



        // grain_sum ¼ÆËã¸öÊý
        count = 1;
        while (count<seed){
                if (grain_size[count] > cut_off){
                        grain_sum++;

                        if (grain_size[count] > max_grain_size){
                                max_grain_size = grain_size[count];
                                std::cout << "max = " << max_grain_size << "\n";
                        }
                }
                count++;
        }
}