【loj6538】烷基计数 加强版 加强版 Burnside引理+多项式牛顿迭代

别问我为啥突然刷了道OI题,也别问我为啥花括号不换行了...

题目描述

求含 $n$ 个碳原子的本质不同的烷基数目模 $998244353$ 的结果。$1\le n\le 10^5$ 。


题解

Burnside引理+多项式牛顿迭代

不考虑同构的话,很容易想到dp方程 $\begin{cases}f_0=1\\f_i=\sum\limits_{j+k+l+1=i}f_jf_kf_l\end{cases}$ 。

考虑同构,可以通过容斥原理,大力讨论一下容斥系数。一个更简单的方法是考虑Burnside引理,即:等价类的数目等于每个置换下不动点数目的平均值。

  • 对于置换 $(1,2,3)$ ,所有组合都是不动点;
  • 对于置换 $(1,3,2)$ 、$(2,1,3)$ 和 $(3,2,1)$ ,“$2+1$” 的组合是不动点;
  • 对于置换 $(2,3,1)$ 和 $(3,1,2)$ ,只有 “$3$” 的组合是不动点。

于是新的dp方程为 $\begin{cases}f_0=1\\f_i=\frac{\sum\limits_{j+k+l+1=i}f_jf_kf_l+\sum\limits_{2j+k+1=i}3f_jf_k+\sum\limits_{3j+1=i}2f_j}6\end{cases}$ 。

$n$ 这么大肯定不能直接dp,考虑多项式解法,则dp方程的多项式形式为 $F(x)=x\cdot\frac{F^3(x)+3F(x)F(x^2)+2F(x^3)}6+1$ 。

由于出现了三次方和 $F(x^2)$ 、$F(x^3)$ 项,因此这个方程难以直接解出。

考虑牛顿迭代,则当我们已知 $F(x)\mod x^n$ 时,$F(x^2)\mod x^{2n}$ 和 $F(x^3)?\text{mod}?x^{2n}$ 就已经是已知量,在迭代时可以当作常量处理。

记 $S(x)=F(x^2)$ ,$C(x)=F(x^3)$ ,则我们要迭代的方程就是 $G(T(x))=x\cdot\frac{T^3(x)+3S(x)T(x)+2C(x)}6-T(x)+1$ 的零点 $G(F(x))=0$ 。

又因为 $G‘(T(x))=x\cdot\frac{3T^2(x)+3S(x)}6-1$ ,代入牛顿迭代公式 $F(x)=F_0(x)-\frac{G(F_0(x))}{G‘(F_0(x))}$ 中可得 $F(x)=F_0(x)-\frac{x(F_0^3(x)+3S(x)F_0(x)+2C(x))-6F_0(x)+6}{3F_0^2(x)+3S(x)-6}$ ,其中 $F_0(x)$ 是上次迭代所得的多项式。

最后的答案就是 $F(x)[n]$ 。

时间复杂度 $O(n\log n)$ 。

