poj 2888 Magic Bracelet <polya定理>
题目:http://poj.org/problem?id=2888
题意:给定n(n <= 10^9)颗珠子,组成一串项链,每颗珠子可以用m种颜色中一种来涂色,如果两种涂色方法通过旋转项链可以得到视为等价。
然后再给定K组限制,每组限制a、b代表颜色a和颜色b不能涂在相邻的珠子上面。问一共有多少种涂色方法。
思路: 如果这题没有后面的限制,就和 poj 2154 一样了:http://www.cnblogs.com/jian1573/p/3234627.html
现在我们要处理的就是 K 种限制, 可以用DP求解。 i为珠子编号, c为颜色编号那么:dp[i][c]=∑dp[i-1][cc] cc 为可以与 c 相邻的颜色编号;
由于N为 1e9 O(N) 会TLE, 所以我们可以用矩阵快速幂来优化为O(lgN):具体用m[i][j]=1,表示合法, m[i][j]=0表示不合法,
那么m^k后的对角线上的元素和即为所求。
#include <iostream>
#include <cmath>
#include <cstdio>
#include <cstring>
using namespace std;
const int Mod=;
int N, M, K, T;
struct Mar
{
int m[][];
inline void zero( ){
memset(m, , sizeof m);
}
inline void one( ){
zero( );
for( int i=; i<M; ++i ){
for( int j=; j<M; ++j ){
m[i][j]=;
}
}
}
inline void unit( ){
zero();
for( int i=; i<M; ++ i )
m[i][i]=;
}
inline Mar operator *(const Mar &a) const {
Mar C;C.zero();
for( int i=; i<M; ++i ){
for(int j=; j<M; ++j ){
for( int k=; k<M; ++k ){
C.m[i][j]+=m[i][k]*a.m[k][j];
C.m[i][j]%=Mod;
}
}
}
return C;
}
inline Mar operator ^ (int t) const{
Mar B=*this, C;
C.unit( );
while(t){
if(t&)C=C*B;
B=B*B;
t>>=;
}
return C;
}
}A;
int a[], p[],cntp=;
void getp( )
{
for( int i=; i<=1e5; ++ i ){
if( !a[i] )p[cntp++]=i;
for( int j=; j<cntp&&i*p[j]<1e5; ++ j ){
a[i*p[j]]=;
if( i%p[j]== )break;
}
}
} int Phi( int x )
{
int res=x;
for( int i=; i<cntp&&p[i]*p[i]<=x; ++i ){
if( x%p[i]== ){
res/=p[i]; res*=(p[i]-);
while(x%p[i]==){
x/=p[i];
}
}
}
if(x>){
res/=x;res*=(x-);
}
return res%Mod;
}
int P_M(int a, int b)
{
int res=;
while(b){
if(b&)res*=a, res%=Mod;
a*=a, a%=Mod;
b>>=;
}
return res; }
int work( int k )
{
Mar C=A^k;
int res=;
for( int i=; i<M; ++i )
res+=C.m[i][i];
return res;
}
int polya( )
{
int ans=;
for( int i=; i*i<=N; ++ i ){
if( N%i== ){
if(i*i==N){
ans+=Phi(i)*work(i);
ans%=Mod;
}
else{
ans+=Phi(i)*work(N/i);
ans%=Mod;
ans+=Phi(N/i)*work(i);
ans%=Mod;
}
}
}
return ans;
}
int main( )
{
getp();
scanf("%d", &T);
while (T--){
scanf("%d%d%d", &N, &M, &K);
A.one();
for( int i=, x, y; i<K; ++i ){
scanf("%d%d", &x, &y);
A.m[x-][y-]=A.m[y-][x-]=;
}
int ans=polya();
int inv=P_M(N%Mod, Mod-);
ans*=inv;ans%=Mod;
printf("%d\n", ans);
}
return ;
}
poj 2888 Magic Bracelet <polya定理>的更多相关文章
- poj 2888 Magic Bracelet(Polya+矩阵快速幂)
Magic Bracelet Time Limit: 2000MS Memory Limit: 131072K Total Submissions: 4990 Accepted: 1610 D ...
- [POJ 2888]Magic Bracelet[Polya Burnside 置换 矩阵]
也许更好的阅读体验 \(\mathcal{Description}\) 大意:给一条长度为\(n\)的项链,有\(m\)种颜色,另有\(k\)条限制,每条限制为不允许\(x,y\)颜色连在一起.要求有 ...
- POJ 2888 Magic Bracelet [Polya 矩阵乘法]
传送门 题意:竟然扯到哈利波特了.... 和上一题差不多,但颜色数很少,给出不能相邻的颜色对 可以相邻的连边建图矩阵乘法求回路个数就得到$f(i)$了.... 感觉这样的环上有限制问题挺套路的...旋 ...
- POJ 2888 Magic Bracelet(Burnside引理,矩阵优化)
Magic Bracelet Time Limit: 2000MS Memory Limit: 131072K Total Submissions: 3731 Accepted: 1227 D ...
- 解题:POJ 2888 Magic Bracelet
题面 这题虽然很老了但是挺好的 仍然套Burnside引理(因为有限制你并不能套Polya定理),思路和这个题一样,问题主要是如何求方案. 思路是把放珠子的方案看成一张图,然后就巧妙的变成了一个经典的 ...
- poj 2888 Magic Bracelet
经典的有限制条件的Burnside计数+矩阵乘法!!! 对于这种限制条件的情况我们可以通过矩阵连乘得到,先初始化矩阵array[i][j]为1.如果颜色a和颜色b不能涂在相邻的珠子, 那么array[ ...
- POJ 2888 Magic Bracelet(burnside引理+矩阵)
题意:一个长度为n的项链,m种颜色染色每个珠子.一些限制给出有些颜色珠子不能相邻.旋转后相同视为相同.有多少种不同的项链? 思路:这题有点综合,首先,我们对于每个n的因数i,都考虑这个因数i下的不变置 ...
- POJ 2888 Magic Bracelet ——Burnside引理
[题目分析] 同样是Burnside引理.但是有几种颜色是不能放在一起的. 所以DP就好了. 然后T掉 所以矩阵乘法就好了. 然后T掉 所以取模取的少一些,矩阵乘法里的取模尤其要注意,就可以了. A掉 ...
- poj 2154 Color【polya定理+欧拉函数】
根据polya定理,答案应该是 \[ \frac{1}{n}\sum_{i=1}^{n}n^{gcd(i,n)} \] 但是这个显然不能直接求,因为n是1e9级别的,所以推一波式子: \[ \frac ...
随机推荐
- 【Vijos1222】等值拉面(DP)
题意:有N个数对(a[i],b[i]) 每次可以把(x,y)变成(x+a[i],y+b[i])或(x+b[i],x+a[i]),后者称为交换一次 求使abs(x-y)最小时的最小交换次数 n<= ...
- 更改 terminal 開啟時,預設的路徑
echo "cd /media" >> ~/.bashrc open a new terminal (ctrl + alt + t)
- jq ajax之beforesend(XHR)
用于在向服务器发送请求前添加一些处理函数.例如:搜寻附近店铺,未搜寻出来的时候显示:正在搜寻中... 常见的一种效果,在用ajax请求时,没有返回前会出现前出现一个转动的loading小图标或者“内容 ...
- loading上下左右居中
.loading { width: 50px; height: 60px; text-align: center; font-size: 10px; position:fixed; left:50%; ...
- 2017CCPC 哈尔滨 B
这题没有考虑到m这个东西,所以就没有往二分答案的方向想 二分答案 check的时候,我们找的是大于等于x的数有多少个被加入到那个数组中.如果 >= m说明这个数可能是答案,否则就不是. 用尺取来 ...
- Network | sk_buff
sk_buff结构可能是linux网络代码中最重要的数据结构,它表示接收或发送数据包的包头信息.它在中定义,并包含很多成员变量供网络代码中的各子系统使用. 这个结构被不同的网络层(MAC或者其他 ...
- python中执行shell命令的几个方法
1.os.system() a=os.system("df -hT | awk 'NR==3{print $(NF-1)}'") 该命令会在页面上打印输出结果,但变量不会保留结果, ...
- Codeforces 235 E Number Challenge
Discription Let's denote d(n) as the number of divisors of a positive integer n. You are given three ...
- down
Description 给出一个数列,求出这个序列的最长下降子序列的长度及方案数,子序列中的权值完全相同视为同一个序列 Input 第一行一个整数n,接下来一行n个整数表示序列的权值 Output 一 ...
- 官方Java编码规范
先由Sun制定,之后Sun把Java卖给了Oracle,最后就成了Oracle制定的了.但是版本比较旧了,停留在1999年. 相比Google的编码规范,罪名写的却别就是Sun采用的是4个空格进行缩进 ...