UVA10601 Cubes - 波利亚定理

Cubes

题意

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

题解

根据波利亚定理公式:

设X是元素集合,G是X的置换群,{u1,u2,...,uk}\{u_1,u_2,...,u_k\}{u1,u2,...,uk}kkk种颜色的集合,CCCXXX的任意着色集.这时,针对各颜色的数目的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 = 933=9个.
    (3)绕对边的中点连成的轴旋转180180180度:6∗1=66*1 = 661=6个.
    (4)绕体对角线旋转120,240120,240120,240度:4∗2=84*2=842=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+6z12z25)

  3. zi=∑u1iu2i...u6iz_i = \sum{u_1^iu_2^i...u_6^i}zi=u1iu2i...u6i代入PGP_GPG式.那么答案就是PGP_GPG函数的u1c1u2c2...ukcku_1^{c_1}u_2^{c_2}...u_k^{c_k}u1c1u2c2...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^zu1au2b...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;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/320993.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

COMCMS_CORE 起步篇,如何运行和部署

前言&#xff1a;关于最近开源后&#xff0c;不少朋友问&#xff0c;怎么我下载下来&#xff0c;运行不了。或者怎么没有左边菜单。货不对板&#xff1f;还是我吃了数据&#xff1f;感言&#xff1a;开源不容易&#xff0c;更不容易的是&#xff0c;明明毫无保留&#xff0c;还…

P4313-文理分科【最小割】

正题 题目链接:https://www.luogu.com.cn/problem/P4313 题目大意 有n∗mn*mn∗m个人&#xff0c;第(i,j)(i,j)(i,j)选择文科就可以获得arti,jart_{i,j}arti,j​的价值&#xff0c;选择理科就可以获得scii,jsci_{i,j}scii,j​的价值。如果一个选择文科的人周围都选择了文科&a…

初一模拟赛总结(3.30)

成绩&#xff1a; rankrankranknamenamenamescorescorescoreT1T1T1T2T2T2T3T3T3T4T4T4111lyflyflyf360360360100100100100100100606060100100100222wjjwjjwjj170170170100100100707070000000333fyfyfy167167167100100100505050000171717444hkyhkyhky11011011010010010000000010…

UVALive7670 Asa's Chess Problem,上下界费用流,另类解法

Asa’s Chess Problem 先阐述一下带上下界的边怎么建. 带上下界的建图方法 设我要建一条边(u→v)(u\rightarrow v)(u→v),流量上界为upupup,下界为downdowndown,费用为costcostcost.则我需要建两条边. 为保证一定会有downdowndown的流量流过去,我们可以建立一条u→vu \right…

ASP.NET Core WebApi使用Swagger生成api说明文档看这篇就够了

引言在使用asp.net core 进行api开发完成后&#xff0c;书写api说明文档对于程序员来说想必是件很痛苦的事情吧&#xff0c;但文档又必须写&#xff0c;而且文档的格式如果没有具体要求的话&#xff0c;最终完成的文档则完全取决于开发者的心情。或者详细点&#xff0c;或者简单…

P3201-[HNOI2009]梦幻布丁【启发式合并,链表】

正题 题目链接:https://www.luogu.com.cn/problem/P3201 题目大意 开始有nnn个布丁&#xff0c;第iii个是cic_ici​颜色的。 每次有操作 将所有颜色为xxx的布丁变为颜色yyy的。询问有多少个布丁颜色段。 解题思路 对于每次修改&#xff0c;我们可以考虑启发式合并&#xff…

【贪心】奶酪厂(jzoj 1285)

奶酪厂 题目大意&#xff1a; 有一个奶酪厂&#xff0c;每个星期&#xff08;共n个星期&#xff09;都有一定的单位生产成本和客户需求量&#xff0c;把奶酪保存一个星期每单位要s元&#xff0c;问一共花的钱最少是多少 Sample Input 4 5 88 200 89 400 97 300 91 500Sampl…

P3768 简单的数学题 [狄利克雷卷积,杜教筛,莫比乌斯反演]

简单的数学题 题目连接 https://www.luogu.org/problemnew/show/P3768 题目描述 输入一个正整数n,n≤1010n,n\le 10^{10}n,n≤1010和p,p≤1.1109p,p \le 1.1 \times 10^9p,p≤1.1109.且ppp为质数. 计算∑i1n∑j1nijgcd(i,j)\sum_{i1}^n\sum_{j1}^nijgcd(i,j)∑i1n​∑j1n​…

.NET Core微服务之基于Exceptionless实现分布式日志记录

一、Exceptionless极简介绍Exceptionless 是一个开源的实时的日志收集框架&#xff0c;它可以应用在基于 ASP.NET&#xff0c;ASP.NET Core&#xff0c;Web API&#xff0c;Web Forms&#xff0c;WPF&#xff0c;Console&#xff0c;ASP.NET MVC 等技术开发的应用程序中&#x…

