bzoj1488 [HNOI2009]图的同构 Burnside 引理

最新推荐文章于 2024-08-02 16:33:41 发布
最新推荐文章于 2024-08-02 16:33:41 发布
阅读量96 收藏
点赞数
文章标签： 数据结构与算法
原文链接：http://www.cnblogs.com/hankeke/p/bzoj1488.html
版权

题目传送门

bzoj1488 - [HNOI2009]图的同构

bzoj1815 - [Shoi2006]color 有色图（双倍经验）

题解

暴力

由于在做题之前已经被告知是 Burnside 引理，貌似思考的时候少了一些乐趣啊。

考虑一个置换 \(p\)，想要求出这个置换下的不动点的个数。对于一个不动点，若存在一条边 \((a, b)\)，一定存在一条边 \((p_a, p_b)\)。

那么考虑一个长度为 \(l\) 的循环，若 \((i, j)\) 是一条 \(i, j\) 均在循环中的点内部的边，那么所有距离等于 \(i, j\) 的距离的点对之间都会有一条边。这样，我们一共可以将这些边分成 \(\frac x2\) 中内部绑定的组。因此，一个长度为 \(l\) 的循环可以产生 \(2^{\frac x2}\) 种方案的贡献。

下面考虑一条两个端点分别在长度为 \(x\) 的循环和 \(y\) 的循环上的边。这样的话，一条边的两个端点在两个循环上走，一直到下一次两个端点都与原来的两个端点重合。显然对于一组点，会走过 \(\operatorname{lcm}(x, y)\) 个点。于是，一共有 \(\frac {xy}{\operatorname{lcm}(x, y)} = \gcd(x, y)\) 组。贡献为 \(2^{\gcd(x, y)}\) 种方案。

优化

但是暴力算法显然会 TLE。

我们一个置换的贡献仅取决于它的每一个循环的长度，因此我们只需要知道它的每一个循环的长度就可以求出贡献了。

于是我们可以使用 Burnside 引理相关题目的常见套路。

我们用 \(B(n)\) 来表示整数划分的方案数。整数划分是指将 \(n\) 划分为若干个无序的正整数的和。当 \(n=100\) 时，\(B(n)\) 大概是 1e6。当 \(n = 40\) 时大概是 \(4e4\) 左右，\(n=45\) 时候大约 \(9e4\)。

因此，由于本题 \(n=45\)，如果暴力枚举划分的话，不会超过 \(9e4\) 次。

那么我们如果枚举出 \(n=\sum \limits_{i=1}^k L_i\)，那么方案相当于是从 \(n\) 个位置中选出 \(L_1\) 个位置作为第一个循环的成分，再从剩下的 \(n-L_1\) 个位置中选出 \(L_2\) 个……

因此对于这样的划分，其方案数是：
\[ \begin{align*} &\binom n{L_1} \binom{n-L_1}{L2} \binom{n-L_1-L_2}{L_2}\cdots\\ =&\frac{n!}{\prod \limits_{i=1}^k Li!} \end{align*} \]
（实际上也可以从可重复元素的全排列的角度考虑

然后求出来划分的方案数以后，我们还需要把每一个数放进这些被安排好的位置，使得每一段真的是一个循环。一个循环的第一个数可以考虑连接到这一段剩余的 \(L-1\) 个位置的任何一个位置。然后对于刚刚选择的那个位置，它又可以连接到剩余的 \(L-2\) 个位置。因此对于一个长度为 \(L\) 的循环，它的合法安排数字的方案数为 \((L-1)!\)。

但是还有一个小问题，如果有两个长度相同的置换，他们是没有什么先后区别的，但是在用之前的组合计算方法的时候，会产生先后的区别。这是不应该产生的。因此令 \(C_i\) 表示等于 \(i\) 的 \(L\) 值的个数。我们还需要乘上 \(C_i!\)。

因此对于一个整数划分的总方案数为：
\[ \frac {n!}{\prod\limits_{i=1}^k L_i \prod\limits_{i=1}^n C_i!} \]
然后就是和暴力一样的计算不动点的个数了。

最终答案为：
\[ \sum_{L} \frac {\prod \limits_{i=1}^k 2^{\frac {L_i}2}\prod \limits_{i=1}^k \prod \limits_{j=i+1}^k 2^{\gcd(L_i, L_j)} }{\prod\limits_{i=1}^k L_i \prod\limits_{i=1}^n C_i!} \]

