UVA11255 Necklace

http://www.elijahqi.win/archives/3401
题目：https://uva.onlinejudge.org/external/112/p11255.pdf

每种颜色有了个数限制怎么办 首先还是按照polya一样做 只不过这次 我们要求有颜色数限制

那么考虑旋转同构的情况 考虑枚举长度 那么要求每个那么要求每个颜色的数量必须是枚举的长度的一个倍数 因为最后我所有颜色都要用完 那么现在问题变成了我有k个循环节 每个循环节只能填一种颜色求方案数 那么就是C_{sum}^{a}*C_{sum-a}^{b}*C_{sum-a-b}^{c}

即我需要分配一下每个循环分别要用哪个颜色来填

然后计算方案之后因为这样一共有phi个再*phi即可

考虑翻转同构 奇数： 每次枚举我钦定的点是什么颜色 然后这个颜色减1 做一遍*n 因为有n种置换

偶数：每次枚举钦定的两个点的颜色 然后分别-1 做一遍*n/2

然后再做一下翻转是以中间空隙作为对称轴的情况

#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if(T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
const int N=44;
ll c[N][N];bool not_prime[N];
int prime[N],tot,a[3],b[3],T,phi[N],n;
inline void init(){
    for (int i=0;i<=40;++i) c[i][0]=1;phi[1]=1;
    for (int i=1;i<=40;++i)
        for (int j=1;j<=i;++j) c[i][j]=c[i-1][j-1]+c[i-1][j];
    for (int i=2;i<=40;++i){
        if (!not_prime[i]) prime[++tot]=i,phi[i]=i-1;
        for (int j=1;prime[j]*i<=40;++j){
            not_prime[prime[j]*i]=1;
            if (i%prime[j]==0){
                phi[prime[j]*i]=phi[i]*prime[j];break;
            }else phi[prime[j]*i]=phi[i]*phi[prime[j]];
        }
    }
}
inline void mv(){memcpy(b,a,sizeof(a));}
inline ll calc(int len){
    ll sum=0,tmp=1;
    for (int i=0;i<3;++i){
        if (b[i]%len) return 0;b[i]/=len;sum+=b[i];
    }
    for (int i=0;i<3;++i)
        tmp=tmp*c[sum][b[i]],sum-=b[i];return tmp;
}
int main(){
    freopen("uva11255.in","r",stdin);
    T=read();init();
    while(T--){n=0;
        for (int i=0;i<3;++i) a[i]=read(),n+=a[i];ll ans=0;
        for (int i=1;i*i<=n;++i){
            if (n%i) continue;
            mv();ans+=calc(i)*phi[i];
            if (i*i==n)continue;
            mv();ans+=calc(n/i)*phi[n/i];
        }
        if (n&1){
            for (int i=0;i<3;++i){
                if (a[i]&1==0) continue;
                mv();--b[i];ans+=calc(2)*n;
            }
        }else{
            for (int i=0;i<3;++i){
                for (int j=0;j<3;++j){
                    mv();--b[i];--b[j];ans+=n/2*calc(2);
                }
            }mv();ans+=n/2*calc(2);
        }
        ans/=n<<1;
        printf("%lld\n",ans);
    }
    return 0;
}