#include <cstdio>
#include <algorithm>
#define N 262155
#define mod 998244353
using namespace std;
typedef long long ll;
ll F[N];
inline ll qpow(ll x , ll y) {
    ll ans = 1;
    while(y) {
        if(y & 1) ans = ans * x % mod;
        x = x * x % mod , y >>= 1;
    }
    return ans;
}
inline void ntt(ll *A , int n , ll flag) {
    int i , j , k;
    for(k = i = 0 ; i < n ; i ++ ) {
        if(i < k) swap(A[i] , A[k]);
        for(j = (n >> 1) ; (k ^= j) < j ; j >>= 1);
    }
    for(k = 2 ; k <= n ; k <<= 1) {
        ll wn = qpow(3 , (mod - 1) / k * flag);
        for(i = 0 ; i < n ; i += k) {
            ll w = 1 , t;
            for(j = i ; j < i + (k >> 1) ; j ++ , w = w * wn % mod)
                t = w * A[j + (k >> 1)] % mod , A[j + (k >> 1)] = (A[j] - t + mod) % mod , A[j] = (A[j] + t) % mod;
        }
    }
    if(flag == mod - 2) {
        ll t = qpow(n , flag);
        for(i = 0 ; i < n ; i ++ ) A[i] = A[i] * t % mod;
    }
}
inline void inv(ll *A , ll *B , int n) {
    static ll T[N];
    int i , j;
    for(i = 0 ; i < (n << 1) ; i ++ ) B[i] = 0;
    B[0] = qpow(A[0] , mod - 2);
    for(i = 2 ; i <= n ; i <<= 1) {
        for(j = 0 ; j < i ; j ++ ) T[j] = A[j] , T[j + i] = 0;
        ntt(T , i << 1 , 1) , ntt(B , i << 1 , 1);
        for(j = 0 ; j < (i << 1) ; j ++ ) B[j] = B[j] * (2 - T[j] * B[j] % mod + mod) % mod;
        ntt(B , i << 1 , mod - 2);
        for(j = i ; j < (i << 1) ; j ++ ) B[j] = 0;
    }
}
inline void solve(int n) {
    static ll G[N] , H[N] , S[N] , C[N] , T[N];
    int i , j;
    F[0] = 1;
    for(i = 2 ; i <= n ; i <<= 1) {
        for(j = 0 ; j < i ; j ++ ) T[j] = F[j] , S[j] = C[j] = T[j + i] = S[j + i] = C[j + i] = 0;
        for(j = 0 ; j < i ; j += 2) S[j] = F[j / 2];
        for(j = 0 ; j < i ; j += 3) C[j] = F[j / 3];
        ntt(T , i << 1 , 1) , ntt(S , i << 1 , 1);
        for(j = 0 ; j < (i << 1) ; j ++ ) G[j] = T[j] * (T[j] * T[j] % mod + 3 * S[j]) % mod , H[j] = 3 * (T[j] * T[j] + S[j]) % mod;
        ntt(G , i << 1 , mod - 2) , ntt(H , i << 1 , mod - 2);
        for(j = i ; j < (i << 1) ; j ++ ) G[j] = H[j] = 0;
        for(j = i - 1 ; j ; j -- ) G[j] = ((G[j - 1] + 2 * C[j - 1] - 6 * F[j]) % mod + mod) % mod , H[j] = H[j - 1];
        G[0] = 0 , H[0] = mod - 6;
        ntt(G , i << 1 , 1) , inv(H , T , i) , ntt(T , i << 1 , 1);
        for(j = 0 ; j < (i << 1) ; j ++ ) G[j] = G[j] * T[j] % mod;
        ntt(G , i << 1 , mod - 2);
        for(j = 0 ; j < i ; j ++ ) F[j] = (F[j] - G[j] + mod) % mod;
    }
}
int main() {
    int n , len = 1;
    scanf("%d" , &n);
    while(len <= n) len <<= 1;
    solve(len);
    printf("%lld\n" , F[n]);
    return 0;
}

原文地址:https://www.cnblogs.com/GXZlegend/p/12056937.html

时间: 2024-10-15 23:46:02

【loj6538】烷基计数 加强版 加强版 Burnside引理+多项式牛顿迭代的相关文章

BZOJ_[HNOI2008]_Cards_(置换+Burnside引理+乘法逆元+费马小定理+快速幂)

描述 http://www.lydsy.com/JudgeOnline/problem.php?id=1004 共n个卡片,染成r,b,g三种颜色,每种颜色的个数有规定.给出一些置换,可以由置换得到的染色方案视为等价的,求等价类计数. 分析 给出置换求等价类计数,用Burnside引理:等价类计数=(每一个置换不动点的和)/置换数.(不知道的建议去看白书) 其中不动点是指一个染色方案经过置换以后染色与之前完全相同. 1.求不动点个数. 不动点的话同一个循环内的每一个点的颜色必须相同(否则不同颜色

Burnside引理与Polya定理

Burnside引理与Polya定理 Burnside引理与Polya定理是有关组合数学的两条十分重要的定理(引理),但是网上的一些资料大多晦涩难懂或者与实际并不相关联,因此在这里做一些浅显的解读,希望通过此文章可以让这两条定理(引理)能够发挥其作用. PS:引理与定理的区别: Ψ引理是数学中为了取得某个更好的定理而作为步骤被证明的命题,其意义并不在于自身被证明,而在于为达成最终定理作出贡献. Ψ一个引理可用于证明多个定理.数学中存在很多著名的引理,这些引理可能对很多问题的解决有帮助.例如欧几里

