\(\mathcal{Description}\)

  Link.

  求含 \(n\) 个点的无标号简单无向图的个数,答案模 \(997\)。

\(\mathcal{Solution}\)

  首先,把题目转化成为有标号 \(K_n\) 的 \(\binom{n}{2}\) 条边染黑(不选)白(选)两种颜色,求本质不同(去除标号)的方案数。想到使用 Pólya 定理求解。设在某个点转置中,循环大小为 \(a_1,a_2,\cdots,a_k\),分别考虑循环内部和循环间的边等价类:

  对于大小为 \(a\) 的循环内部,仍旧是一个 \(K_a\),尝试求其中边等价类的个数,例如 \(a=6\),如图(来自 @pythoner713 的博客,下图同):

可见共红色、绿色、蓝色三个等价类。进一步,若把 \(K_a\) 画作类似的正多边形,发现边属于同一等价类当且仅当它们长度相等,继而得出大小为 \(a\) 的循环的等价类个数为 \(\lfloor\frac{a}{2}\rfloor\)。

  对于两个大小分别为 \(a,b\) 的循环之间,考虑转置它们直到与初始状态重合,如图:

循环大小为 \(\operatorname{lcm}(a,b)\),那么共 \(\gcd(a,b)\) 个等价类。

  综上,符合 \(\{a_k\}\) 所描述的转置的等价类个数 \(c(a)\) 为

\[\sum_{i=1}^k\lfloor\frac{a_i}{2}\rfloor+\sum_{1\le i<j\le k}\gcd(a_i,a_j)
\]

  考虑直接枚举 \(\{a_k\}\),则符合条件的转置方案数为

\[\frac{n!}{\prod_{i=1}^ka_i\prod_{i=1}^nb_i}
\]

其中 \(b_i=\sum_{j=1}^k[a_j=i]\),即大小为 \(i\) 的循环个数。组合意义为:在所有转置中,固定每个循环的第一个元素,再去除等大循环的顺序。

  最终方案数为

\[\frac{1}{n!}\sum_{\{a_k\}}2^{c(a)}\frac{n!}{\prod_{i=1}^ka_i\prod_{i=1}^nb_i}\\=\sum_{\{a_k\}}2^{c(a)}\left(\prod_{i=1}^ka_i\prod_{i=1}^nb_i\right)^{-1}
\]

拆分 \(n\) 求解即可,复杂度见 A296010

\(\mathcal{Code}\)

/* Clearink */

#include <cstdio>

#define rep( i, l, r ) for ( int i = l, rep##i = r; i <= rep##i; ++i )
#define per( i, r, l ) for ( int i = r, per##i = l; i >= per##i; --i ) const int MOD = 997, MAXN = 60;
int n, fac[MAXN + 5]; inline int gcd( const int a, const int b ) { return b ? gcd( b, a % b ) : a; }
inline int mul( const long long a, const int b ) { return int( a * b % MOD ); }
inline int add( int a, const int b ) { return ( a += b ) < MOD ? a : a - MOD; }
inline void addeq( int& a, const int b ) { ( a += b ) >= MOD && ( a -= MOD ); }
inline int mpow( int a, int b ) {
int ret = 1;
for ( ; b; a = mul( a, a ), b >>= 1 ) ret = mul( ret, b & 1 ? a : 1 );
return ret;
} int a[MAXN + 5], ans, t; inline int calc() {
t += a[0] * a[0];
int s = 1;
for ( int l = 1, r; l <= a[0]; l = r ) {
for ( r = l; r <= a[0] && a[l] == a[r]; ++r );
s = mul( s, fac[r - l] );
} int idx = 0;
rep ( i, 1, a[0] ) {
idx += a[i] >> 1, s = mul( s, a[i] );
rep ( j, i + 1, a[0] ) idx += gcd( a[i], a[j] );
} return mul( mpow( s, MOD - 2 ), mpow( 2, idx ) );
} inline void split( const int x, const int rest ) {
if ( !rest ) return addeq( ans, calc() );
if ( rest < x ) return ;
split( x + 1, rest );
a[++a[0]] = x, split( x, rest - x ), --a[0];
} int main() {
scanf( "%d", &n ); fac[0] = 1;
rep ( i, 1, n ) fac[i] = mul( i, fac[i - 1] ); split( 1, n ); printf( "%d\n", t );
return 0;
}