P2801-教主的魔法【分块,二分】

正题 题目链接:https://www.luogu.com.cn/problem/P2801 题目大意 nnn个数字&#xff0c;要求支持 区间加上一个数字www询问一个区间内不小于www的数的个数 解题思路 考虑分块&#xff0c;对于块内我们维护一个排序后的数组&#xff0c;查询时直接在整个块中二分答案即可。修…

【多重背包】太空电梯(jzoj 1286)

太空电梯 题目大意&#xff1a; 有n&#xff08;1<n<400&#xff09;种石头&#xff0c;每种石头有它的数量c&#xff08;1<c<10&#xff09;&#xff0c;高度h&#xff08;1<h<100&#xff09;&#xff0c;可搭到的最高高度a&#xff08;1<a<40000&…

P2522 HAOI2011 Problem b [莫比乌斯反演,数论分块]

P2522 HAOI2011 题意 对于给出的n个询问&#xff0c;每次求有多少个数对(x,y)(x,y)(x,y)&#xff0c;满足a≤x≤ba≤x≤ba≤x≤b&#xff0c;c≤y≤dc≤y≤dc≤y≤d&#xff0c;且gcd(x,y)kgcd(x,y) kgcd(x,y)k&#xff0c;gcd(x,y)gcd(x,y)gcd(x,y)函数为xxx和yyy的最大公约…

.netcore 整合 log4net

1.背景前两天&#xff0c;曾经的一个同事咨询我&#xff0c;怎样将log4net以中间件的形式整合到core里边去。我不假思索的回答&#xff0c;这种问题应该有人做过吧&#xff0c;他说没有。于是&#xff0c;我去博客园搜了下&#xff0c;发现还真没有&#xff0c;全部都是传统.NE…

P5459-[BJOI2016]回转寿司【树状数组】

正题 题目链接:https://www.luogu.com.cn/problem/P5459 题目大意 nnn个数&#xff0c;求有多少个区间和在[L,R][L,R][L,R]范围内。 解题思路 显然我们做了前缀和之后&#xff0c;枚举右端点就只需要找到有多少个左端点满足在[x−R,x−L][x-R,x-L][x−R,x−L]这个范围内就好了…

【暴力】穹妹的求助

穹妹的求助 题目大意&#xff1a; 输入两个数&#xff0c;输出这两个数之间因数最多的数&#xff0c;和这个数的的因数个数 原题&#xff1a; 题目描述 由于穹妹很聪明&#xff0c;她的数学老师给她布置了一个作业&#xff0c;让她求出L到R之间不同因子数最多的那个数和这…

Asp.Net Core中利用Seq组件展示结构化日志功能

在一次.Net Core小项目的开发中&#xff0c;掌握的不够深入&#xff0c;对日志记录并没有好好利用&#xff0c;以至于一出现异常问题&#xff0c;都得跑动服务器上查看&#xff0c;那时一度怀疑自己肯定没学好&#xff0c;不然这一块日志不可能需要自己扒服务器日志来查看&…

P2633-Count on a tree【主席树,LCA】

正题 题目链接:https://www.luogu.com.cn/problem/P2633 题目大意 nnn个点的树&#xff0c;每个点有点权&#xff0c;求u∼vu\sim vu∼v的路径上第kkk小的权值&#xff0c;强制在线。 解题思路 考虑在树上维护主席树&#xff0c;我们不难发现如果sizxsiz_xsizx​表示根节点到…

P3327 约数的个数和 [约数函数性质,数论分块]

P3327 约数的个数和 题意 d(x)d(x)d(x)为约数的个数,对于每个询问,回答∑i1n∑j1md(ij)\sum_{i1}^n\sum_{j1}^md(ij)∑i1n​∑j1m​d(ij). 题解 这个题推得我头皮发麻,然后还没推出来,后来发现要做这题的先知道一个性质: d(ij)∑x∣i∑y∣j[gcd(x,y)1]d(ij)\sum_{x|i}\sum_{…

【结论】小X的矩阵

小X的矩阵 题目大意&#xff1a; 有一个nn的矩阵&#xff0c;要你执行g此操作&#xff0c;然后根据操作输出&#xff08;详情见原题&#xff09; 原题&#xff1a; 题目描述 小X最近迷上了矩阵&#xff0c;他定义了一个对于一种特殊矩阵的特征函数G。对于N∗NN*NN∗N的矩阵…

Apache SkyWalking的架构设计【译文】

Apache SkyWalking提供了一个功能强大并且很轻量级的后端。在此&#xff0c;将介绍为什么采用以下方式来设计它&#xff0c;以及它又是如何工作的。架构图对于APM而言&#xff0c;agent或SDKs仅是如何使用libs的技术细节。手动或自动的形式与架构无关&#xff0c;因此在本文中&…