Problem A: Sgu282 Isomorphism
Time Limit: 15 Sec Memory Limit: 64 MBSubmit: 172 Solved: 88
[Submit][Status][Discuss]
Description
给 定一个N 个结点的无向完全图( 任意两个结点之间有一条边), 现在你可以用 M 种颜色对这个图的每条边进行染色,每条边必须染一种颜色。 若两个已染色的图,其中一个图可以通过结点重新编号而与另一个图完全相同, 就称这两个染色方案相同。 现在问你有多少种本质不同的染色方法,输出结果 mod P。P 是一个大于N 的质数。
Input
仅一行包含三个数,N、M、P。
Output
仅一行,为染色方法数 mod P 的结果。
Sample Input
3 4 97
Sample Output
20
HINT
题解:
一眼可以看出这是置换群吧,但是它要求的是边置换,开始感觉没什么思路,但是想想一条边由(u,v)两个点构成
于是我们有了新的思路:考虑将点置换转换为边置换
我们可以发现点置换转化为的边置换同样具有相应的循环节
于是考虑使用polya定理解决这个问题
L:表示一个循环的大小; C:表示循环节的个数;
首先对于一条边(u,v)它要分为两种情况:
(1)u,v不在在同一个点循环,于是对于这条边所在的循环的大小为Lu-v=lcm(Lu,Lv),Cu-v=(Lu*Lv)/lcm(Lu,Lv)=gcd(Lu,Lv);
(2) u,v在同一个点循环,于是分奇数和偶数进行讨论:
一共C(L,2)条边 注:此处C为组合数
1.奇数:Li是奇数,每个循环覆盖Li条边;循环节个数:
2.偶数:Li是偶数,一个特例:覆盖Li/2条边的循环;循环节个数:
(3)由上可知循环节个数:
•实际N!个点置换中,有多少个 结构呢?
•一个循环看成一个圆排列,现在要把
N个人分配到
k个长度分别为
的独立不相关圆排列中
![](https://i-blog.csdnimg.cn/blog_migrate/eda3c843e80b1b587af7922042b1f98f.png)
•因为有 L
i==L
j
的情况,设Bi为有多少个Lj=i
•总分配数为
![](https://i-blog.csdnimg.cn/blog_migrate/a355b058627722d354fb51a2820e9810.png)
(4由上求出答案即可:
![](https://i-blog.csdnimg.cn/blog_migrate/8f900a89c6347c561fdf2122f13be562.gif)
![](https://i-blog.csdnimg.cn/blog_migrate/961ddebeb323a10fe0623af514929fc1.gif)
#include<iostream> #include<algorithm> #include<cstring> #include<cmath> #include<cstdio> #define ll long long #define N 70 using namespace std; ll n,m,p; ll num[N],l[N],bin[N]; ll ans; ll read() { ll x=0,f=1; char ch; while (ch=getchar(),ch<'0'||ch>'9') if (ch=='-') f=-1; while (x=x*10+ch-'0',ch=getchar(),ch>='0'&&ch<='9'); return x*f; } int gcd(int a,int b){return b?gcd(b,a%b):a; } ll ksm(int x,int k){ll res=1; for (int i=k; i; i>>=1,x=1ll*x*x%p) if (i&1) res=1ll*res*x%p; return res; } void cal(int k) { ll s=1; for (int i=1; i<=n; i++) num[i]=0; for (int i=1; i<=k ;i++) num[l[i]]++; for (int i=1; i<=k; i++) s=s*l[i]%p; for (int i=1; i<=n; i++) s=s*bin[num[i]]%p; int c=0; for (int i=1; i<=k; i++) { c+=l[i]/2; for (int j=i+1; j<=k; j++) c+=gcd(l[i],l[j]); } ans=(ans+ksm(m,c)*bin[n]%p*ksm(s,p-2)%p)%p; } void dfs(int x,int k,int last) { if (x==n+1) {cal(k-1); return;} for (int i=1; i<=last && x+i<=n+1; i++) { l[k]=i; dfs(x+i,k+1,i); } } int main() { n=read(); m=read(); p=read(); bin[0]=1; for (int i=1; i<=n; i++) bin[i]=bin[i-1]*i%p; dfs(1,1,n); printf("%lld\n",ans*ksm(bin[n],p-2)%p); return 0; }