polya/burnside 学习
参考链接:
http://www.cnblogs.com/hankers/archive/2012/08/03/2622231.html
http://blog.csdn.net/raalghul/article/details/51767941
首先来说说burnside引理是什么。
一天你正在刷题,看到一道关于染色的问题,你认为是一个傻逼题,然后认真一看题目上面写着旋转、翻转后相同的计算一次......你立刻就傻眼了。
接下来是科普时间。
首先我们考虑什么东西叫置换,例如(1,2,3,4,5)->(2,1,4,5,3)就是一个置换,1喂进去会变成2,3喂进去会变成4,喂进去一个(5,4,3,2,1)会变成(3,5,4,1,2)。
老司机们可能会把它写成这样:
$\begin{pmatrix}1&2&3&4&5\\2&1&4&5&3\end{pmatrix}\\$
例如现在三角形有三个顶点123,那么旋转本质上就是置换:(1,2,3)->(2,3,1)和(1,2,3)->(3,1,2)。
那么burnside引理就是对于k个置换,设在每个置换的作用下保持不变的染色方案为z1...zk,那么在置换意义下不同的染色方案一共有(z1+z2+...+zk)/k种。
应用的时候需要注意,例如转一次和转两次这两个置换都需要考虑,以及本身也是一个置换。
举个老例子,4个格子排成2*2的正方形,用两种颜色染色,旋转之后相同的算同一种,求不同染色方案数。
置换有四个:不动、转90°、转180°、转270°,在这些置换下不变的方案数共有16、2、4、2种,由burnside原理可以知道有(16+2+4+2)/4=6种方案。
接下来讲polya定理。首先,每个置换都可以被分解成若干不相交循环的乘积,什么意思呢?
例如上面的那个置换,我们注意到1->2->1,3->4->5->3,所以就可以写成(12)(345),(12)、(345)就叫循环。
我们发现在每个置换的作用下保持不变的染色方案,每个循环必须染一样的颜色,所以可以发现z就等于颜色数^循环节个数。
对应到上面那个题目,我们先把正方形四个格子按从左到右从上到下1234编号。
不动:(1)(2)(3)(4)。转90°:(1234)。转180°:(13)(24)。转270°:(1432)。
那么一共也是有(2^4+2^1+2^2+2^1)/4=6种方案。
例1 poj2409
用c种颜色染一个长度为s的环,旋转翻转之后相同的算同一种,求方案数,cs<=32(似乎只是为了保证答案不会太大)。
暴力硬上吧。由于某些原因,置换的循环节我都是用并查集找的你来打我啊
#include <iostream>
#include <stdio.h>
#include <math.h>
#include <string.h>
#include <time.h>
#include <stdlib.h>
#include <string>
#include <bitset>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <algorithm>
#include <sstream>
#include <stack>
#include <iomanip>
using namespace std;
#define pb push_back
#define mp make_pair
typedef pair<int,int> pii;
typedef long long ll;
typedef double ld;
typedef vector<int> vi;
#define fi first
#define se second
#define fe first
#define FO(x) {freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);}
#define Edg int M=0,fst[SZ],vb[SZ],nxt[SZ];void ad_de(int a,int b){++M;nxt[M]=fst[a];fst[a]=M;vb[M]=b;}void adde(int a,int b){ad_de(a,b);ad_de(b,a);}
#define Edgc int M=0,fst[SZ],vb[SZ],nxt[SZ],vc[SZ];void ad_de(int a,int b,int c){++M;nxt[M]=fst[a];fst[a]=M;vb[M]=b;vc[M]=c;}void adde(int a,int b,int c){ad_de(a,b,c);ad_de(b,a,c);}
#define es(x,e) (int e=fst[x];e;e=nxt[e])
#define VIZ {printf("digraph G{\n"); for(int i=1;i<=n;i++) for es(i,e) printf("%d->%d;\n",i,vb[e]); puts("}");}
#ifdef LOCAL
#define TIMER cerr<<clock()<<"ms\n"
#else
#define TIMER
#endif
#define SZ 666666
int ff[SZ];
void init(int n)
{
for(int i=0;i<=n;i++) ff[i]=-1;
}
int gf(int x)
{
return (ff[x]==-1)?x:ff[x]=gf(ff[x]);
}
void uni(int a,int b)
{
int ga=gf(a),gb=gf(b);
if(ga^gb) ff[ga]=gb;
}
ll qpl(ll a,ll b)
{
ll ans=1;
while(b)
{
if(b&1) ans=ans*a;
a=a*a; b>>=1;
}
return ans;
}
ll query(int c,int s)
{
ll ans=0;
for(int i=0;i<s;i++)
{
init(s);
for(int j=0;j<s;j++) uni(j,(j+i)%s);
int xh=0;
for(int j=0;j<s;j++) xh+=gf(j)==j;
ans+=qpl(c,xh);
}
for(int i=1;i<=s;i++)
{
init(s);
for(int j=0;j<s;j++) uni(j,((i-j)%s+s)%s);
int xh=0;
for(int j=0;j<s;j++) xh+=gf(j)==j;
ans+=qpl(c,xh);
}
return ans/(s*2);
}
int main()
{
int c,s;
while(1)
{
cin>>c>>s;
if(c==0&&s==0) break;
cout<<query(c,s)<<"\n";
}
}
然后这种做法显然过于暴力了,我们发现旋转k个位置之后的循环节就是gcd(s,k),翻转的话,如果奇数就是每个点到对边这样s个置换,每个置换循环节都是(s+1)/2,如果偶数就是对顶点和对边这样各s/2个,对顶点的循环节是s/2+1,对边的是s/2+1。
#include <iostream>
#include <stdio.h>
#include <math.h>
#include <string.h>
#include <time.h>
#include <stdlib.h>
#include <string>
#include <bitset>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <algorithm>
#include <sstream>
#include <stack>
#include <iomanip>
using namespace std;
#define pb push_back
#define mp make_pair
typedef pair<int,int> pii;
typedef long long ll;
typedef double ld;
typedef vector<int> vi;
#define fi first
#define se second
#define fe first
#define FO(x) {freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);}
#define Edg int M=0,fst[SZ],vb[SZ],nxt[SZ];void ad_de(int a,int b){++M;nxt[M]=fst[a];fst[a]=M;vb[M]=b;}void adde(int a,int b){ad_de(a,b);ad_de(b,a);}
#define Edgc int M=0,fst[SZ],vb[SZ],nxt[SZ],vc[SZ];void ad_de(int a,int b,int c){++M;nxt[M]=fst[a];fst[a]=M;vb[M]=b;vc[M]=c;}void adde(int a,int b,int c){ad_de(a,b,c);ad_de(b,a,c);}
#define es(x,e) (int e=fst[x];e;e=nxt[e])
#define VIZ {printf("digraph G{\n"); for(int i=1;i<=n;i++) for es(i,e) printf("%d->%d;\n",i,vb[e]); puts("}");}
#ifdef LOCAL
#define TIMER cerr<<clock()<<"ms\n"
#else
#define TIMER
#endif
#define SZ 666666
int ff[SZ];
void init(int n)
{
for(int i=0;i<=n;i++) ff[i]=-1;
}
int gf(int x)
{
return (ff[x]==-1)?x:ff[x]=gf(ff[x]);
}
void uni(int a,int b)
{
int ga=gf(a),gb=gf(b);
if(ga^gb) ff[ga]=gb;
}
ll qpl(ll a,ll b)
{
ll ans=1;
while(b)
{
if(b&1) ans=ans*a;
a=a*a; b>>=1;
}
return ans;
}
ll query(int c,int s)
{
ll ans=0;
for(int i=1;i<=s;i++)
ans=ans+qpl(c,__gcd(s,i));
int a,b;
if(s&1) a=s,b=0;
else a=b=s/2;
ans+=a*qpl(c,s/2+1);
ans+=b*qpl(c,s/2);
return ans/(s*2);
}
int main()
{
int c,s;
while(1)
{
cin>>c>>s;
if(c==0&&s==0) break;
cout<<query(c,s)<<"\n";
}
}
(gcd懒得写了,就用__gcd了)
例2 uestc75/tju2795 The Queen's New Necklaces
有一串项链用若干种颜色染色,第i种颜色有mi个,旋转相同的不计,问有多少种不同的项链。
我们考虑枚举旋转i个珠子之后不动的方案数,那么就会有gcd(n,i)个循环节,然后我们设p=n/gcd(n,i)为每个循环节的长度,那么如果有颜色数量不是p的倍数显然不可能有不动的,否则每个循环节颜色要一样,就随便统计一下就是了。
高精度懒得写,就不贴代码了。
例3 poj2154 color
用不超过n种颜色给一个长度为n的环染色,旋转相同的不计,求方案数mod P。n<=10^9,3500组数据。
似乎比例1还水,不过n这么大,显然需要优化。
ll ans=0;
for(int i=1;i<=n;i++)
ans=ans+qpl(n,__gcd(n,i));
return ans/n;
这是例1关于旋转那部分的代码,那么我们需要优化的实际上就是:$\sum_{i=1}^nn^{(i,n)}$
查询oeis可以发现:http://oeis.org/A228640
我们枚举(i,n)的gcd,那么(i/gcd,n/gcd)=1,所以对答案的贡献就是$\varphi(\frac{n}{gcd})*n^{gcd}$。
那么我们就得到了oeis上的这个公式:$\sum_{d|n}{\varphi(d)*n^{\frac{n}{d}}}$
然后我们发现这个式子似乎不是很好搞......
那么我们不妨枚举n的因数d,暴力算欧拉函数(拿所有质数筛一筛),这样效率应该还挺高的。
因为最后我们没办法除以n,所以算次方的时候要人工-1。
(由于懒没有写qaq有空写吧)
polya/burnside 学习的更多相关文章
- Burnside引理与Polya定理 学习笔记
原文链接www.cnblogs.com/zhouzhendong/p/Burnside-Polya.html 问题模型 有一个长度为 $n$ 的序列,序列中的每一个元素有 $m$ 种取值. 如果两个序 ...
- 【BZOJ】1004: [HNOI2008]Cards(置换群+polya+burnside)
http://www.lydsy.com/JudgeOnline/problem.php?id=1004 学习了下polya计数和burnside引理,最好的资料就是:<Pólya 计数法的应用 ...
- [POJ 2888]Magic Bracelet[Polya Burnside 置换 矩阵]
也许更好的阅读体验 \(\mathcal{Description}\) 大意:给一条长度为\(n\)的项链,有\(m\)种颜色,另有\(k\)条限制,每条限制为不允许\(x,y\)颜色连在一起.要求有 ...
- Polya 定理 学习笔记
群 群的定义 我们定义,对于一个集合 \(G\) 以及二元运算 \(\times\),如果满足以下四种性质,那我们就称 \((G,\times)\) 为一个群. 1. 封闭性 对于 \(a\in G, ...
- polya burnside 专题
polya题目:uva 11077 Find the Permutationsuva 10294 Arif in DhakaLA 3641 Leonardo's Notebookuva 11077 F ...
- 博客索引and题目列表
目录 笔记整理 计划 要学的东西 缺省源 要做的题 搜索 高斯消元 矩阵 排列组合 2019.7.9 2019.7.10 kmp ac自动机 2019.7.11 2019.7.15 笔记整理 1.同余 ...
- THUSC 2017 游记
Day0 早上在家里整理东西. 下午坐飞机去北京.(怎么又去北京,上周刚去的北京) 一开始飞机爬升的时候太无聊就睡着了.醒了以后就开始吃东西.吐槽一句:厦航的飞机就是好啊.上面的点心也比上次海航的好吃 ...
- 等价类计数(Polya定理/Burnside引理)学习笔记
参考:刘汝佳<算法竞赛入门经典训练指南> 感觉是非常远古的东西了,几乎从来没有看到过需要用这个的题,还是学一发以防翻车. 置换:排列的一一映射.置换乘法相当于函数复合.满足结合律,不满足交 ...
- 等价类计数:Burnside引理 & Polya定理
提示: 本文并非严谨的数学分析,有很多地方是自己瞎口胡的,仅供参考.有错误请不吝指出 :p 1. 群 1.1 群的概念 群 \((S,\circ)\) 是一个元素集合 \(S\) 和一种二元运算 $ ...
随机推荐
- windows 环境下nginx + tomcat群 + redis 实现session共享
nginx作为负载均衡根据定义将不同的用户请求分发到不同的服务器,同时也解决了因单点部署服务器故障导致的整个应用不能访问的问题 在加入nginx之后,如果多个服务器中的一个或多个(不是全部)发生故障, ...
- android 视频录制 混淆打包 之native层 异常的解决
原文地址:http://www.cnblogs.com/linguanh/ (滑至文章末,直接看解决方法) 问题起因: 前5天,因为项目里面有个类似 仿微信 视频录制的功能, 先是上网找了个 开 ...
- javaScript笔记
两边符号是tab键上面的那个键,不是单引号 即在sort内置的函数中先将各字符串转化为统一的大写或者小写,再进行比较即可. -------------------------------------- ...
- scikit-learn预处理实例之一:使用FunctionTransformer选择列
本例展示怎样在一个管道中使用FunctionTransformer.如果你知道你的数据集的第一主成分与分类任务无关,你可以使用FunctionTransformer选取除PCA转化的数据的第一列之外的 ...
- USB设备(移动硬盘、鼠标)掉电掉驱动的两种解决方案
症状: 当你发现"移动硬盘图标"经常无故消失,又自己出现时. 你可以把这个现象称之为"掉电" or "掉驱动". 遇到这种情况,相当不爽. ...
- ADO.NET存取数据库数据
步骤: //数据库连接串 string conStr = "......" //创建连接对象 SqlConnection connection = new SqlConnectio ...
- 如何实现一个php框架系列文章【4】url路由管理
直接通过url参数访问业务模块($app)中控制器($ctl)里的函数($act) 我们支持3种路由模式 普通模式 _a=$app, _u=$ctl.$act 最简单的方式,专注实现业务$ac ...
- 关于fefo函数
feof是C语言标准库函数函数,其原型在stdio.h中,其功能是检测流上的文件结束符. 函数原型: int feof(FILE *stream); 返回值:如果文件结束,则返回非0值,否则返回0 在 ...
- 设计模式03备忘录(java)
先贴代码有空来写内容. 备忘录1 //简单的备忘录,只可以记录上一次修改前的状态,实现撤回一次的操作. class Student{ private String name; private Stri ...
- mui 手势事件配置
在开发中监听双击屏幕事件时不起作用,需要在mui.init方法的gestureConfig参数中设置需要监听的手势事件 手势事件配置: 根据使用频率,mui默认会监听部分手势事件,如点击.滑动事件:为 ...