P4128 [SHOI2006]有色图
数学渣渣看题解看得想死Ծ‸Ծ
首先发现这玩意儿看着很像polya定理
\]
然而polya定理只能用来求点的置换,边的置换是布星的
于是我们考虑一个点的置换,把它写成若干循环的乘积\((a_1,a_2,..)(b_1,b_2,...)...\)
1.对于不在同一个循环里的点,比方说一条边\((a_1,b_1)\),那么和它在同一个循环的边有\(((a_1,b_1),(a_2,b_2),...)\)设\(a\)的循环节为\(l_1\),\(b\)的循环节为\(l_2\),那么这个边的循环的循环节长度就是\(lcm(l_1,l_2)\),而总共的边数为\(l_1*l_2\),那么循环的个数就是\(\frac{l_1*l_2}{lcm(l_1,l_2)}=gcd(l_1,l_2)\)
2.对于在同一个循环内的点,设\(a\)的循环节长度为\(l_1\)
如果\(l_1\)长度为奇数,那么循环的长度就是\(l_1\),总共有\(C_{l_1}^{2}\)条边,那么循环的个数就是\(\frac{l-1}2\)
如果\(l_1\)长度为偶数,除了上面的情况之外,还有一种很gg的情况就是比方说\((a_1,a_{l_1/2+1})\)的边所在的循环,这个循环的长度是\(l_1/2\),占的边数为\(l_1/2\),除此之外其他的情况都是一样的,所以循环的个数就是\(\frac{\frac{l(l-1)}{2}-\frac{l}2}{l}+1=\frac l 2\)
那么,总共的边的置换就是
\]
然而如果直接枚举的话复杂度是带一个感叹号的……然而我们只需要知道所有的\(l_i\)就行了,而不需要知道循环里具体是什么数字……所以会T就是我们知道的太多了
于是我们可以枚举\(l_i\),为了不重不漏保证\(l_i\)不降,先考虑如果\(l_i\)互不相同的话有多少种方案。我们可以这样理解,枚举\(n\)个数字的全排列,然后按\(l_i\)从左到右依次分组,那么这样肯定就是一个置换了……然后对于其中的每一个循环\((a_1,a_2,...,a_l)\)来说,\((a_2,a_3,...,a_l,a_1)...\)之类的其实是跟它一样的,也就是说每个循环有\(l_i\)个同构的,所以要除掉,那么方案数就是
\]
然而现在问题是\(l_i\)有可能会相等,如果按上面那样考虑的话有可能会有两个\(l_i\)相等的循环被算到不同的里面了……所以还要设\(B_i\)为\(l_j==i\)的个数,然后除掉他们中间排列的个数,那么方案数应该是
\]
然后dfs暴力找\(l_i\)即可
最后就是
\]
之后就是抄代码了
//minamoto
#include<bits/stdc++.h>
#define int long long
#define fp(i,a,b) for(register int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(register int i=a,I=b-1;i>I;--i)
using namespace std;
const int N=105;
int ans,P,n,m,fac[N],rec[N];
int gcd(int x,int y){return y?gcd(y,x%y):x;}
int ksm(int x,int y){
int res=1;
for(;y;y>>=1,x=x*x%P)if(y&1)res=res*x%P;
return res;
}
void calc(int x){
int sum=0,mul=1,now=1;
fp(i,1,x)sum+=rec[i]/2;
fp(i,1,x)fp(j,i+1,x)sum+=gcd(rec[i],rec[j]);
fp(i,1,x)(mul*=rec[i])%=P;
fp(i,2,x){
if(rec[i]!=rec[i-1])(mul*=fac[now])%=P,now=0;
++now;
}(mul*=fac[now])%=P,mul=fac[n]*ksm(mul,P-2)%P;
(ans+=mul*ksm(m,sum)%P)%=P;
}
void dfs(int k,int x,int s){
if(!x)calc(k-1);if(x<s)return;
fp(i,s,x)rec[k]=i,dfs(k+1,x-i,i);
}
signed main(){
// freopen("testdata.in","r",stdin);
scanf("%lld%lld%lld",&n,&m,&P),fac[0]=1;
fp(i,1,n)fac[i]=fac[i-1]*i%P;
dfs(1,n,1);(ans*=ksm(fac[n],P-2))%=P;
printf("%lld\n",ans);return 0;
}
P4128 [SHOI2006]有色图的更多相关文章
- 洛谷 P4128 [SHOI2006]有色图 解题报告
P4128 [SHOI2006]有色图 题目描述 如果一张无向完全图(完全图就是任意两个不同的顶点之间有且仅有一条边相连)的每条边都被染成了一种颜色,我们就称这种图为有色图.如果两张有色图有相同数量的 ...
- 洛谷 P4128: bzoj 1815: [SHOI2006]有色图
题目传送门:洛谷 P4128. 计数好题,原来是 13 年前就出现了经典套路啊.这题在当年应该很难吧. 题意简述: \(n\) 个点的完全图,点没有颜色,边有 \(m\) 种颜色,问本质不同的图的数量 ...
- [SHOI2006] 有色图
Description 给一张 \(n\) 个点的无向完全图,同时还有 \(m\) 种颜色.要求给每条边染色,问有多少种不同的染色方案.两种方案不同当且仅当顶点标号任意重排后不同.\(n\leq 53 ...
- BZOJ1815 SHOI2006有色图(Polya定理)
置换数量是阶乘级别的,但容易发现本质不同的点的置换数量仅仅是n的整数拆分个数,OEIS(或者写个dp或者暴力)一下会发现不是很大,当n=53时约在3e5左右. 于是暴力枚举点的置换,并且发现根据点的置 ...
- bzoj 1815: [Shoi2006]color 有色图 置换群
1815: [Shoi2006]color 有色图 Time Limit: 4 Sec Memory Limit: 64 MBSubmit: 136 Solved: 50[Submit][Stat ...
- BZOJ1815: [Shoi2006]color 有色图
BZOJ1815: [Shoi2006]color 有色图 Description Input 输入三个整数N,M,P 1< = N <= 53 1< = M < = 1000 ...
- BZOJ 1815: [Shoi2006]color 有色图(Polya定理)
题意 如果一张无向完全图(完全图就是任意两个不同的顶点之间有且仅有一条边相连)的每条边都被染成了一种颜色,我们就称这种图为有色图. 如果两张有色图有相同数量的顶点,而且经过某种顶点编号的重排,能够使得 ...
- BZOJ 1815: [Shoi2006]color 有色图 [Polya DFS 重复合并]
传送门 题意: 染色图是无向完全图,且每条边可被染成k种颜色中的一种.两个染色图是同构的,当且仅当可以改变一个图的顶点的编号,使得两个染色图完全相同.问N个顶点,k种颜色,本质不同的染色图个数(模质数 ...
- [SHOI2006]color 有色图[群论、组合计数]
题意 用 \(m\) 种颜色,给 \(n\) 个点的无向完全图的 \(\frac{n(n-1)}{2}\) 条边染色,两种方案相同当且仅当一种方案交换一些点的编号后可以变成另一种方案.问有多少本质不同 ...
随机推荐
- css伪类实现文字两侧划线效果
css伪类实现文字两侧划线效果,效果如下: 代码如下: <!DOCTYPE HTML> <html> <head> <title> css伪类的学习 & ...
- POJ-1274The Perfect Stall,二分匹配裸模板题
The Perfect Stall Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 23313 Accepted: 103 ...
- Bone Collector II(hdu 2639)
题意:求01背包的第k最优值 输入:第一行为T,下面是T组数据,每组数据有n,m,k 代表n件物品,m容量,和题目要求的k,下一行是n个物品的价值,再一行是n个物品的体积 输出:T行答案 /* 类似于 ...
- 主席树初探--BZOJ2588: Spoj 10628. Count on a tree
n<=100000的点权树,有m<=100000个询问,每次问两个点间的第k小点权,保证有解,强制在线. 主席上树啦!类似于之前的序列不带修改询问的前缀表示法,现在只要把前缀当成某点到根的 ...
- linux 命令练习 2018-08-27
linux 命令练习 2018-08-27 uname 显示系统名字 [test@localhost ~]$ uname Linux uname -a 即列出linux的内核版本号 [test@ ...
- springmvc json 数据
这里是controllor层 @RequestMapping("/traceupdatestatus") @ResponseBody public boolean traceupd ...
- linux下crontab安装和使用(定时任务)
在Unix和Linux的操作系统之中,cron可以让系统在指定的时间,去执行某个指定的任务,crontab命令常用于管理对应的cron. 一.crontab在线安装 yum -y install vi ...
- list去掉重复元素
需求: 有list 里面含有重复元素,要求去掉重复元素: solution 1: >>> a [1, 2, 2, 1, 3, 4, 5, 6, 5] >>> set ...
- spark开发环境配置
以后spark,mapreduce,mpi可能三者集于同一平台,各自的侧重点有所不用,相当于云计算与高性能计算的集合,互补,把spark的基础看了看,现在把开发环境看看,主要是看源码,最近Apache ...
- Nginx/Spring:增加上传文件尺寸限制
Nginx: 1. vi /etc/nginx/nginx.config 2. 在server中增加如下行 # set client body size to 20M # client_max_bod ...