eclazz.c

来自「c数值算法的源码」· C语言 代码 · 共 16 行

void eclazz(nf,n,equiv)
int (*equiv)(),n,nf[];
{
	int kk,jj;

	nf[1]=1;
	for (jj=2;jj<=n;jj++) {
		nf[jj]=jj;
		for (kk=1;kk<=(jj-1);kk++) {
			nf[kk]=nf[nf[kk]];
			if ((*equiv)(jj,kk)) nf[nf[nf[kk]]]=jj;
		}
	}
	for (jj=1;jj<=n;jj++) nf[jj]=nf[nf[jj]];
}