Solution -「HNOI 2009」「洛谷 P4727」图的同构计数的更多相关文章

  1. 「区间DP」「洛谷P1043」数字游戏

    「洛谷P1043」数字游戏 日后再写 代码 /*#!/bin/sh dir=$GEDIT_CURRENT_DOCUMENT_DIR name=$GEDIT_CURRENT_DOCUMENT_NAME ...

  2. 洛谷 P1596 [USACO10OCT]湖计数Lake Counting

    题目链接 https://www.luogu.org/problemnew/show/P1596 题目描述 Due to recent rains, water has pooled in vario ...

  3. 洛谷P1144 最短路计数(SPFA)

    To 洛谷.1144 最短路计数 题目描述 给出一个N个顶点M条边的无向无权图,顶点编号为1-N.问从顶点1开始,到其他每个点的最短路有几条. 输入输出格式 输入格式: 输入第一行包含2个正整数N,M ...

  4. 洛谷 1144 最短路计数 bfs

    洛谷1144 最短路计数 传送门 其实这道题目的正解应该是spfa里面加一些处理,,然而,,然而,,既然它是无权图,,那么就直接bfs了,用一个cnt记录一下每一个点的方案数,分几种情况讨论一下转移, ...

  5. 洛谷 P4017 最大食物链计数

    洛谷 P4017 最大食物链计数 洛谷传送门 题目背景 你知道食物链吗?Delia生物考试的时候,数食物链条数的题目全都错了,因为她总是重复数了几条或漏掉了几条.于是她来就来求助你,然而你也不会啊!写 ...

  6. 动态规划 洛谷P4017 最大食物链计数——图上动态规划 拓扑排序

    洛谷P4017 最大食物链计数 这是洛谷一题普及/提高-的题目,也是我第一次做的一题 图上动态规划/拓扑排序 ,我认为这题是很好的学习拓扑排序的题目. 在这题中,我学到了几个名词,入度,出度,及没有环 ...

  7. Solution -「JSOI 2019」「洛谷 P5334」节日庆典

    \(\mathscr{Description}\)   Link.   给定字符串 \(S\),求 \(S\) 的每个前缀的最小表示法起始下标(若有多个,取最小的).   \(|S|\le3\time ...

  8. Solution -「洛谷 P4372」Out of Sorts P

    \(\mathcal{Description}\)   OurOJ & 洛谷 P4372(几乎一致)   设计一个排序算法,设现在对 \(\{a_n\}\) 中 \([l,r]\) 内的元素排 ...

  9. Solution -「POI 2010」「洛谷 P3511」MOS-Bridges

    \(\mathcal{Description}\)   Link.(洛谷上这翻译真的一言难尽呐.   给定一个 \(n\) 个点 \(m\) 条边的无向图,一条边 \((u,v,a,b)\) 表示从 ...

随机推荐

  1. Win10+Java7环境配置

    原文链接: https://www.toutiao.com/i6487838676326810125/ 安装包: jdk-7u79-windows-x64-20151024 打开安装目录 双击运行程序 ...

  2. 获取js代码运行的时间

    console.time('yue') //代码部分 console.timeEnd('yue')

  3. kali linux2020 虚拟机改root密码

    kali在2020版的更新中,好多小伙伴登不进root账号,这里来教大家怎样改root账户的密码 1.当我们打开虚拟机看到这个界面的时候,按e进入编辑模式 2.在编辑模式中,"quite s ...

  4. AVD模拟器怎么配置上网

    转自:http://blog.csdn.net/you_jinjin/article/details/7228303 方法一 首先,Windows下,配置Adroid环境变量(Win7为例) 1.桌面 ...

  5. 《剑指offer》面试题44. 数字序列中某一位的数字

    问题描述 数字以0123456789101112131415-的格式序列化到一个字符序列中.在这个序列中,第5位(从下标0开始计数)是5,第13位是1,第19位是4,等等. 请写一个函数,求任意第n位 ...

  6. porcupine语音唤醒python实现

    note it is not for arm pyaudio <= 3.6 version porcupine 3.5 3.6 not 3.7 code import struct import ...

  7. 推荐召回--基于物品的协同过滤:ItemCF

    目录 1. 前言 2. 原理&计算&改进 3. 总结 1. 前言 说完基于用户的协同过滤后,趁热打铁,我们来说说基于物品的协同过滤:"看了又看","买了又 ...

  8. python29day

    内容回顾 网络编程 概念 B/S C/S架构 B/S browser server C/S client 装客户端使用的 server远程服务器的 osi七层协议 今日内容 tcp协议的编程 如何在连 ...

  9. 【Vue源码学习】依赖收集

    前面我们学习了vue的响应式原理,我们知道了vue2底层是通过Object.defineProperty来实现数据响应式的,但是单有这个还不够,我们在data中定义的数据可能没有用于模版渲染,修改这些 ...

  10. MySQL数据库索引介绍

    一.什么是索引 索引是mysql数据库中的一种数据结构,就是一种数据的组织方式,这种数据结构又称为key 表中的一行行数据按照索引规定的结构组织成了一种树型结构,该树叫B+树 二.为何要用索引 优化查 ...