下面是代码。由于 \(\gcd\) 和 \(2\) 的正整数幂均可以预处理，因此总的时间复杂度为 \(O(B_n\cdot n^2)\).

#include<bits/stdc++.h>

#define fec(i, x, y) (int i = head[x], y = g[i].to; i; i = g[i].ne, y = g[i].to)
#define dbg(...) fprintf(stderr, __VA_ARGS__)
#define File(x) freopen(#x".in", "r", stdin), freopen(#x".out", "w", stdout)
#define fi first
#define se second
#define pb push_back

template<typename A, typename B> inline char smax(A &a, const B &b) {return a < b ? a = b , 1 : 0;}
template<typename A, typename B> inline char smin(A &a, const B &b) {return b < a ?  a = b , 1 : 0;}

typedef long long ll; typedef unsigned long long ull; typedef std::pair<int, int> pii;

template<typename I>
inline void read(I &x) {
    int f = 0, c;
    while (!isdigit(c = getchar())) c == '-' ? f = 1 : 0;
    x = c & 15;
    while (isdigit(c = getchar())) x = (x << 1) + (x << 3) + (c & 15);
    f ? x = -x : 0;
}

const int N = 200 + 7;
const int P = 997;

int n, ans, hkk;
int inv[N], fac[N], ifac[N << 1], pw[N], gcd[N][N];
int l[N];

inline void sadd(int &x, int y) { x += y; x >= P ? x -= P : x; }
inline int smod(int x) { return x >= P ? x - P : x; }

inline void ycl() {
    fac[0] = 1; for (int i = 1; i <= n; ++i) fac[i] = (ll)fac[i - 1] * i % P;
    inv[1] = 1; for (int i = 2; i <= n; ++i) inv[i] = (ll)(P - P / i) * inv[P % i] % P;
    ifac[0] = 1; for (int i = 1; i <= n; ++i) ifac[i] = (ll)ifac[i - 1] * inv[i] % P;
    pw[0] = 1; for (int i = 1; i <= n; ++i) pw[i] = smod(pw[i - 1] << 1);
    for (int i = 0; i <= n; ++i) {
        for (int j = 0; j <= i; ++j) {
            if (!i || !j) gcd[i][j] = i ^ j;
            else gcd[i][j] = gcd[j][i % j];
            gcd[j][i] = gcd[i][j];
//          dbg("%d %d %d\n", i, j, gcd[i][j]);
        }
    }
}

inline void calc(int div) {
    int cnt = 1;
    for (int i = 1; i <= l[0]; ++i) cnt = (ll)cnt * pw[l[i] / 2] % P;//, dbg("l[%d] = %d, %d\n", i, l[i], pw[l[i] / 2]);
    for (int i = 1; i <= l[0]; ++i)
        for (int j = 1; j < i; ++j) cnt = (ll)cnt * pw[gcd[l[i]][l[j]]] % P;//, dbg("*****************");
    sadd(ans, (ll)cnt * div % P);
}

inline void dfs(int n, int lim, int div) {
//  dbg("n = %d, lim = %d, div = %d\n", n, lim, div);
    if (!n) return calc(div);
    for (int i = std::min(lim, n); i; --i) {
        int df = div;
        for (int j = 1; i * j <= n; ++j) {
            df = (ll)df * inv[i] % P, l[++l[0]] = i;
            dfs(n - i * j, i - 1, (ll)df * ifac[j] % P);
        }
        while (l[l[0]] == i) --l[0];
    }
}

inline void work() {
    ycl();
    ans = 0;
    dfs(n, n, 1);
    printf("%d\n", ans);
}

inline void init() {
    read(n);
}

int main() {
#ifdef hzhkk
    freopen("hkk.in", "r", stdin);
#endif
    init();
    work();
    fclose(stdin), fclose(stdout);
    return 0;
}

附录 - bzoj1815 - [Shoi2006]color 有色图 代码

发现只需要把上面的有/无的状态改为 \(k\) 中颜色就可以了。

#include<bits/stdc++.h>

#define fec(i, x, y) (int i = head[x], y = g[i].to; i; i = g[i].ne, y = g[i].to)
#define dbg(...) fprintf(stderr, __VA_ARGS__)
#define File(x) freopen(#x".in", "r", stdin), freopen(#x".out", "w", stdout)
#define fi first
#define se second
#define pb push_back

