Cubes

题意

给出121212根长度相等的木棒,颜色最多有666种,问能构成的本质不同的正方体数量.

题解

根据波利亚定理公式:

设X是元素集合,G是X的置换群,{u1,u2,...,uk}\{u_1,u_2,...,u_k\}{u1​,u2​,...,uk​}是kkk种颜色的集合,CCC是XXX的任意着色集.这时,针对各颜色的数目的C的非等价着色数的生成函数是:
PG(u1+...+uk,u12+...+uk2,...,u1n+...+ukn)P_G{(u_1+...+u_k,u_1^2+...+u_k^2,...,u_1^n+...+u_k^n)}PG​(u1​+...+uk​,u12​+...+uk2​,...,u1n​+...+ukn​)

  1. 我们先求出121212根木棒的颜色数,假设有kkk中,每种颜色假设有{c1,c2,...,ck}\{c_1,c_2,...,c_k\}{c1​,c2​,...,ck​}个,注意c1+...+ck=12c_1+...+c_k=12c1​+...+ck​=12.

  2. 求出置换群.
    众所周知的,正方体的置换群大小是242424.
    (1)恒等变换:111个
    (2)绕对面的中心旋转90,180,27090,180,27090,180,270度:3∗3=93*3 = 93∗3=9个.
    (3)绕对边的中点连成的轴旋转180180180度:6∗1=66*1 = 66∗1=6个.
    (4)绕体对角线旋转120,240120,240120,240度:4∗2=84*2=84∗2=8个.
    根据上述四种情况得到的关于线元素的置换平均值为:
    PG=124(z112+6z43+3z26+8z34+6z12z25)P_G=\frac{1}{24}(z_1^{12}+6z_4^3+3z_2^6+8z_3^4+6z_1^2z_2^5)PG​=241​(z112​+6z43​+3z26​+8z34​+6z12​z25​)

  3. 将zi=∑u1iu2i...u6iz_i = \sum{u_1^iu_2^i...u_6^i}zi​=∑u1i​u2i​...u6i​代入PGP_GPG​式.那么答案就是PGP_GPG​函数的u1c1u2c2...ukcku_1^{c_1}u_2^{c_2}...u_k^{c_k}u1c1​​u2c2​​...ukck​​前的系数

  4. PGP_GPG​函数的每一项都是形如(u1+...+uk)a(u12+...+uk2)b...(u1n+...+ukn)z(u_1+...+u_k)^{a}(u_1^2+...+u_k^2)^{b}...(u_1^n+...+u_k^n)^{z}(u1​+...+uk​)a(u12​+...+uk2​)b...(u1n​+...+ukn​)z的,我要对每个这样的式子,求u1au2b...ukzu_1^au_2^b...u_k^zu1a​u2b​...ukz​前的系数之和.每个项单独进行dfsdfsdfs就可以了.

代码

