eclazz.cpp

来自「工程算法 这是一个很有用的工程数值算法集锦」· C++ 代码 · 共 34 行

int equiv(int i, int j)
{
    if ((i % 4) == (j % 4))
	{
		return -1;
	}
	else
	{
		return 0;
	}
}

void eclazz(int nf[], int n)
{
	int jj;
    nf[1] = 1;
    for (jj = 2; jj<=n; jj++)
	{
        nf[jj] = jj;
        for (int 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]];
    }
}