xeclass.c

Numerical Recipes in C的源代码

/* Driver for routine eclass */

#include <stdio.h>
#include "nr.h"
#include "nrutil.h"

#define M 11
#define N 15

main()
{
	int i,j,k,lclas=0,nclass,*nf,*nflag,*nsav;
	static int lista[]={0,1,1,5,2,6,2,7,11,3,4,12},
		   listb[]={0,5,9,13,6,10,14,3,7,15,8,4};

	nf=ivector(1,N);
	nflag=ivector(1,N);
	nsav=ivector(1,N);
	eclass(nf,N,lista,listb,M);
	for (i=1;i<=N;i++) nflag[i]=1;
	printf("\nNumbers from 1-%d divided according to\n",N);
	printf("their value modulo 4:\n\n");
	for (i=1;i<=N;i++) {
		nclass=nf[i];
		if (nflag[nclass]) {
			nflag[nclass]=0;
			lclas++;
			k=0;
			for (j=i;j<=N;j++)
				if (nf[j] == nf[i]) nsav[++k]=j;
			printf("Class %2d:      ",lclas);
			for (j=1;j<=k;j++) printf("%3d",nsav[j]);
			printf("\n");
		}
	}
	free_ivector(nsav,1,N);
	free_ivector(nflag,1,N);
	free_ivector(nf,1,N);
	return 0;
}