template<typename A, typename B> inline char smax(A &a, const B &b) {return a < b ? a = b , 1 : 0;}
template<typename A, typename B> inline char smin(A &a, const B &b) {return b < a ? a = b , 1 : 0;}

typedef long long ll; typedef unsigned long long ull; typedef std::pair<int, int> pii;

template<typename I>
inline void read(I &x) {
    int f = 0, c;
    while (!isdigit(c = getchar())) c == '-' ? f = 1 : 0;
    x = c & 15;
    while (isdigit(c = getchar())) x = (x << 1) + (x << 3) + (c & 15);
    f ? x = -x : 0;
}

const int N = 53 + 7;

int n, k, P, ans;
int inv[N], fac[N], ifac[N << 1], pw[N], gcd[N][N];
int l[N];

inline void sadd(int &x, int y) { x += y; x >= P ? x -= P : x; }
inline int smod(int x) { return x >= P ? x - P : x; }

inline void ycl() {
    fac[0] = 1; for (int i = 1; i <= n; ++i) fac[i] = (ll)fac[i - 1] * i % P;
    inv[1] = 1; for (int i = 2; i <= n; ++i) inv[i] = (ll)(P - P / i) * inv[P % i] % P;
    ifac[0] = 1; for (int i = 1; i <= n; ++i) ifac[i] = (ll)ifac[i - 1] * inv[i] % P;
    pw[0] = 1; for (int i = 1; i <= n; ++i) pw[i] = (ll)pw[i - 1] * k % P;
    for (int i = 0; i <= n; ++i) {
        for (int j = 0; j <= i; ++j) {
            if (!i || !j) gcd[i][j] = i ^ j;
            else gcd[i][j] = gcd[j][i % j];
            gcd[j][i] = gcd[i][j];
//          dbg("%d %d %d\n", i, j, gcd[i][j]);
        }
    }
}

inline void calc(int div) {
    int cnt = 1;
    for (int i = 1; i <= l[0]; ++i) cnt = (ll)cnt * pw[l[i] / 2] % P;//, dbg("l[%d] = %d, %d\n", i, l[i], pw[l[i] / 2]);
    for (int i = 1; i <= l[0]; ++i)
        for (int j = 1; j < i; ++j) cnt = (ll)cnt * pw[gcd[l[i]][l[j]]] % P;//, dbg("*****************");
    sadd(ans, (ll)cnt * div % P);
}

inline void dfs(int n, int lim, int div) {
//  dbg("n = %d, lim = %d, div = %d\n", n, lim, div);
    if (!n) return calc(div);
    for (int i = std::min(lim, n); i; --i) {
        int df = div;
        for (int j = 1; i * j <= n; ++j) {
            df = (ll)df * inv[i] % P, l[++l[0]] = i;
            dfs(n - i * j, i - 1, (ll)df * ifac[j] % P);
        }
        while (l[l[0]] == i) --l[0];
    }
}

inline void work() {
    ycl();
    ans = 0;
    dfs(n, n, 1);
    printf("%d\n", ans);
}

inline void init() {
    read(n), read(k), read(P);
}

int main() {
#ifdef hzhkk
    freopen("hkk.in", "r", stdin);
#endif
    init();
    work();
    fclose(stdin), fclose(stdout);
    return 0;
}

转载于:https://www.cnblogs.com/hankeke/p/bzoj1488.html

