【算法讲27：Polya定理】BurnSide 引理 Polya定理 | 洛谷 P4980

最新推荐文章于 2021-08-25 20:21:48 发布

溢流眼泪

最新推荐文章于 2021-08-25 20:21:48 发布

阅读量518

点赞数 1

分类专栏：【算法/知识点浅谈】文章标签：数论 polya定理

本文链接：https://blog.csdn.net/weixin_45775438/article/details/118945017

版权

【算法/知识点浅谈】专栏收录该内容

65 篇文章 10 订阅

订阅专栏

【算法讲27：Polya定理】洛谷 P4980

前置
引入
前置： $B u r n S i d e$ 引理
Polya 定理
回到题目
代码

注：下文证明都略，只记录重要的性质结论

前置

群，置换
下文的记号：
$g$ 置换
$G$ 置换集
$∣ G ∣$ 置换集大小
$(a,b,c,\cdots)$ 一个置换

引入

洛谷 P4980
一个串，上面有 $N$ 个珠子，用 $N$ 种颜色给它染色
求有多少个本质不同的染色方案
本质相同：可以通过旋转，让两串珠子相同
样例数 $1\le T\le 1000$
$1\le N\le 10^9$

前置： $B u r n S i d e$ 引理

不过貌似讲讲都挺麻烦的…
把祖传图给薅过来
问这样的 $2\times2$ 的矩形，每个块染色 $0 / 1$ ，有多少种本质不同的染色
本质相同的染色：通过旋转得到相同的块算本质相同
虽然貌似有 $16$ 种染色方案，但是我们通过如下的操作，某种方案会转换成某种方案：
$\begin{aligned} &顺时针旋转0\degree：&(1)(2)\cdots(16)\\ &逆时针旋转90\degree：&(1)(2)(3,4,5,6)(7,8,9,10)(11,12)(13,14,15,16)\\ &顺时针旋转90\degree：&(1)(2)(6,5,4,3)(10,9,8,7)(12,11)(16,15,14,13)\\ &旋转180\degree：&(1)(2)(3,5)(4,6)(7,9)(8,10)(11)(12)(13,15)(14,16)\\ \end{aligned}$
其中 $(a_1,a_2,\cdots,a_k)$ 表示一个置换群
$B u r n S i d e$ 引理：
$|X/G|=\frac{1}{|G|}\sum_{g\in G}|X^g|$
正规表述：轨道数 $=$ 被 $G$ 中一个元素保持不动的点的个数的平均值
见：百度百科
在这里 $2\times2$ 矩阵染色问题中， $|X/G|=\frac{16+2+2+4}{4}=6$ ，故只有 $6$ 种合法染色方案

Polya 定理

由于我们太难把所有的置换等价操作给求出来（太多了…）故有更优化的 $P o l y a$ 定理
它的关键就是给出计算不动点的个数的方法
$P o l y a$ 定理：
$Ans=\frac{1}{|G|}\sum_{g\in G}m^{c(g)}$
表述：对于有 $m$ 种颜色染色的问题， $X^g|=m^{c(g)}$
其中 $c (g)$ 表示在置换 $g$ 的情况下，不动点的个数
在这里 $2\times2$ 矩阵染色问题中，查看四个格子的置换关系：
$\begin{aligned} &顺时针旋转0\degree：&(1)(2)(3)(4)\\ &顺时针旋转90\degree：&(4,3,2,1)\\ &顺时针旋转180\degree：&(1,3)(2,4)\\ &顺时针旋转270\degree：&(1,2,3,4)\\ \end{aligned}$
于是 $Ans=\frac{1}{4}\times(2^4+2^1+2^1+2^2)=6$ ，答案也是一样的

回到题目

我们直接使用 $P o l y a$ 定理，对于本质相同的操作，只有旋转操作
旋转操作有 $n$ 种，就是顺时针旋转 $(0\sim n-1)$ 个珠子，当然也可以写成旋转 $(1\sim n)$ 个珠子
考虑，如果我们顺时针旋转 $k$ 个珠子，会获得多少个不动点
根据我们的数学知识，答案为 $\gcd(k,n)$
那么答案就是
$Ans=\frac{1}{|n|}\sum_{k=1}^n n^{\gcd(k,n)}$
考虑到 $n$ 非常大，我们不能直接枚举，我们像做莫比乌斯反演一样，套路化做：
$\begin{aligned} Ans&=\frac{1}{n}\sum_{d|n}n^d\sum_{i=1}^n\Big[ \gcd(k,n)=d \Big]\\ &=\frac{1}{n}\sum_{d|n}n^d\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\Big[ \gcd(k,\frac{n}{d})=1 \Big]\\ &=\frac{1}{n}\sum_{d|n}n^d\times\varphi(\frac{n}{d}) \end{aligned}$
暴力枚举，时间复杂度为 $O(T\times n^{3/4})$ ，差不多也能过

代码

#include <bits/stdc++.h>
#define IOS ios::sync_with_stdio(false);cin.tie(NULL);cout.tie(NULL);
using namespace std;
typedef long long ll;
void show(){std::cerr << endl;}template<typename T,typename... Args>void show(T x,Args... args){std::cerr << "[ " << x <<  " ] , ";show(args...);}

const int MAX = 2e5+50;
const int MOD = 1e9+7;
const int INF = 0x3f3f3f3f;
const ll LINF = 0x3f3f3f3f3f3f3f3f;
const double EPS = 1e-5;

ll qpow(ll a,ll n){/* */ll res = 1LL;while(n){if(n&1)res=res*a%MOD;a=a*a%MOD;n>>=1;}return res;}
ll qpow(ll a,ll n,ll p){a%=p;ll res = 1LL;while(n){if(n&1)res=res*a%p;a=a*a%p;n>>=1;}return res;}
ll npow(ll a,ll n){/* */ll res = 1LL;while(n){if(n&1)res=res*a;a=a*a;n>>=1;if(res<0||a<0)return 0;}return res;}
ll inv(ll a){/* */return qpow(a,MOD-2);}
ll inv(ll a,ll p){return qpow(a,p-2,p);}


ll phi(ll n){
    ll res = n;
    for(int i = 2;i * i <= n;++i){
        if(n % i == 0){
            res = res * (1 - inv(i) + MOD) % MOD;
            while(n % i == 0)n /= i;
        }
    }
    if(n > 1)res = res * (1 - inv(n) + MOD) % MOD;
    return (res + MOD) % MOD;
}


int main()
{
    IOS;
    int T;cin >> T;
    while(T--){
        int n;cin >> n;
        ll res = 0;
        for(int i = 1;i * i <= n;++i){
            if(n % i == 0){
                res = (res + qpow(n,i) * phi(n / i) % MOD) % MOD;
                if(i * i != n)res = (res + qpow(n,n/i) * phi(i) % MOD) % MOD;
            }
        }
        res = res * inv(n) % MOD;
        cout << res << endl;
    }
    return 0;
}