传送门.

不妨设\(A(x)\)表示答案。

对于一个点,考虑它的三个子节点,直接卷起来是\(A(x)^3\),但是这样肯定会计重,因为我们要的是无序的子节点。

那么用burnside引理,枚举一个排列,一个环的选择要相同,如果环的大小是y,则对应\(A(x^y)\)。

最后可以得到:

\(A(x)=x{A(x)^3+3A(x^2)A(x)+2A(x^3)\over 6}+1\)

分治NTT可以解这个方程,不过因为有3次的,比较复杂,考虑用牛顿迭代:

\(F(A(x))=x{A(x)^3+3A(x^2)A(x)+2A(x^3)\over 6}+1-A(x)=0\)

\(new A(x)=A(x)-{F(A(x))\over F(A(x))'}\)

问题在于\(F(A(x))\)中有\(A(x^2)、A(x^3)、x\)这样的项,如何求导。

注意牛顿迭代的倍增中,已经求出了\(mod~x^{n/2}\)的答案,那么\(A(x^2)、A(x^3)、x\)是已知的,可以视作常数。

所以\(F(A(x))'=x{3A(x)^2 + 3A(x^2)\over 6} - 1\)。

Code:

#include<bits/stdc++.h>
#define fo(i, x, y) for(int i = x, B = y; i <= B; i ++)
#define ff(i, x, y) for(int i = x, B = y; i < B; i ++)
#define fd(i, x, y) for(int i = x, B = y; i >= B; i --)
#define ll long long
#define pp printf
#define hh pp("\n")
using namespace std; const int mo = 998244353; ll ksm(ll x, ll y) {
ll s = 1;
for(; y; y /= 2, x = x * x % mo)
if(y & 1) s = s * x % mo;
return s;
} typedef vector<ll> V;
#define pb push_back
#define si size()
#define re resize namespace ntt {
const int nm = 1 << 18;
int r[nm]; ll w[nm], a[nm];
void build() {
for(int n = 1; n < nm; n *= 2) {
ll v = ksm(3, (mo - 1) / 2 / n) ;
w[n] = 1;
ff(i, 1, n) w[n + i] = w[n + i - 1] * v % mo;
}
}
void dft(V &c, int f) {
int n = c.si;
ff(i, 0, n) a[i] = c[i];
ff(i, 0, n) {
r[i] = r[i / 2] / 2 + (i & 1) * n / 2;
if(i < r[i]) swap(a[i], a[r[i]]);
} ll b;
for(int i = 1; i < n; i *= 2) for(int j = 0; j < n; j += 2 * i) ff(k, 0, i)
b = a[i + j + k] * w[i + k], a[i + j + k] = (a[j + k] - b) % mo, a[j + k] = (a[j + k] + b) % mo;
if(f == -1) {
reverse(a + 1, a + n);
b = ksm(n, mo - 2);
ff(i, 0, n) a[i] = (a[i] + mo) * b % mo;
}
ff(i, 0, n) c[i] = a[i];
}
}
using ntt :: dft; V operator *(V a, V b) {
int n0 = a.si + b.si - 1, n = 1;
while(n < n0) n *= 2;
a.re(n); b.re(n);
dft(a, 1); dft(b, 1);
ff(i, 0, n) a[i] = a[i] * b[i] % mo;
dft(a, -1);
a.re(n0); return a;
} V operator +(V a, V b) {
a.re(max(a.si, b.si));
ff(i, 0, b.si) a[i] = (a[i] + b[i]) % mo;
return a;
} V operator -(V a, V b) {
a.re(max(a.si, b.si));
ff(i, 0, b.si) a[i] = (a[i] - b[i] + mo) % mo;
return a;
} V operator * (V a, int b) {
ff(i, 0, a.si) a[i] = a[i] * b % mo;
return a;
} V qni(V a) {
V b; b.resize(1); b[0] = ksm(a[0], mo - 2);
for(int n = 2; n < a.si * 2; n *= 2) {
V c = a; c.re(n); c.re(2 * n);
V d = b; b.re(2 * n); d.re(n);
dft(b, 1); dft(c, 1);
ff(i, 0, b.si) b[i] = b[i] * b[i] % mo * c[i] % mo;
dft(b, -1); b.re(n);
ff(i, 0, n) b[i] = (2 * d[i] - b[i] + mo) % mo;
}
b.re(a.si); return b;
} V yy(V a) {
a.insert(a.begin(), 0);
return a;
}
V stay(V a, int y, int n) {
V b; b.resize(n);
int mx = a.si - 1; mx = min(mx, (n - 1) / y);
fo(i, 0, mx) b[i * y] = a[i];
return b;
} const ll ni6 = ksm(6, mo - 2); V newton(int n0) {
V b; b.re(1); b[0] = 1;
for(int n = 2; n <= n0 * 2; n *= 2) {
V c = b * b * b + stay(b, 2, n) * b * 3 + stay(b, 3, n) * 2;
c = yy(c * ni6); c.re(n);
c[0] ++; c = c - b;
V d = b * b * 3 + stay(b, 2, n) * 3;
d = yy(d * ni6); d[0] --; d.re(n);
b = b - c * qni(d); b.re(n);
}
return b;
} int n; int main() {
ntt :: build();
scanf("%d", &n);
V a = newton(n);
pp("%lld\n", a[n]);
}

LOJ #6538. 烷基计数 加强版 加强版(生成函数,burnside引理,多项式牛顿迭代)的更多相关文章

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

    别问我为啥突然刷了道OI题,也别问我为啥花括号不换行了... 题目描述 求含 $n$ 个碳原子的本质不同的烷基数目模 $998244353$ 的结果.$1\le n\le 10^5$ . 题解 Bur ...

  2. [LOJ 6185]烷基计数

    Description 众所周知,大连 24 中是一所神奇的学校,在那里,化竞的同学很多都擅长写代码. 有一天,化学不及格的胡小兔向化竞巨佬晴岚请教化学题: “n 个碳原子的烷基共有多少种同分异构体? ...

  3. 等价类计数(Polya定理/Burnside引理)学习笔记

    参考:刘汝佳<算法竞赛入门经典训练指南> 感觉是非常远古的东西了,几乎从来没有看到过需要用这个的题,还是学一发以防翻车. 置换:排列的一一映射.置换乘法相当于函数复合.满足结合律,不满足交 ...

  4. BZOJ5119 生成树计数(prufer+生成函数+分治FFT+多项式exp)

    https://www.luogu.org/problemnew/solution/P4002 神树的题解写的很清楚了.稍微补充: 1.[x^i]ln(A(ax))=a^i[x^i]ln(A(x)), ...

  5. 洛谷 P4708 - 画画(Burnside 引理+组合数学)

    洛谷题面传送门 神仙题 %%%%%%%%%%%%%%%%%%%% 题解搬运人来了 首先看到本质不同(无标号)的图计数咱们可以想到 Burnside 引理,具体来说,我们枚举一个排列 \(p\),并统计 ...

  6. Luogu P5564 [Celeste-B]Say Goodbye (多项式、FFT、Burnside引理、组合计数)

    题目链接 https://www.luogu.org/problem/P5564 题解 这题最重要的一步是读明白题. 为了方便起见下面设环长可以是\(1\), 最后统计答案时去掉即可. 实际上就相当于 ...

  7. 【等价的穿越】Burnside引理&Pólya计数法

    Problem 起源: SGU 294 He's Circle 遗憾的是,被吃了. Poj有道类似的: Mission 一个长度为n(1≤n≤24)的环由0,1,2组成,求有多少本质不同的环. 实际上 ...

  8. 等价类计数:Burnside引理 & Polya定理

    提示: 本文并非严谨的数学分析,有很多地方是自己瞎口胡的,仅供参考.有错误请不吝指出 :p 1. 群 1.1 群的概念 群 \((S,\circ)\) 是一个元素集合 \(S\) 和一种二元运算 $ ...

  9. BZOJ 1488 Luogu P4727 [HNOI2009]图的同构 (Burnside引理、组合计数)

    题目链接 (Luogu) https://www.luogu.org/problem/P4727 (BZOJ) https://www.lydsy.com/JudgeOnline/problem.ph ...

随机推荐

  1. nc临时开启端口并监听

    port="6379 3306 27017 4505 4506 24007 24008 49152" #while true #do for i in $port do isexi ...

  2. leetcode-163周赛-1263-推箱子*

    题目描述: 自己的提交: class Solution: def minPushBox(self, grid: List[List[str]]) -> int: driction = [(0,1 ...

  3. js 邮箱,手机电话验证

    var mail=/^([a-zA-Z0-9_\.\-])+\@(([a-zA-Z0-9\-])+\.)+([a-zA-Z0-9]{2,4})+$/; var tel=/^(\d3,4|\d{3,4} ...

  4. 【代码工具】Orika JavaBean映射工具探秘

    转自:https://www.cnblogs.com/albert1024/articles/8434741.html Orika是一个简单.快速的JavaBean拷贝框架,Orika使用字节代码生成 ...

  5. luoguP3806 【模板】点分治1 [点分治]

    题目背景 感谢hzwer的点分治互测. 题目描述 给定一棵有n个点的树 询问树上距离为k的点对是否存在. 输入输出格式 输入格式: n,m 接下来n-1条边a,b,c描述a到b有一条长度为c的路径 接 ...

  6. kafka为什么快?

    为什么Kafka那么快 原创 2016-07-18 fireflyc 写程序的康德 网上有很多Kafka的测试文章,测试结果通常都是"吊打"其他MQ.感慨它的牛B之余我觉得必要仔细 ...

  7. 一:unittest框架配合selenium工具之CSS_selector定位。

    做了自动化测试这么久了,一直没有梳理到元素定位这一块的内容,其重要性不言而喻.趁着周末有时间,梳理一下. 1,通过id定位 driver.find_element_by_css_selector(&q ...

  8. fatal error C1076: compiler limit : internal heap limit reached; use /Zm to specify a higher limit

    最近想用一下Xtreme ToolkitPro 界面库,安装后用VC6根据向导 产生一个工程,编译时出现如下的错误: fatal error C1076: compiler limit : inter ...

  9. html01. <!DOCTYPE html>

    解释 在HTML文档初,往往会有这么一句话<!DOCTYPE html>,它是html5标准网页声明,全称为Document Type HyperText Mark-up Language ...

  10. PHP-FPM 远程代码执行漏洞(CVE-2019-11043)复现-含EXP

    搭建容器 安装golang 利用程序 https://github.com/neex/phuip-fpizdam 安装git Cobra包安装 go get -v github.com/spf13/c ...