weixin_30455023
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
【BZOJ1488】【HNOI2009】图的同构
cz_xuyixuan的博客
01-11 710
【题目链接】 点击打开链接 【思路要点】 考虑Burnside引理或Polya定理，求解总共\(N!\)种点置换中边的染色方案等价类的个数总和，再将答案除以N！。考虑一个点集的置换\(P_{i}\)，若将一个置换看做一张由\(N\)个点、\(N\)条边的无向图，那么显然，该图由若干个环组成。记其中每一个环长为\(A_{i}\)，\(A\)可以看做一个\(N\)
色情图片识别
山水云的博客
11-18 4740
foreach ($posts as $k => $v) { $pics = $v->vPics->get_val(); $filename = $pics[0]->val; $arr = getimagesize($filename); list($width, $height, $type) = getimagesize($filename); switch
MATLAB中运用imagesc函数显示图像，怎么出来的图像与原图像比例不一样？
热门推荐
lcj_cjfykx的专栏
05-15 1万+
本人使用imagesc函数显示图像，因为该参数可以加colorbar功能，但是显示的图像和原图像比例不一样，不知道什么原因？而且，colorbar上的数字字体大小虽然可以改变但是一存储，改大的字体又还原了。难道是此功能的matlab中bug？若有哪位知道怎么解决该问题，不甚感激！ kk_2010 (站内联系TA) imagesc(im2); axis image%图释%保持图
BZOJ 1488 Luogu P4727 [HNOI2009]图的同构 (Burnside引理、组合计数)
weixin_30439131的博客
08-03 109
题目链接 (Luogu) https://www.luogu.org/problem/P4727 (BZOJ) https://www.lydsy.com/JudgeOnline/problem.php?id=1488 题解 Burnside引理经典题。 首先考虑一个\(O(n!\times poly(n))\)暴力: 枚举点的置换，然后计算在置换下保持不变的图的个数。 把置换拆成若干个...
BZOJ 1488 Luogu P4727 [HNOI2009]图的同构 (群论、Burnside引理、组合计数)
suncongbo's blog
08-03 114
题目链接 (Luogu) https://www.luogu.org/problem/P4727 (BZOJ) https://www.lydsy.com/JudgeOnline/problem.php?id=1488 题解 Burnside引理经典题。 首先考虑一个\(O(n!\times poly(n))\)暴力: 枚举点的置换，然后计算在置换下保持不变的图的个数。 把置换拆成若干个轮换。 (...
【bzoj1004】[HNOI2008]Cards Burnside引理+背包dp
weixin_30293079的博客
01-16 90
题目描述 用三种颜色染一个长度为 $n=Sr+Sb+Sg$ 序列，要求三种颜色分别有 $Sr,Sb,Sg$ 个。给出 $m$ 个置换，保证这 $m$ 个置换和置换 ${1,2,3,...,n\choose 1,2,3,...,n}$ 构成一个置换群，求置换后不同构的序列个数模 $p$ 。 $0\le Sr,Sb,Sg\le 20,0\le m\le 60,m+1\le p\le 10...
【BZOJ1488】[HNOI2009]图的同构计数
weixin_30692143的博客
03-19 75
题目链接 题意 求 n 个点的同构意义下不同的图的数量。\((n\leq 60)\) Sol \(Polya\) 定理的练手题。 我们这里先把边的存在与否变成对边进行黑白染色,白色代表不存在,这样就变成了一个对完全图中的边进行染色的问题,于是可以使用 对\(Burnside\)引理 进行优化后的 \(Polya\) 定理。 显然总的置换群大小是 \(n!\) 对应了每一种对边进行重新编号的方案,关...
即将退役选手最后的挣扎
weixin_30911451的博客
03-29 162
动态dp 原理 对于树上一类支持动态修改，动态查询dp结果的问题，可以用动态dp解决。这种dp状态不能太多，不然复杂度会有问题。 其核心思想就是，我们把转移看做一个线性变化，对于重链维护线段树保存矩阵乘的答案，然后对于每个点的轻儿子的贡献当做一个常数矩阵，我们要查某个点的信息就可以在当前这个点向下延伸到底的那条重链上询问。 然后考虑修改的时候，跳到 top 然后改父亲的轻儿子的常数矩阵，大部分题都...
群论——置换群
weixin_30501857的博客
03-14 220
最近研究了一下有关置换群的东西……群论这个东西博大精深,我也就大概知道一下群的概念(网上随处可见)……置换这个东西博大精深,我也就大概该了解了一下相关概念:·置换:我们所说的置换是指集合论中的置换,并不是组合数学中的置换,所以其概念就是一个集合从自身到自身的双射·轮换、对换见http://www.doc88.com/p-4794548487547.html·还有一些其他的操作https://wen...
数据结构与算法--【链表1】力扣练习 || 链表 / 移除链表元素
qq_44578108的博客
08-02 792
将线性表L=(a0,a1,……,an-1)中各元素分布在存储器的不同存储块，称为结点，通过地址或指针建立元素之间的联系。每一个结点由两部分组成：数据域和指针域。结点的data域存放数据元素ai，而next域是一个指针，指向ai的直接后继ai+1所在的结点。// 单链表int val;// 节点上存储的元素// 指向下一个节点的指针ListNode(int x) : val(x), next(NULL) {} // 节点的构造函数// 单链表int val;// 节点上存储的元素。
算法训练.
cxj0695的博客
08-01 358
这应该是个二分题，但是我用的是贪心暴力写的，思想就是循坏每次都让湿度最高的使用一次烘衣机，要是湿度最高的可以在自然条件下都能晒干就结束循环，这样内部我第一想法就是每次都排个降序；但是交上去发现后面数据超时了，也就是每次都排序太麻烦了，于是边思考到每次只需要把刚刚使用过烘衣机过后的湿度在插入进去，就只要一个for循环即可了这样就会快很多；随后也是AC了。
【Story】《程序员面试的“八股文”辩论：技术基础与实际能力的博弈》
EleganceJiaBao的博客
07-30 1200
“八股文”在程序员面试中扮演了重要的角色，但它并不是评价一个程序员能力的唯一标准。面试过程应结合理论知识、实际编程能力、项目经验和解决问题的能力等多个方面，综合评估候选人的综合素质。通过合理的面试考察方式，能够确保招聘到的程序员具备扎实的技术基础和实际工作能力，满足企业的技术需求。
C语言程序设计之链表篇2
最新发布
繁华落尽心殇
08-02 905
建立单向链表的一般步骤是：建立头指针 -> 建立第一个结点 -> 头指针指向第一个结点 -> 建立第二个结点 -> 头指针指向第二个结点。考查链表的数据结构，需利用指针变量才能实现，一个结点中应包含一个指针变量，用来存放下一个结点的地址。的功能是：将不带头结点的单向链表结点数据中的数据从小到大排序。的功能是：将不带头结点的单向链表逆置，即若原链表中从头自尾结点数据域依次为。的功能是 ：统计带头结点的单向链表中结点的个数，并存放在形参。在带头结点的单向链表中，查找数据域中值为。的结点，函数返回 0 值。
【前端面试3+1】20 css三栏布局6种实现方式、多行文本溢出怎么实现、token过期了怎么处理、【二叉树的中序遍历】
qq_61798306的博客
07-30 764
本文包含八股文：css三栏布局6种实现方式、多行文本溢出怎么实现、token过期了怎么处理，算法【二叉树的中序遍历】
常用的机器学习模型的不同调用方式
weixin_71458119的博客
08-01 961
不同机器学习方法的调用
剑指offer题解合集——Week7day1[滑动窗口的最大值]
qq_51931826的博客
08-01 442
给定一个数组和滑动窗口的大小，请找出所有滑动窗口里的最大值。例如，如果输入数组 [2,3,4,2,6,2,5,1]及滑动窗口的大小 3，那么一共存在 6个滑动窗口，它们的最大值分别为 [4,4,6,6,6,5]注意：数据保证 k大于 0，且 k小于等于数组长度。数据范围数组长度 [1,1000]样例输入：[2, 3, 4, 2, 6, 2, 5, 1] , k=3输出: [4, 4, 6, 6, 6, 5]
数据结构：哈希表与二叉搜索树详解
hong161688的博客
08-02 519
哈希表以其高效的查找、插入和删除操作而著称，这使得它在处理大量数据时特别有用。通过精心设计的哈希函数和冲突解决策略，哈希表能够在常数时间内完成这些基本操作（平均情况下），尽管在最坏情况下可能会退化到线性时间复杂度。优点快速查找：通过哈希函数直接定位数据，查找效率高。动态伸缩：支持动态扩容，可以适应不同规模的数据集。实现简单：相比于其他数据结构，哈希表的实现相对简单。缺点空间浪费：为了减少冲突，哈希表可能会分配比实际需要更多的空间。敏感于哈希函数：哈希函数的质量直接影响哈希表的性能。不支持有序遍历。
bzoj.org俄罗斯方块
05-06
bzoj.org是一个在线的编程题库，其中包含了很多经典算法问题和实际应用中的编程问题。其中包括了很多有趣的游戏算法问题，比如Puzzle，Tetris等。 而bzoj.org俄罗斯方块就是在该平台上的一个经典的俄罗斯方块游戏。...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值