时钟 51Nod - 1282

https://www.51nod.com/onlineJudge/questionCode.html#!problemId=1282

每个钟表有m个状态 每个状态hash为一个三进制数 只用最小的一个来表示该状态即可 因为两个钟表同构的话 每个状态hash出的值都是一样的 最后并查集求出答案即可 复杂度n*n*logn

#include <bits/stdc++.h>
using namespace std;
typedef unsigned long long ull;
const ull N=1000000000;

ull val[510];
int per[510][510];
int book[510],f[510],cnt[510];
int n,m,p;

int getf(int p)
{
    if(f[p]==p) return p;
    else
    {
        f[p]=getf(f[p]);
        return f[p];
    }
}

void unite(int u,int v)
{
    int fu,fv;
    fu=getf(u),fv=getf(v);
    if(fu!=fv) f[fv]=fu,cnt[fu]+=cnt[fv];
}

ull quickpow(ull a,int b)
{
    ull res;
    res=1;
    while(b>0)
    {
        if(b%2) res=(res*a);
        a*=a,b/=2;
    }
    return res;
}

int main()
{
    ull gou,res,tmp;
    int i,j,k,l,ans;
    scanf("%d%d%d",&n,&m,&p);
    tmp=quickpow(3,p-1);
    //printf("***%llu***\n",tmp);
    for(i=0;i<n;i++)
    {
        for(j=0;j<m;j++)
        {
            scanf("%d",&per[i][j]);
            per[i][j]%=p;
        }
        sort(per[i],per[i]+m);
        for(j=m-1;j>=0;j--) per[i][j]-=per[i][0];
        gou=0;
        for(j=0;j<m;j++)
        {
            if(j!=0) res=quickpow(3,per[i][j]-per[i][j-1]);
            else res=1;
            gou=res*gou+1;
        }
        if(per[i][m-1]<p-1) gou*=quickpow(3,p-per[i][m-1]-1);

        /*
        printf("***%llu***\n",gou);
        for(j=0;j<m;j++) printf("%d ",per[i][j]);
        printf("\n\n");
        */

        val[i]=gou;
        for(j=1;j<m;j++)
        {
            gou=(gou-tmp)*quickpow(3,per[i][j]-per[i][j-1])+quickpow(3,per[i][j]-per[i][j-1]-1);
            val[i]=min(val[i],gou);
        }
    }

    //for(i=0;i<n;i++) printf("%llu ",val[i]);
    //printf("\n");

    for(i=0;i<n;i++) f[i]=i,cnt[i]=1;
    for(i=0;i<n;i++)
    {
        for(j=i+1;j<n;j++)
        {
            if(val[i]==val[j]) unite(i,j);
        }
    }
    ans=0;
    for(i=0;i<n;i++)
    {
        if(f[i]==i) ans+=(cnt[i]*(cnt[i]-1))/2;
    }
    printf("%d\n",ans);
    return 0;
}