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

题目大意

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

题解

这个题是学习了Polya定理和群论以后的练手题，但是推了好久并没有推出来。。。。真的是太难辣。。。
首先我先说一下我错误的想法：
很容易就把这个题转化成了给(K_n)的完全图上的边进行二着色的问题，然后，由于在组合数学课程中经常接触到多边形着色，所以我就把这个题错误的转化成了在一个正(frac{n(n-1)}{2})边形的顶点上进行二着色的问题。然而对于n=1,2,3这种方法都是可行的，但是到了n=4的情况，这种方法就不可行了。我仔细观察了一下，发现这个转化不符合满足纯粹性和完备性。。。
然后就说一下正解吧。
首先我们考虑n = 4的情况，对于(K_4)进行二着色，我们很容易发现，由于图是可以任意扭转的，所以它的置换群实际上是一个对称群！
那么对于点的每一个置换我们要计算对应的边的置换。
在一个置换中，考察一条边，如果这条边的两个节点位于相同的循环中，那么我们可以得出边的循环个数是点的循环个数的一半。
如果这条边的两个节点位于不同的循环中，那么我们画一画图就可以知道如果点的循环个数分别是a, b,那么边的循环个数就是gcd(a,b)。
根据这样的方法，我们就可以把点的置换转化为边的置换了。
下面的任务就是要枚举置换。
如果直接暴力，复杂度很高。
我们考虑这样的枚举（回溯）方法：
依次考虑每一种阶的循环的个数，然后暴力dfs即可。
现在假设我们已经枚举好了一个置换，那么这种置换的个数根据一些基本的排列组合知识，可以知道是：

[frac{n!}{prod_{i = 1}^{cnt} Val_i * Num_i !} ]

稍微解释一下这个式子。除以(Val_i)是因为圆形排列，除以(Num_i)是因为同阶循环的重复排列。

根据Polya定理，等价类的个数就是：

[l = frac{1}{N!} * sum 2^m ]

参考题解

事实上，这个题还有一个变态的做法：
就是上OEIS上查询通项公式。。。。

代码

#include <bits/stdc++.h>
const int mod = 997;
const int maxn = 1010;
using namespace std;
int n, cnt, ans;
int two[maxn], factor[maxn], val[maxn], num[maxn];
int pow(int n, int m) {
    int ans = 1;
    int b = m;
    while(b) {
        if(b & 1) ans = (ans * n) % mod;
        b >>= 1;
        n = (n*n) % mod;
    }
    return ans;
}
int inv(int n) {
    return pow(n, mod-2);
}
int gcd(int a, int b) {
    if(b == 0) return a;
    else return gcd(b, a%b) % mod;
}
void init() {
    factor[0] = factor[1] = two[0] = 1;
    for(int i = 2; i <= 1000; i++) {
        factor[i] = ((i % mod) * factor[i-1]) % mod;
    }
    for(int i = 1; i <= 1000; i++) {
        two[i] = (two[i-1] * 2) % mod;
    }
}
void dfs(int now_num, int left) {
    if(left == 0) {
        int sum1 = 0, sum2 = 1;
        //sum1:这一种置换的循环个数
        //sum2:这一种置换的个数
        for(int i = 1; i <= cnt; i++) {
            sum1 += (num[i] * (num[i] - 1) / 2 * val[i]) + (val[i]/2 * num[i]); 
            //前一部分:对于同一种循环中的不同循环的边的处理
            for(int j = i + 1; j<= cnt; j++) {
                sum1 += num[i] * num[j] * gcd(val[i], val[j]);
            }
        }
        for(int i = 1; i <= cnt; i++) {
            sum2 = (sum2 * pow(val[i], num[i])%mod*factor[num[i]])%mod;
        }
        sum2 = inv(sum2) * factor[n] % mod;
        ans = (ans + pow(2, sum1) * sum2 % mod) % mod;
    }
    if(now_num > left) return;
    dfs(now_num+1, left);
    //这里的dfs放到外面可以降低常数：
    //如果放在for循环里面，那么num数组中会多出许多0
    //浪费时间。
    for(int i = 1; i * now_num <= left; i++) {
        val[++cnt] = now_num, num[cnt] = i;
        dfs(now_num+1, left - i * now_num);
        cnt--; //回溯
    }
}

int main() {
    scanf("%d", &n);
    init();
    dfs(1, n);
    ans = (ans * inv(factor[n])) % mod;
    printf("%d", ans);
    return 0;
}