bzoj4330

最新推荐文章于 2017-04-23 20:42:42 发布

la1la1la_

最新推荐文章于 2017-04-23 20:42:42 发布

阅读量400

点赞数

分类专栏： polya bzoj

本文链接：https://blog.csdn.net/u010600261/article/details/53001971

版权

bzoj 同时被 2 个专栏收录

24 篇文章 0 订阅

订阅专栏

polya

2 篇文章 0 订阅

订阅专栏

不会polya。。用的是burnside引理。
假设有n种置换法，本质不同方案就等于在每一种置换下不变的情况数之和除n
即： $\frac{1}{n}\sum_{i=1置换pi}^n有多少种方案在pi变换前和变换后一样$
题意：
n个戒指串成一个项链，每个项链都是一个有m个宝石的环，宝石共有r种颜色，问有多少种不同方案。要求项链上相邻戒指不同。戒指在旋转下相同视为相同，项链不用考虑旋转（即长得不一样就算两种方案）。N<=10^15，M<=10^9，R<=10^6。

#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<iostream>
#define mmod 3214567
#define N 1100000
#define LL long long
using namespace std;
struct mat{
    LL a[2][2];
    mat(){memset(a,0,sizeof(a));}
    friend mat operator *(mat a,mat b)
    {
        mat c;
        for(LL i=0;i<2;i++)
            for(LL j=0;j<2;j++)
                for(LL k=0;k<2;k++)
                    c.a[i][j]=(c.a[i][j]+a.a[i][k]*b.a[k][j])%mmod;
        return c;
    }
}a,g;
LL n,m,c[N],ans;
LL gcd(LL a,LL b)
{
    if(b==0) return a;
    return gcd(b,a%b);
}
LL qmod(LL a,LL b)
{
    LL t=1;
    while(b)
    {
        if(b%2) t=t*a%mmod;
        b/=2;
        a=a*a%mmod;
    }
    return t;
}
LL phi(LL x)
{
    int q=sqrt(x),res=x;
    for(int i=2;i<=q;i++)
    {
        if(x%i) continue;
        res=res/i*(i-1);
        while(x%i==0) x/=i;
        if(i*i>x) break;
    }
    if(x>1) res=res/x*(x-1);
    return res;
}
void init()
{
    LL r,t=0;
    scanf("%lld%lld%lld",&n,&m,&r);
    int q=sqrt(m);
    for(LL i=1;i<=q;i++)
    {
        if(m%i) continue;
        LL g=i;
        t=(t+qmod(r,g)*phi(m/g))%mmod;
        if(i*i==m) break;
        g=m/i;
        t=(t+qmod(r,g)*phi(m/g))%mmod;
    }
    LL ny=qmod(m,mmod-2);
    t=t*ny%mmod;
    m=t;
}
void solve()
{
    a.a[0][0]=m-2;a.a[0][1]=1;
    a.a[1][0]=m-1;a.a[1][1]=0;
    for(LL i=0;i<=1;i++) g.a[i][i]=1;
    n--;
    while(n)
    {
        if(n%2) g=g*a;
        n/=2;
        a=a*a;
    }
    ans=m*g.a[1][0]%mmod;
    ans=(ans+mmod)%mmod;
    printf("%lld\n",ans);
}
int main()
{
    init();
    solve();
    return 0;
}

题解：
先求出有多少种不同的戒指。因为戒指可以旋转，所以要用polya去重。在这里置换就是顺时针转1到m个单位。当我们转k个单位的时候，相当于i向(i+k)%n连边。这样每个点都是一入度一出度，形成了若干个环。要使旋转前后一样，显然同一个环要染同种颜色。若有g个环，那方案数就是 $r^g$ 。g=gcd(m,k)。证明就是我们想象一个点x，如果x要从环上绕回x，那走的步数为n的倍数，每次加k，第一次碰到m的倍数就是lcm(m,k)，步数即为 $\frac{lcm(m,k)}{k}$ 。也就是说有 $\frac{m}{\frac{lcm(m,k)}{k}}=g$ 个环。显然枚举k会超时，但注意我们关注的只是环的数量，那么就枚举gcd(k,m)=g，能和m取到gcd为g的k显然有 $\phi(\frac{m}{g})$ 个。
然后设f[i][0,1]表示放了前i个戒指，第i个戒指和第一个是否相同方案数，矩乘即可。