#include <iostream>
#include <algorithm>
#include <map>
#include <cstring>
#define pr(x) std::cout << #x << ':' << x << std::endl
#define rep(i,a,b) for(int i = a;i <= b;++i)int T;
int tc[7],x[7],ox[7];
int C[13][13];void init() {C[0][0] = 1;rep(i,1,12) {C[i][0] = 1;rep(j,1,i) {C[i][j] = C[i-1][j] + C[i-1][j-1];}}
}
int ans = 0;
void dfs(int dep,int mul) {if(dep == 5) {int f = 0;rep(i,1,6) f += x[i];if(!f) {ans += mul;}return ;}int tx[7];rep(i,1,6) tx[i] = x[i];for(int x1 = 0;x1 <= tc[dep] && dep * x1 <= tx[1];++x1) {int res1 = C[tc[dep]][x1];for(int x2 = 0;x2 <= tc[dep] - x1 && dep*x2 <= tx[2];++x2) {int res2 = res1 * C[tc[dep]-x1][x2];for(int x3 = 0;x3 <= tc[dep]-x1-x2 && dep*x3 <= tx[3];++x3) {int res3 = res2 * C[tc[dep]-x1-x2][x3];for(int x4 = 0;x4 <= tc[dep]-x1-x2-x3 && dep*x4 <= tx[4];++x4) {int res4 = res3 * C[tc[dep]-x1-x2-x3][x4];for(int x5 = 0;x5 <= tc[dep]-x1-x2-x3-x4 && dep*x5 <= tx[5];++x5) {int res5 = res4 * C[tc[dep]-x1-x2-x3-x4][x5];int x6 = tc[dep]-x1-x2-x3-x4-x5;//pr(x1);pr(x2);pr(x3);pr(x4);pr(x5);pr(x6);if(dep*x6 <= tx[6]) {x[1] = tx[1] - dep*x1;x[2] = tx[2] - dep*x2;x[3] = tx[3] - dep*x3;x[4] = tx[4] - dep*x4;x[5] = tx[5] - dep*x5;x[6] = tx[6] - dep*x6;dfs(dep+1,res5*mul);}}}}}}
}int main() {init();std::ios::sync_with_stdio(false);std::cin >> T;while(T--) {memset(ox,0,sizeof(ox));rep(i,1,12) {int tmp;std::cin >> tmp;ox[tmp] ++;}int res = 0;rep(i,1,6) x[i] = ox[i];ans = 0;tc[1] = 12,tc[2] = tc[3] = tc[4] = 0;dfs(1,1);res += ans;rep(i,1,6) x[i] = ox[i];ans = 0;tc[1] = tc[2] = tc[3] = 0,tc[4] = 3;dfs(1,1);res += 6*ans;rep(i,1,6) x[i] = ox[i];ans = 0;tc[1] = tc[3] = tc[4] = 0,tc[2] = 6;dfs(1,1);res += 3*ans;rep(i,1,6) x[i] = ox[i];ans = 0;tc[1] = tc[2] = tc[4] = 0,tc[3] = 4;dfs(1,1);res += 8*ans;rep(i,1,6) x[i] = ox[i];ans = 0;tc[1] = 2,tc[2] = 5,tc[3] = tc[4] = 0;dfs(1,1);res += 6*ans;std::cout << res/24 << std::endl;}return 0;
}

UVA10601 Cubes - 波利亚定理相关推荐

  1. 用数学方法计算井字棋合法局面数——波利亚定理的简单应用

    用数学方法计算井字棋合法局面数--波利亚定理的简单应用 0.前言 记不得那天在B站发现一个互动下井字棋智商普查(BV1JE411G71J),想当年被小学同学评价为无敌破战士(显然到现在都不知道这称号是 ...

  2. 【转】波利亚酒鬼回家定理的证明

    证明定理:喝醉的酒鬼总能找到回家的路,喝醉的小鸟则可能永远也回不了家. 假设有一条水平直线,从某个位置出发,每次有 50% 的概率向左走1米,有50%的概率向右走1米.按照这种方式无限地随机游走下去, ...

  3. 《怎样解题》-波利亚

    看一篇论文时发现这样一个有用的表<怎样解题表>,它出自波利亚的<怎样解题>一书 感觉很有感触啊啊... 1. 弄清问题 首先你必须弄清问题 z 未知数是什么?已知数是什么?条件 ...

  4. 跟波利亚学解题(rev#3)

    跟波利亚学解题(rev#3) By 刘未鹏(pongba) C++的罗浮宫(http://blog.csdn.net/pongba) TopLanguage(http://groups.google. ...

  5. 波利亚罐子模型及其应用

    1 引言 第1章介绍基本的离散分布和一些现代工具,例如随机过程和可交换性.为了符合本书的目标和范围,我们会尽快切换到骨灰盒语言.分布的实现通过urn参数得以证实,所有概念的示例均以urns来表示. 第 ...

  6. 转--总结波利亚在书中提到的思维方法

    总结波利亚在书中提到的思维方法,尤其是<How To Solve It>中的启发式思考方法,有这样一些: 时刻不忘未知量(即时刻别忘记你到底想要求什么,问题是什么.)莱布尼兹曾经将人的解题 ...

  7. 波利亚对教师日常工作的看法:〈教师十诫〉

    本文由何耿旭.陈彦宏翻译,洪志阳校订,原刊于<HPM通讯>.兹将该文略加修饰,转载于此(http://highscope.ch.ntu.edu.tw/wordpress/?p=18386) ...

  8. 罐子模型以及波利亚罐子模型(附加对论文的分析)

    1  罐子模型简介 1.1 基本模型 基本罐子模型中,罐子包含x个白色和y个黑色的球,它们混合在一起.从中中随机抽取一个球,观察其颜色:然后将其放回缸中(或不放回缸中),并重复选择过程. 在此模型中可 ...

  9. 波利亚(Polya)罐子模型

    (波利亚(Polya)罐子模型)罐中有a个白球,b个黑球,每次从罐中随机抽取一球,观察其颜色后,连同附加的c个同色球 (波利亚(Polya)罐子模型)罐中有a个白球,b个黑球,每次从罐中随机抽取一球, ...

最新文章

  1. python主辅线程_python主线程捕获子线程的方法
  2. 001_jQuery简介
  3. 自由自在意式手工冰淇淋式的生活方式
  4. 关于 To B,我欠粉丝们一个热搜
  5. 从人生(论文)的数量追求人生的质量
  6. java+解析未知json_在Java中解析JSON时如何忽略未知属性– Jackson @JsonIgnoreProperties注释示例...
  7. ubuntu下vim语法高亮问题
  8. Tensorflow入门之运算篇
  9. STL学习系列五:Queue容器
  10. 20155320 Exp3 免杀原理与实践
  11. python基础篇——列表与列表算法(下)
  12. 微信小程序自动定位城市
  13. 那是飞扬的青春—小评《DOOM启示录》
  14. C# 实现Excel导出图片
  15. 数据查询和业务流分开_基于大数据的舆情分析系统架构 - 架构篇
  16. 不只卖电脑的联想,还有什么“新看点”?
  17. Spring-IoC概念
  18. 网关技术详解及常见网关对比
  19. 武汉大学计算机学院 推免,2012年武汉大学计算机学院推免研究生计分细则
  20. PMP试题(435道)

热门文章

  1. linux删除第二次出现的字符,linux下 怎样删除文件名中包含特殊字符的文件
  2. c语言怎么让图形界面单独显示,「分享」C语言如何编写图形界面
  3. leetcode28. 实现 strStr(KMP详解)
  4. 7-1 内存分区分配--首次适应算法 (100 分)
  5. 用计算机探索 按一定规律排列,探索规律练习题1
  6. [JavaWeb-HTML]HTML标签_图片标签
  7. [Java基础]Set集合概述和特点
  8. [剑指offer]面试题10:二进制中1的个数
  9. 最长公共子序列-dp
  10. [蓝桥杯2016初赛]平方怪圈-数论,模拟