/*
	切乞 朔亀 繕紫 覗稽益轡

	脊径 督析 : KS 刃失莫 努什闘 督析
	
	厩肯企俳嘘 陳濃斗俳採 悪渋縦
*/
#include <stdio.h>
#include "ksc2kssm.h"
#include "kssm2ksc.h"

char *ch[] = {
	"",
	"@",
	"ぁ",
	"あ",
	"い",
	"ぇ",
	"え",
	"ぉ",
	"け",
	"げ",
	"こ",
	"さ",
	"ざ",
	"し",
	"じ",
	"す",
	"ず",
	"せ",
	"ぜ",
	"そ",
	"ぞ"
};

char *ju[30] = {
	"",
	"",
	"@",
	"た",
	"だ",
	"ち",
	"ぢ",
	"っ",
	"",
	"",
	"つ",
	"づ",
	"て",
	"で",
	"と",
	"ど",
	"",
	"",
	"な",
	"に",
	"ぬ",
	"ね",
	"の",
	"は",
	"",
	"",
	"ば",
	"ぱ",
	"ひ",
	"び"
};

char *jo[30] = {
	"",
	"@",
	"ぁ",
	"あ",
	"ぃ",
	"い",
	"ぅ",
	"う",
	"ぇ",
	"ぉ",
	"お",
	"か",
	"が",
	"き",
	"ぎ",
	"く",
	"ぐ",
	"け",
	"",
	"げ",
	"ご",
	"さ",
	"ざ",
	"し",
	"じ",
	"ず",
	"せ",
	"ぜ",
	"そ",
	"ぞ"
};

unsigned short ksc2kssm(unsigned char c, unsigned char d) {
	if (c < 0xB0 || c > 0xC8 || d < 0xA1 || d == 0xFF)
		return 0;
	return KSSMCODE[c-0xB0][d-0xA1];
}

unsigned short kssm2ksc(unsigned kssm) {
	if (kssm < 0x8861) return 0;
	return KSC_CODE[kssm-0x8861];
}

/*
	data[]税 壱朔亀 授生稽 社特吉 授辞研 index[]拭 煽舌
*/
void set_sort_index(int data[], int index[], int n)
{
	int i, j, temp;

	for (i = 0; i < n; i++) index[i] = i;

	for (i = 0; i < n-1; i++) {	/* 壱朔亀 授生稽 sorting */
		for (j = i+1; j < n; j++) {
			if (data[index[i]] < data[index[j]]) {
				temp = index[i];
				index[i] = index[j];
				index[j] = temp;
			}
		}
	}
}

