题意

给你一个长为\(n\)的项链和手镯,每个珠子有\(m\)种颜色.

两个手镯定义为相同,即它们通过翻转和旋转得到一样的手镯.

两个项链定义为相同,即它们只能通过旋转得到一样的项链.

求出有多少种本质不同的项链和手镯.

\((1 \le n \le 50, 1 \le m \le 10)\)

题解

(参考了一下这篇大佬博客)

大白书上的原题,一个裸的Polya定理(逃

Polya定理 : $$L=\frac{1}{|G|}\sum \limits _{i=1} ^{|G|} m^{c(g_i)}$$

其中\(G=\{g_1, ..., g_s\}\) \(c(g_i)\)为置换\(g_i\)的循环节个数, \(L\)为本质不同的方案数, \(m\) 是可以选择的颜色种数.

  • 首先考虑旋转 :

    假设当前旋转\(i\)颗珠子,那么就有\(\gcd (n,i)\)个循环,每个循环长度则为\(\displaystyle \frac{n}{\gcd(n,i)}\).

    这个证明同样参考了之前那篇博客....(我用 \(\LaTeX\) 再打一下..)

    将珠子从\(0\)到\(n-1\)标号,那么对于旋转\(i\)位的置换,在以\(0\)号为起点,长度为\(t\)的一个循环节,

    元素标号就为\(0,i \bmod n, 2i \bmod n, ... , (t-1)i \bmod n\).

    所以就有\(t \cdot i \bmod n = 0\),即有\(t \cdot i = n \cdot k\). 使等式左右成立的最小正整数就为\(\mathrm{lcm} (n, i)\).

    那么\(t \cdot i = \mathrm{lcm}(n,i)\)所以\(\displaystyle t=\frac{\mathrm{lcm}(n,i)}{i}=\frac{n}{\gcd(n,i)}\). 那么循环节个数就是\(\displaystyle \frac{n}{t}=\gcd(n,i)\).

    所以这些置换的贡献就是\(a=\sum \limits_{i=0}^{n-1} m^{\gcd(i,n)}\).

  • 再考虑一下翻转 :

    这个要分序列奇偶性分别考虑,应该是比较好考虑的.

    1. 如果长度\(n\)为奇数,那么我们就只能沿着一个珠子翻转,那么就有\(\frac{n-1}{2}\)个长度为\(2\)的等价类 和 \(1\)个长度为\(1\)的等价类.
    2. 如果长度\(n\)为偶数,那么我们就有两种方式.沿着两个珠子翻转,那么就有\(\frac{n-2}{2}\)个长度为\(2\)的等价类和\(2\)个长度为\(1\)的等价类. 否则我们沿着两个珠子中间来翻转,那么就有\(\frac{n}{2}\)个长度为\(2\)的等价类.

    所以这些置换的贡献就是 \(\displaystyle b=\begin{cases}n \cdot m^{\frac{n+1}{2}} & (m \ is \ odd)\\ \frac{n}{2}\cdot (m ^ {\frac{n}{2}+1}+m^{\frac{n}{2}}) &(m\ is\ even)\end{cases}\)

这样的话似乎会漏算情况(hany01大大问了下我)然后还是前面那篇大佬博客上有解释.

就是旋转再翻转的情况,肯定是其中另一种翻转.这也是因为群有封闭性,不管怎样\(*\)(二元运算)都是群内的元素.

最后要除以一个\(|G|\),项链只有\(n\)个置换群,手镯有\(2n\)个啦.

代码

#include <bits/stdc++.h>
#define For(i, l, r) for(register int i = (l), _end_ = (int)(r); i <= _end_; ++i)
#define Fordown(i, r, l) for(register int i = (r), _end_ = (int)(l); i >= _end_; --i)
#define Set(a, v) memset(a, v, sizeof(a))
using namespace std; typedef long long ll;
bool chkmin(int &a, int b) {return b < a ? a = b, 1 : 0;}
bool chkmax(int &a, int b) {return b > a ? a = b, 1 : 0;} void File() {
#ifdef zjp_shadow
freopen ("10294.in", "r", stdin);
freopen ("10294.out", "w", stdout);
#endif
} const int N = 55;
ll n, t, Pow[N], a, b; int main () {
File();
while (~scanf("%lld%lld", &n, &t) && n) {
Pow[0] = 1;
a = b = 0;
for (int i = 1; i <= n; ++i) Pow[i] = Pow[i - 1] * t;
for (int i = 0; i < n; ++ i) a += Pow[__gcd(i, (int)n)];
if (n & 1) b = n * Pow[n / 2 + 1];
else b = n / 2 * (Pow[n / 2 + 1] + Pow[n / 2]);
printf ("%lld %lld\n", a / n, (a + b) / 2 / n);
}
return 0;
}

UVA10294 Arif in Dhaka (群论,Polya定理)的更多相关文章

  1. [Uva10294]Arif in Dhaka

    [Uva10294]Arif in Dhaka 标签: 置换 Burnside引理 题目链接 题意 有很多个珠子穿成环形首饰,手镯可以翻转和旋转,项链只能旋转.(翻转过的手镯相同,而项链不同) 有n个 ...

  2. Uva10294 Arif in Dhaka (置换问题)

    扯回正题,此题需要知道的是置换群的概念,这点在刘汝佳的书中写的比较详细,此处不多做赘述.此处多说一句的是第二种手镯的情况.在下图中“左图顺时针转1个位置”和“右图顺时针旋转5个位置”是相同的,所以在最 ...

  3. UVA10294 Arif in Dhaka (First Love Part 2) —— 置换、poyla定理

    题目链接:https://vjudge.net/problem/UVA-10294 题解: 白书P146~147. 为什么旋转i个间距,就有gcd(i,n)个循环,且每个循环有n/gcd(i,n)个元 ...

  4. UVA10294 Arif in Dhaka (First Love Part 2)

    题意 PDF 分析 用n颗宝石串成项链和手镯, 每颗宝石的颜色可以t种颜色中的一种,当A类项链经过旋转得B类项链时,A和B属于一类项链, 而手镯不仅可以旋转还可以翻转,当A类手镯经过翻转得得到B类手镯 ...

  5. 【uva 10294】 Arif in Dhaka (First Love Part 2) (置换,burnside引理|polya定理)

    题目来源:UVa 10294 Arif in Dhaka (First Love Part 2) 题意:n颗珠子t种颜色 求有多少种项链和手镯 项链不可以翻转 手镯可以翻转 [分析] 要开始学置换了. ...

  6. Arif in Dhaka (First Love Part 2) UVA - 10294(Polya定理)

    这题和POJ-1286一样 题意: 给出t种颜色的n颗珠子 (每种颜色的珠子个数无限制,但总数必须是n), 求能制作出项链和手镯的个数 注意手镯可以翻转和旋转  而 项练只能旋转 解析: 注意Poly ...

  7. UVa 10294 Arif in Dhaka (First Love Part 2) (Polya定理)

    题意:给定 n 和 m 表示要制作一个项链和手镯,项链和手镯的区别就是手镯旋转和翻转都是相同的,而项链旋转都是相同的,而翻转是不同的,问你使用 n 个珠子和 m 种颜色可以制作多少种项链和手镯. 析: ...

  8. 【群论】polya定理

    对Polya定理的个人认识     我们先来看一道经典题目:     He's Circles(SGU 294)         有一个长度为N的环,上面写着“X”和“E”,问本质不同的环有多少个(不 ...

  9. [bzoj1488][HNOI2009]图的同构——Polya定理

    题目大意 求两两互不同构的含n个点的简单图有多少种. 简单图是关联一对顶点的无向边不多于一条的不含自环的图. a图与b图被认为是同构的是指a图的顶点经过一定的重新标号以后,a图的顶点集和边集能完全与b ...

随机推荐

  1. 乞丐版servlet容器第1篇

    本系列参照pkpk1234大神的BeggarServletContainer,具体请访问:https://github.com/pkpk1234/BeggarServletContainer. 一步一 ...

  2. 调试 smallcorgi/Faster-RCNN_TF 的demo过程遇到的问题

    最近在调试faster R-CNN时,遇到了各种各样的问题.使用的算法库为https://github.com/smallcorgi/Faster-RCNN_TF 注:本文使用的是通过virtuale ...

  3. win7(windows 7)系统下安装SQL2005(SQL Server 2005)图文教程

    操作系统:Microsoft Windows 7 旗舰版(32位) 数据库版本:SQL Server 2005 简体中文开发板 数据库下载链接: https://pan.baidu.com/s/1cq ...

  4. MySQL复制同一个服务器的表结构和表数据

    例如,现在服务器上有数据库 dbx 和 dby,dbx中有很多表,要把dbx中的表全部复制到dby,如下操作: 首先: use dby; [复制表结构] CREATE TABLE user LIKE ...

  5. SqlBulkCopy 参数配置示例

    SqlBulkCopy  做为SQL Server 官方 批量入库类,性能不会太差.针对其参数做了一些测试.   A. 先准备测试场景 ,关于SqlBulkCopyOptions.KeepIdenti ...

  6. haproxy下X-Frame-Options修复方法

    X-Frame-Options是网页常见漏洞之一,在其他的frame下加载你网页的内容,达到劫持的目的. 修复方法很简单,在http 请求投中添加X-Frame-Options参数.网上可以查一下Ap ...

  7. 2015年百度实习生前端笔试题上海卷a

    1.写出javascript运行结果:alert(‘5’+5); 结果:’55’ 2.写出javascript运行结果:for(var i=0; i<10; i++){} alert(i); 结 ...

  8. JS标签的各种事件的举例

    1.鼠标单击事件( onclick ) <!DOCTYPE HTML> <html> <head> <meta http-equiv="Conten ...

  9. 业余草分享100套精选1000G架构师资料课程(超1T的IT学习资料免费送)

    业余草分享100套精选1000G架构师资料课程(超1T的IT学习资料免费送). 超过1024G的IT学习资料免费领取,你值得拥有! 领取资源方式,关注“业余草”公众号,回复对应的关键字 01.回复”我 ...

  10. 计蒜客 无脑博士 bfs

    题目链接无脑博士的试管们 思路:直接模拟倒水过程即可,但是需要记忆判断当前的情况是否已经处理过.dfs和bfs都ok AC代码 #include <cstdio> #include < ...