hdu 5868 2016 ACM/ICPC Asia Regional Dalian Online 1001 (burnside引理 polya定理)

Different Circle Permutation Time Limit: 3000/1500 MS (Java/Others)    Memory Limit: 262144/262144 K (Java/Others)Total Submission(s): 208    Accepted Submission(s): 101 Problem Description You may not know this but it's a fact that Xinghai Square is

bzoj1004: [HNOI2008]Cards Burnside引理+01背包

三维01背包算出在每一个置换下不变的染色方案数,Burnside引理计算答案. PS:数据太水所以只算恒等置换也是可以过的. #include<bits/stdc++.h> using namespace std; int n,m,p,x,y,z; bool u[61]; int f[21][21][21],s[61],v[61]; int power(int u,int v){ int d=1; for(;v;v>>=1){ if(v&1) d=d*u%p; u=u*u%

Polya定理,Burnside引理(转)

设G是一个集合,*是G上的二元运算,如果(G,*)满足下面的条件: 封闭性:对于任何a,b∈G,有a*b∈G; 结合律:对任何a,b,c∈G有(a*b)*c=a*(b*c); 单位元:存在e∈G,使得对所有的a∈G,都有a*e=e*a=a; 逆元:对于每个元素a∈G,存在x∈G,使得a*x=x*a=e,这个时候记x为a-1,称为a的逆元,那么则称(G,*)为一个群. 例:G={0,1,2,3,4....n-1}那么它在mod n加法下是一个群. 群元素的个数有限,称为有限群,且其中元素的个数称为

bzoj1004: [HNOI2008]Cards(burnside引理+DP)

题目大意:3种颜色,每种染si个,有m个置换,求所有本质不同的染色方案数. 置换群的burnside引理,还有个Pólya过几天再看看... burnside引理:有m个置换k种颜色,所有本质不同的染色方案数就是每种置换的不变元素的个数的平均数. 求每种置换的不变元素的个数用背包解决.因为置换之后元素不变,所以对于每个循环节我们要染一个颜色,于是先处理出循环节作为背包中的"物体",然后一个三维背包解决.f[i][j][k]的i j k表示三种颜色分别还可以染多少次. 除m%p用费马小定

[bzoj1004][HNOI2008][Cards] (置换群+Burnside引理+动态规划)

Description 小春现在很清闲,面对书桌上的N张牌,他决定给每张染色,目前小春只有3种颜色:红色,蓝色,绿色.他询问Sun有多少种染色方案,Sun很快就给出了答案.进一步,小春要求染出Sr张红色,Sb张蓝色,Sg张绝色.他又询问有多少种方案,Sun想了一下,又给出了正确答案. 最后小春发明了M种不同的洗牌法,这里他又问Sun有多少种不同的染色方案.两种染色方法相同当且仅当其中一种可以通过任意的洗牌法(即可以使用多种洗牌法,而每种方法可以使用多次)洗成另一种.Sun发现这个问题有点难度,决

BZOJ 1004 HNOI 2008 Cards Burnside引理

题目大意:给出一个置换群,求有多少种本质不同的染色方案. 思路:Burnside引理:置换群的等价类数目=所有置换的不动点数目的平均值. 有了这个引理,我们只需要求出所有不动点的数目求一个平均值就可以的到等价类的数目了. 要使一种染色的方案在一种置换的意义下是不动点,需要让这个置换的每个循环节中的颜色都相同.先求出所有置换的循环,然后用一个背包就可以初解了.最后乘法逆元搞一下除法. CODE: #define _CRT_SECURE_NO_DEPRECATE #include <cstdio>

BZOJ 1004 HNOI2008 Cards Burnside引理

题目大意:给定n张卡牌和m个置换,求等价类个数 数据保证这m个置换加上自身置换后构成一个置换群 BZOJ坑爹0.0 这么重要的条件不给出来尼玛怎么做 Burnside引理--昨晚为了做这题硬啃了一晚上白书0.0 都快啃吐了0.0 Burnside引理:一个置换群下的等价类个数等于所有置换的不动点个数的平均值 没有接触过群论的建议去啃白书-- 网上的东西看不懂的 最后那个除法要用乘法逆元 我懒得写EXGCD写了费马小定理0.0 #include<cstdio> #include<cstri