int main(int argc, char *argv[])
{
	FILE *ifp, *ofp=stdout;
	int c1, c2, cnext1=0, cnext2=0;
	unsigned syl;	// 繕杯莫 製箭

	long sum, nSyl = 0;	// 恥 製箭 鯵呪
	long nSyl2 = 0;		// 恥 尻紗製箭 鯵呪 -- 曽失-段失 朔亀 域至遂
	long cho[21] = {0};
	long jung[30]= {0};
	long jong[30]= {0};

	long chojung[21*30] = {0};
	long jungjong[30*30]= {0};
	long jongcho[30*21] = {0};

	int i, i1, i2, i3, j;
	int si[999];	// 朔亀紺 sorting遂 index 煽舌

	ifp = (argc > 1) ? fopen(argv[1], "r") : stdin;
	if (!ifp) { fprintf(stderr, "File <%s> open failed.\n", argv[1]); return 0; }

	if (argc > 2) {
		if (ofp = fopen(argv[2], "r")) {
			fprintf(stderr, "File <%s> exists. Overwrite it? ", argv[2]); fclose(ofp);
			if (getchar() != 'y') return 0;
		}
		ofp = fopen(argv[2], "w");
		if (!ofp) { fprintf(stderr, "File <%s> open failed.\n", argv[2]); return 0; }
	}

	while (!feof(ifp)) {
		c1 = cnext1 ? cnext1 : fgetc(ifp); if (feof(ifp)) break; cnext1 = 0;

		if (c1 & 0x80) {
			c2  = cnext2 ? cnext2 : fgetc(ifp); if (feof(ifp)) break; cnext2 = 0;

			if (!(syl = ksc2kssm((unsigned char) c1, (unsigned char) c2)))
				continue;

			i1 = (syl >> 10) & 0x1f;	// 段失
			i2 = (syl >>  5) & 0x1f;	// 掻失
			i3 = syl & 0x1f;			// 曽失

			nSyl++;

			cho[i1]++;	// 段失 朔亀
			jung[i2]++;	// 掻失 朔亀
			jong[i3]++;	// 曽失 朔亀

			chojung[i1*30 + i2]++;	// 段失-掻失 朔亀
			jungjong[i2*30 + i3]++;	// 掻失-曽失 朔亀

			cnext1 = fgetc(ifp); if (feof(ifp)) break;
			if (!(cnext1 & 0x80)) { cnext1 = 0; continue; }
			cnext2 = fgetc(ifp); //if (feof(ifp)) break;
			if (!(cnext2 & 0x80)) { cnext2 = 0; continue; }

			if (syl = ksc2kssm((unsigned char) cnext1, (unsigned char) cnext2)) {
				nSyl2++;	// 曽失-段失 尻紗 製箭 鯵呪
				jongcho[i3*21 + ((syl >> 10) & 0x1f)]++;	// 曽失-段失 朔亀
			} else { cnext1 = cnext2 = 0; }
		}
	}

	fprintf(ofp, "切乞 朔亀 繕紫 : 恥 製箭呪 = %ld", nSyl);
	fprintf(ofp, " ( @ : 切乞 蒸製 妊獣 )\n");

// 段失
	fprintf(ofp, "\n***段失 朔亀 : 亜蟹陥 授***\n");
	for (i = sum = 0; i < 21; i++) {
		sum += cho[i];
		if (cho[i]) fprintf(ofp, "%s\t:%8d\t%6.2f\t%6.2f\n",
			ch[i], cho[i], (double) cho[i] / nSyl * 100., (double) sum / nSyl * 100.);
	}

	fprintf(ofp, "\n***段失 朔亀 : 壱朔亀 授***\n");
	set_sort_index(cho, si, 21);	// 壱朔亀 授生稽 sorting
	for (j = sum = 0; j < 21; j++) {
		i = si[j];
		sum += cho[i];
		if (cho[i]) fprintf(ofp, "%s\t:%8d\t%6.2f\t%6.2f\n",
			ch[i], cho[i], (double) cho[i] / nSyl * 100., (double) sum / nSyl * 100.);
	}

// 掻失
	fprintf(ofp, "\n***掻失 朔亀 : 亜蟹陥 授***\n");
	for (i = sum = 0; i < 30; i++) {
		sum += jung[i];
		if (jung[i]) fprintf(ofp, "%s\t:%8d\t%6.2f\t%6.2f\n",
			ju[i], jung[i], (double) jung[i] / nSyl * 100., (double) sum / nSyl * 100.);
	}

	fprintf(ofp, "\n***掻失 朔亀 : 壱朔亀 授***\n");
	set_sort_index(jung, si, 30);	// 壱朔亀 授生稽 sorting
	for (j = sum = 0; j < 30; j++) {
		i = si[j];
		sum += jung[i];
		if (jung[i]) fprintf(ofp, "%s\t:%8d\t%6.2f\t%6.2f\n",
			ju[i], jung[i], (double) jung[i] / nSyl * 100., (double) sum / nSyl * 100.);
	}

// 曽失
	fprintf(ofp, "\n***曽失 朔亀 : 亜蟹陥 授***\n");
	for (i = sum = 0; i < 30; i++) {
		sum += jong[i];
		if (jong[i]) fprintf(ofp, "%s\t:%8d\t%6.2f\t%6.2f\n",
			jo[i], jong[i], (double) jong[i] / nSyl * 100., (double) sum / nSyl * 100.);
	}

	fprintf(ofp, "\n***曽失 朔亀 : 壱朔亀 授***\n");
	set_sort_index(jong, si, 30);	// 壱朔亀 授生稽 sorting
	for (j = sum = 0; j < 30; j++) {
		i = si[j];
		sum += jong[i];
		if (jong[i]) fprintf(ofp, "%s\t:%8d\t%6.2f\t%6.2f\n",
			jo[i], jong[i], (double) jong[i] / nSyl * 100., (double) sum / nSyl * 100.);
	}

// 段失-掻失
	fprintf(ofp, "\n尻紗 切製-乞製 朔亀 繕紫\n");
	fprintf(ofp, "\n***段失-掻失 朔亀 : 亜蟹陥 授***\n");
	for (i = sum = 0; i < 21*30; i++) {
		sum += chojung[i];
		syl = kssm2ksc(0x8000 | (i/30) << 10 | (i%30) << 5 | 0x01);
		if (chojung[i]) fprintf(ofp, "[%s%s:%c%c]\t:%8d\t%6.2f\t%6.2f\n",
				ch[i/30], ju[i%30], syl >> 8 & 0xff, syl & 0xff, chojung[i],
				(double) chojung[i] / nSyl * 100., (double) sum / nSyl * 100.);
	}

	set_sort_index(chojung, si, 21*30);		// 壱朔亀 授生稽 sorting
	fprintf(ofp, "\n***段失-掻失 朔亀 : 壱朔亀 授***\n");
	for (j = sum = 0; j < 21*30; j++) {
		i = si[j];
		sum += chojung[i];
		syl = kssm2ksc(0x8000 | (i/30) << 10 | (i%30) << 5 | 0x01);
		if (chojung[i]) fprintf(ofp, "[%s%s:%c%c]\t:%8d\t%6.2f\t%6.2f\n",
				ch[i/30], ju[i%30], syl >> 8 & 0xff, syl & 0xff, chojung[i],
				(double) chojung[i] / nSyl * 100., (double) sum / nSyl * 100.);
	}

// 掻失-曽失
	fprintf(ofp, "\n***掻失-曽失 朔亀 : 亜蟹陥 授***\n");
	for (i = sum = 0; i < 30*30; i++) {
		sum += jungjong[i];
		if (jungjong[i]) fprintf(ofp, "[%s%s]\t:%8d\t%6.2f\t%6.2f\n",
			ju[i/30], jo[i%30], jungjong[i],
			(double) jungjong[i] / nSyl * 100., (double) sum / nSyl * 100.);
	}

	set_sort_index(jungjong, si, 30*30);	// 壱朔亀 授生稽 sorting
	fprintf(ofp, "\n***掻失-曽失 朔亀 : 壱朔亀 授***\n");
	for (j = sum = 0; j < 30*30; j++) {
		i = si[j];
		sum += jungjong[i];
		if (jungjong[i]) fprintf(ofp, "[%s%s]\t:%8d\t%6.2f\t%6.2f\n",
			ju[i/30], jo[i%30], jungjong[i],
			(double) jungjong[i] / nSyl * 100., (double) sum / nSyl * 100.);
	}

// 曽失-段失
	fprintf(ofp, "\n恥 尻紗 製箭 鯵呪 : %d\n", nSyl2);
	fprintf(ofp, "\n***曽失-段失 朔亀 : 亜蟹陥 授***\n");
	for (i = sum = 0; i < 30*21; i++) {
		sum += jongcho[i];
		if (jongcho[i]) fprintf(ofp, "[%s%s]\t:%8d\t%6.2f\t%6.2f\n",
			jo[i/21], ch[i%21], jongcho[i],
			(double) jongcho[i] / nSyl2 * 100., (double) sum / nSyl2 * 100.);
	}

	set_sort_index(jongcho, si, 30*21);		// 壱朔亀 授生稽 sorting
	fprintf(ofp, "\n***曽失-段失 朔亀 : 壱朔亀 授***\n");
	for (j = sum = 0; j < 30*21; j++) {
		i = si[j];
		sum += jongcho[i];
		if (jongcho[i]) fprintf(ofp, "[%s%s]\t:%8d\t%6.2f\t%6.2f\n",
			jo[i/21], ch[i%21], jongcho[i],
			(double) jongcho[i] / nSyl2 * 100., (double) sum / nSyl2 * 100.);
	}

	if (ifp != stdin) fclose(ifp);
	if (ofp != stdout) fclose(ofp);
	return 0;
}