#include <stdio.h>
#include <string.h>
#include <stdlib.h>

int main(void)
{
	// 入出力ファイルに利用 
	char inputf[255];
	char outf[255];
	FILE *fpi;
	FILE *fpo;

	/*
	* データを保存する配列、構造体の宣言
	* 詳細はそれぞれのコメント文を参照の事
	* なお、それぞれの項目は、FORMMBY.TXTの用語を使用している。
	*/

	struct mesh_head {	// メッシュヘッダレコード用
		int mesh;	// 2次メッシュコード
		int node;	// ノード数
		int line;	// ライン数
		int area;	// エリア数
		int nokori;	// レコード数
	};

	struct area_record {				// エリアレコード用
		int gyousei;			// データ項目コード    
		char todouhuken[10];	// 都道府県名
		char gunshi[20];		// 郡市・北海道の支庁名
		char kutyou[20];		// 政令指定都市・区町村名
	};
	
	struct mesh_head syousai[100];	// メッシュヘッダレコード用の構造体
	struct area_record naiyo[1000];	// エリアレコード用の構造体
	
	int rinsetsu[5000][2];	// ラインレコード用の配列
 
	char str[80];  // 読み込み用
	char buf[80];  // 処理用
	
	int i, j, k;  // カウント用

	int countMESH;	// 表示位置の指定
	int countLINE;	// 表示位置の指定
	int countAREA;	// 表示位置の指定
	int nagasa;		// kuwariの整列用
	
	int min, t, max; //ソート時に利用

	int kuwari[5000][30];
	int a[100][30];

	/* 配列kuwari[][]の初期化 */
	
	for(i=0; i<5000; i++){
		for(j=0; j<30; j++){
			kuwari[i][j] = 0;
		}
	}
	
	for(i=0; i<100; i++){
		for(j=0; j<30; j++){
		a[i][j]=0;
		}
	}
	
	/* ファイルの処理 */
	printf("入力ファイル名を入力してください\n");
	scanf("%s",inputf);
	printf("出力ファイル名を入力してください\n");
	scanf("%s",outf);

	fpi = fopen(inputf, "r");
	if ( fpi == NULL) {
		printf("ファイルをオープンできませんでした。\n");
		return 0;
	}

	fpo = fopen(outf, "w");
	if ( fpo == NULL) {
		printf("ファイルをオープンできませんでした。\n");
		return 0;
	}
	
	i = 0;
	j = 0;
	k = 0;
	
	/* データの読み込み */
	while ( fgets (str, 80, fpi) != NULL) {
		strncpy(buf, str, 2);
		buf[2] = '\0';

		if (strcmp(buf,"M ") == 0){  //メッシュヘッドレコードの場合
			strncpy(buf, str+2, 6);
			buf[6] = '\0';
			syousai[i].mesh = atoi(buf);
			
			strncpy(buf, str+2+6+20+3, 5);
			buf[5] = '\0';
			syousai[i].node = atoi(buf);

			strncpy(buf, str+2+6+20+3+5, 5);
			buf[5] = '\0';
			syousai[i].line = atoi(buf);
      
			strncpy(buf, str+2+6+20+3+5+5, 5);
			buf[5] = '\0';
			syousai[i].area = atoi(buf);
      
			strncpy(buf, str+2+6+20+3+5+5+5+5, 5);
			buf[5] = '\0';
			syousai[i].nokori = atoi(buf);

			i++;
		}
		else if (strcmp(buf,"L ") == 0){  //ラインレコードの場合
			strncpy(buf, str+2+2+2+5+6+5+1+5+1, 5);
			buf[5] = '\0';
			rinsetsu[j][0] = atoi(buf);
			rinsetsu[j+1][1] = atoi(buf);

			strncpy(buf, str+2+2+2+5+6+5+1+5+1+5+5, 5);
			buf[5] = '\0';	 
			rinsetsu[j][1] = atoi(buf);
			rinsetsu[j+1][0] = atoi(buf);
			
			j=j+2;
		}
		else if (strcmp(buf,"A ") == 0){  //エリアレコードの場合
			strncpy(buf, str+2+2, 5);
			buf[5] = '\0';
			naiyo[k].gyousei = atoi(buf);

			strncpy(buf, str+2+2+5+5+5+5+4+4, 8);
			buf[8] = '\0';
			strcpy(naiyo[k].todouhuken, buf);

			strncpy(buf, str+2+2+5+5+5+5+4+4+8, 16);
			buf[16] = '\0';
			strcpy(naiyo[k].gunshi, buf);

			strncpy(buf, str+2+2+5+5+5+5+4+4+8+16, 16);
			buf[16] = '\0';
			strcpy(naiyo[k].kutyou, buf);

			k++;
		}
	}

	/* 画面上への出力 */
	countMESH = i;
	countLINE = j/2;
	countAREA = k;

	printf("出力情報\n");
	printf("メッシュヘッダレコード数は%dです。\n",countMESH);
	printf("ラインレコード数は%dです。\n",countLINE);
	printf("エリアレコード数は%dです。\n\n",countAREA);
	printf("-----\n");

	countLINE=0;
	countAREA=0;

	for(i=0; i<countMESH; i++){
		/* メッシュヘッドレコードの出力 */
		printf("Mesh:%d\n",syousai[i].mesh);
		printf("Node:%d\n",syousai[i].node);
		printf("Line:%d\n",syousai[i].line);
		printf("Area:%d\n",syousai[i].area);

		/* ラインコードの出力 */
		printf("\n");
		for (j=0; j<syousai[i].line; j++){
			printf("左:%d\t",rinsetsu[countLINE][0]);
			printf("右:%d\n",rinsetsu[countLINE][1]);
			countLINE = countLINE+2;
		}
		printf("\n");

		/* エリアレコードの出力 */
		for(k=0; k<syousai[i].area; k++){
			printf("データ項目コード:%d\n",naiyo[countAREA].gyousei);
			printf("都道府県名:%s\n",naiyo[countAREA].todouhuken);
			printf("群市:%s\n",naiyo[countAREA].gunshi);
			printf("区町村:%s\n",naiyo[countAREA].kutyou);
			countAREA++;
		}
		printf("\n");
	}
	/* 画面上への出力　終了 */

	/* 
	隣接関係の読み取り
	rinsetsuからkuwariにデータを移す。
	その際に、不必要なデータ及び、重複が起こらないようにする。
	*/
	nagasa = -1;
	
	for ( i=0; i<5000; i++){
		k = 0;
		if((rinsetsu[i][0] != 0) && (rinsetsu[i][0] != 88888) && (rinsetsu[i][0] != 99999)){
			nagasa++;
			kuwari[nagasa][k] = rinsetsu[i][0];
			if((rinsetsu[i][1] != 88888) && (rinsetsu[i][1] != 99999)){
				k++;
				kuwari[nagasa][k] = rinsetsu[i][1];
			}
		}

		j=i+1;
		for ( j; j<5000; j++){
			if(kuwari[nagasa][0] == rinsetsu[j][0]){
				if((rinsetsu[j][1] != 88888) && (rinsetsu[j][1] != 99999)){
					k++;
					kuwari[nagasa][k] = rinsetsu[j][1];
				}
				rinsetsu[j][0] = 0;
			}
		}
	}
	/* 隣接関係の読み取り 終了 */
	
	/* 並び替えの処理を以下で行う。
	並び替えた結果はa[][]に記憶される。
	はじめにkuwari[][0](列)の項目で並び替え、
	その後にkuwair[x][項目数](行)の項目で並び替える。
	*/

	// 事前準備
	for(i=0; i< nagasa; i++){
		a[i][1]=kuwari[i][0];
	}

	// 行のソート
	for (i=0; i < nagasa-1; i++){
		min = i;
		for (j = i+1; j <= nagasa-1; j++){
			if (a[j][1] < a[min][1]){
				min =j;
			}
		}
		t = a[min][1]; a[min][1] = a[i][1]; a[i][1] = t;
	}
	
	// 列の項目数の計算、及びa[][]へのデータコピー
	for(i=0; i<nagasa; i++){
		for(j=0; j<nagasa; j++){
			if(a[i][1] == kuwari[j][0]){
				for(k=0; k<30; k++){
					if(kuwari[j][k] == 0){
						a[i][0] = k;
						break;
					}
					else{
						a[i][k+1] = kuwari[j][k];
					}
				}
			}
		}
	}

	// 列のソート
	for(k=0; k<nagasa; k++){
		for (i=2; i < a[k][0]; i++){
			min = i;
			for (j = i+1; j <= a[k][0]; j++){
				if (a[k][j] < a[k][min]){
					min =j;
				}
			}
			t = a[k][min]; a[k][min] = a[k][i]; a[k][i] = t;
		}
	}

	//同一列内での同一している数字の削除
	for(i=0; i<nagasa; i++){
		for(j=2; j<a[i][0]; j++){
			if(a[i][j] == a[i][j+1]){
				for(k=j+1;k<a[i][0];k++){
					a[i][k]=a[i][k+1];
				}
				a[i][0]--;
			}
		}
	}

	// 最大表示幅
	max = a[i][0];
	for (i=1; i <= nagasa; i++){
		if(max<a[i][0]){
			max = a[i][0];
		}
	}
	printf("max:%d\n",max);
	/* 隣接関係の解析 終了 */

	/* 隣接関係の解析の結果*/
	for(i=0; i<nagasa; i++){
		for(j=0; j<30; j++){
			printf("%d,",kuwari[i][j]);
		}
		printf("\n");
	}

	printf("-----\n");

	for(i=0; i<nagasa; i++){
		printf("%d,==,",a[i][1]);
		for(j=2; j<=a[i][0]; j++){
			printf("%d,",a[i][j]);
		}
		if(max-a[i][0] > 0){
			for(k=0; k<max-a[i][0]; k++){
				printf("0,");
			}
		}
		printf("\n");
	}
	/* 隣関係の解析の結果 終了 */
	
	for(i=0; i<nagasa; i++){
		fprintf(fpo,"%d,==,",a[i][1]);
		for(j=2; j<=a[i][0]; j++){
			fprintf(fpo,"%d,",a[i][j]);
		}
		if(max-a[i][0] > 0){
			for(k=0; k<max-a[i][0]; k++){
				fprintf(fpo,"0,");
			}
		}
		fprintf(fpo,"\n");
	}

	/* 終了処理 */
	fclose(fpi);
	fclose(fpo);
	return 0;
}