#3456. 城市规划(生成函数,多项式求逆)

#3456. 城市规划

fnf_nfnnnn个点的的点的简单无向连通图数目,gng_ngnnnn个点的简单无向图个数(不要求联通)。

对于gng_ngn显然有gn=2n(n−1)2g_n = 2 ^{\frac{n(n - 1)}{2}}gn=22n(n1),共有n(n+1)2\frac{n(n + 1)}{2}2n(n+1)条边,然后每条边可选可不选。

我们枚举111所在的点的连通块可得:
gn=∑i=1nCn−1i−1fign−i2(2n)=∑i=1n(i−1n−1)fi2(2n−i)2(2n)=∑i=1n(n−1)!(i−1)!(n−i)!fi2(2n−i)2(2n)(n−1)!=∑i=1nfi(i−1)!2(2n−i)(n−i)!设G(x)=∑n=1∞2(2n)(n−1)!xnF(x)=∑n=1∞fn(n−1)!H(x)=∑n=0∞2(2n)n!G(x)=F(x)H(x)F(x)=G(x)H−1(x)构造多项式,多项式求逆,把[xn]项系数乘上(n−1)!即是答案g_n = \sum\limits_{i = 1} ^{n}C_{n - 1} ^{i - 1}f_ig_{n - i}\\ 2 ^{(_2 ^ n)} = \sum_{i = 1} ^{n}(_{i - 1} ^{n - 1})f_i 2 ^{(_2 ^{n - i})}\\ 2 ^{(_2 ^ n)} = \sum_{i = 1} ^{n} \frac{(n - 1)!}{(i - 1)!(n - i)!} f_i 2 ^{(_2 ^{n - i})}\\ \frac{2 ^{(_2 ^ n)}}{(n - 1)!} = \sum_{i = 1} ^{n} \frac{f_i}{(i - 1)!} \frac{2^{(_2 ^{n - i})}}{(n - i)!}\\ 设G(x) = \sum_{n = 1} ^{\infty} \frac{2 ^{(_2 ^n)}}{(n - 1)!} x ^ n\\ F(x) = \sum_{n = 1} ^{\infty} \frac{f_n}{(n - 1)!}\\ H(x) = \sum_{n = 0} ^{\infty} \frac{2 ^{(_2 ^n)}}{n!}\\ G(x) = F(x) H(x)\\ F(x) = G(x)H^{-1}(x)\\ 构造多项式,多项式求逆,把[x ^ n]项系数乘上(n - 1)!即是答案\\ gn=i=1nCn1i1figni2(2n)=i=1n(i1n1)fi2(2ni)2(2n)=i=1n(i1)!(ni)!(n1)!fi2(2ni)(n1)!2(2n)=i=1n(i1)!fi(ni)!2(2ni)G(x)=n=1(n1)!2(2n)xnF(x)=n=1(n1)!fnH(x)=n=0n!2(2n)G(x)=F(x)H(x)F(x)=G(x)H1(x)[xn](n1)!

#include <bits/stdc++.h>using namespace std;const int mod = 1004535809, inv2 = mod + 1 >> 1;namespace Quadratic_residue {struct Complex {int r, i;Complex(int _r = 0, int _i = 0) : r(_r), i(_i) {}};int I2;Complex operator * (const Complex &a, Complex &b) {return Complex((1ll * a.r * b.r % mod  + 1ll * a.i * b.i % mod * I2 % mod) % mod, (1ll * a.r * b.i % mod + 1ll * a.i * b.r % mod) % mod);}Complex quick_pow(Complex a, int n) {Complex ans = Complex(1, 0);while (n) {if (n & 1) {ans = ans * a;}a = a * a;n >>= 1;}return ans;}int get_residue(int n) {mt19937 e(233);if (n == 0) {return 0;}if(quick_pow(n, (mod - 1) >> 1).r == mod - 1) {return -1;}uniform_int_distribution<int> r(0, mod - 1);int a = r(e);while(quick_pow((1ll * a * a % mod - n + mod) % mod, (mod - 1) >> 1).r == 1) {a = r(e);}I2 = (1ll * a * a % mod - n + mod) % mod;int x = quick_pow(Complex(a, 1), (mod + 1) >> 1).r, y = mod - x;if(x > y) swap(x, y);return x;}
}const int N = 1e6 + 10;int r[N], inv[N], b[N], c[N], d[N], e[N], t[N];int quick_pow(int a, int n) {int ans = 1;while (n) {if (n & 1) {ans = 1ll * a * ans % mod;}a = 1ll * a * a % mod;n >>= 1;}return ans;
}void get_r(int lim) {for (int i = 0; i < lim; i++) {r[i] = (i & 1) * (lim >> 1) + (r[i >> 1] >> 1);}
}void get_inv(int n) {inv[1] = 1;for (int i = 2; i <= n; i++) {inv[i] = 1ll * (mod - mod / i) * inv[mod % i] % mod;}
}void NTT(int *f, int lim, int rev) {for (int i = 0; i < lim; i++) {if (i < r[i]) {swap(f[i], f[r[i]]);}}for (int mid = 1; mid < lim; mid <<= 1) {int wn = quick_pow(3, (mod - 1) / (mid << 1));for (int len = mid << 1, cur = 0; cur < lim; cur += len) {int w = 1;for (int k = 0; k < mid; k++, w = 1ll * w * wn % mod) {int x = f[cur + k], y = 1ll * w * f[cur + mid + k] % mod;f[cur + k] = (x + y) % mod, f[cur + mid + k] = (x - y + mod) % mod;}}}if (rev == -1) {int inv = quick_pow(lim, mod - 2);reverse(f + 1, f + lim);for (int i = 0; i < lim; i++) {f[i] = 1ll * f[i] * inv % mod;}}
}void polyinv(int *f, int *g, int n) {if (n == 1) {g[0] = quick_pow(f[0], mod - 2);return ;}polyinv(f, g, n + 1 >> 1);for (int i = 0; i < n; i++) {t[i] = f[i];}int lim = 1;while (lim < 2 * n) {lim <<= 1;}get_r(lim);NTT(t, lim, 1);NTT(g, lim, 1);for (int i = 0; i < lim; i++) {int cur = (2 - 1ll * g[i] * t[i] % mod + mod) % mod;g[i] = 1ll * g[i] * cur % mod;t[i] = 0;}NTT(g, lim, -1);for (int i = n; i < lim; i++) {g[i] = 0;}
}void polysqrt(int *f, int *g, int n) {if (n == 1) {g[0] = Quadratic_residue::get_residue(f[0]);return ;}polysqrt(f, g, n + 1 >> 1);polyinv(g, b, n);int lim = 1;while (lim < 2 * n) {lim <<= 1;}get_r(lim);for (int i = 0; i < n; i++) {t[i] = f[i];}NTT(g, lim, 1);NTT(b, lim, 1);NTT(t, lim, 1);for (int i = 0; i < lim; i++) {g[i] = (1ll * inv2 * g[i] % mod + 1ll * inv2 * b[i] % mod * t[i] % mod) % mod;b[i] = t[i] = 0;}NTT(g, lim, -1);for (int i = n; i < lim; i++) {g[i] = 0;}
}void derivative(int *a, int *b, int n) {for (int i = 0; i < n; i++) {b[i] = 1ll * a[i + 1] * (i + 1) % mod;}
}void integrate(int *a, int n) {for (int i = n - 1; i >= 1; i--) {a[i] = 1ll * a[i - 1] * inv[i] % mod;}a[0] = 0;
}void polyln(int *f, int *g, int n) {polyinv(f, b, n);derivative(f, g, n);int lim = 1;while (lim < 2 * n) {lim <<= 1;}get_r(lim);NTT(g, lim, 1);NTT(b, lim, 1);for (int i = 0; i < lim; i++) {g[i] = 1ll * g[i] * b[i] % mod;b[i] = 0;}NTT(g, lim, -1);for (int i = n; i < lim; i++) {g[i] = 0;}integrate(g, n);
}void polyexp(int *f, int *g, int n) {if (n == 1) {g[0] = 1;return ;}polyexp(f, g, n + 1 >> 1);int lim = 1;while (lim < 2 * n) {lim <<= 1;}polyln(g, d, n);for (int i = 0; i < n; i++) {t[i] = (f[i] - d[i] + mod) % mod;}t[0] = (t[0] + 1) % mod;get_r(lim);NTT(g, lim, 1);NTT(t, lim, 1);for (int i = 0; i < lim; i++) {g[i] = 1ll * g[i] * t[i] % mod;t[i] = d[i] =  0;}NTT(g, lim, -1);for (int i = n; i < lim; i++) {g[i] = 0;}
}/*b存放多项式逆,c存放多项式开根,d存放多项式对数ln,e存放多项式指数exp,t作为中间转移数组,如果要用到polyln,得提前调用get_inv(n)先预先得到我们想要得到的逆元范围。
*/int h[N], g[N], fac[N], ifac[N], n;void init() {fac[0] = 1;for (int i = 1; i < N; i++) {fac[i] = 1ll * i * fac[i - 1] % mod;}ifac[N - 1] = quick_pow(fac[N - 1], mod - 2);for (int i = N - 2; i >= 0; i--) {ifac[i] = 1ll * ifac[i + 1] * (i + 1) % mod;}
}int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);// ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);scanf("%d", &n);init();for (int i = 1; i <= n; i++) {g[i] = 1ll * quick_pow(2, 1ll * i * (i - 1) / 2 % (mod - 1)) * ifac[i - 1] % mod;}for (int i = 0; i <= n; i++) {h[i] = 1ll * quick_pow(2, 1ll * i * (i - 1) / 2 % (mod - 1)) * ifac[i] % mod;}polyinv(h, b, n + 1);for (int i = 0; i <= n; i++) {h[i] = b[i];b[i] = 0;}int lim = 1;while (lim <= 2 * n) {lim <<= 1;}get_r(lim);NTT(g, lim, 1);NTT(h, lim, 1);for (int i = 0; i < lim; i++) {g[i] = 1ll * g[i] * h[i] % mod;h[i] = 0;}NTT(g, lim, -1);printf("%d\n", 1ll * g[n] * fac[n - 1] % mod);return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/313937.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【盛派内部分享资料】本期主题:使用JavaScript HTML CSS构建跨平台桌面应用

先上录制视频友情提示&#xff1a;如果公众号内视频无法显示高清视频&#xff0c;您也可以在小程序内观看高清视频&#xff0c;点击下方按钮观看&#xff1a;【盛派内部分享资料】本期主题&#xff1a;使用JavaScript HTML CSS构建跨平台桌面应用本次活动由盛派技术人员伏允坤主…

负载均衡及负载均衡器

负载均衡在分布式系统中&#xff0c;负载均衡&#xff08;load balancing&#xff09;是一种有效的将网络请求分配到多个服务器的过程。通过将负载进行负载均衡&#xff0c;可以有效地改进系统响应时间&#xff0c;提高系统的可用性。随着系统变的愈发复杂&#xff0c;用户增多…

Sequence II (HDU 5919)(主席树)

Sequence II 题目大意是有mmm次询问&#xff0c;每次询问一段区间[l,r][l, r][l,r]&#xff0c;从左到右&#xff0c;如果这个数是在这个区间第一次出现&#xff0c;则记录下其下标&#xff0c; 我们会得到一个新的数组&#xff0c;要求这个数组的中位数是什么。 考虑使用主…

VS Code 摸鱼插件开发小记

插件地址&#xff1a;https://marketplace.visualstudio.com/items?itemNamebugbreeder.vscode-readhub原文链接&#xff1a;https://github.com/alex-yh99/vscode-readhub/blob/master/docs/develop-note.md一、关于插件在 VSCode 中快捷查阅科技动态、开发者资讯等 Readhub …

帕斯卡三角形杨辉三角

给定一个非负整数 numRows&#xff0c;生成杨辉三角的前 numRows 行。 在杨辉三角中&#xff0c;每个数是它左上方和右上方的数的和。 示例: 输入: 5 输出: [[1],[1,1],[1,2,1],[1,3,3,1],[1,4,6,4,1] ]解法 思路1&#xff1a;循环遍历 public List<List<Integer>…

ICPC 南昌现场赛 K:Tree(dsu on tree + 动态开点线段树)

Tree 让我们找满足一下五个条件的(x,y(x, y(x,y)点对有多少&#xff1a; x≠yx \neq yx​yxxx不是yyy的祖先yyy不是xxx的祖先dis(x,y)≤kdis(x, y)\leq kdis(x,y)≤kzzz是x,yx, yx,y的最近公共祖先&#xff0c;valuexvaluey2valuezvalue_x value_y 2value_zvaluex​valuey…

缺失数字

题目描述 给定一个包含 0, 1, 2, …, n 中 n 个数的序列&#xff0c;找出 0 … n 中没有出现在序列中的那个数。 示例 1: 输入: [3,0,1] 输出: 2示例 2: 输入: [9,6,4,2,3,5,7,0,1] 输出: 8解法 思路1&#xff1a;使用额外数组存储0~n&#xff0c;遍历数组&#xff0c;缺失…

9月数据库排行:Microsoft SQL Server分数罕见下滑

DB-Engines 数据库流行度排行榜 9 月更新已发布&#xff0c;排名前二十如下&#xff1a;明显能看出来&#xff0c;整体排名和上个月的保持一致。而上个月出现分数下跌的 PostgreSQL 和 MongoDB 终于“回归正道”&#xff0c;分数有了稍微的上升。与上个月相比&#xff0c;Postg…

P4175 [CTSC2008]网络管理(树状数组套主席树)

P4175 [CTSC2008]网络管理 每次询问两点间的第kkk大&#xff0c;考虑建立nnn棵主席树&#xff0c;每棵主席树记录的是从根节点到当前节点的状态&#xff0c; 由于主席树是一种类似前缀和的数据结构&#xff0c;可以使用差分的方式得到(u,v)(u, v)(u,v)点对之间的主席树为uv−…

Exceptionless使用介绍

一.Exceptionless介绍Exceptionless专注于.net平台提供实时错误和日志报告。主要包括&#xff1a;错误通知、智能分组异常、详细错误报告堆栈跟踪、支持离线、UI查看重要错误和确定优先级、仪表板上的统计信息和趋势、对异常标记为已修复&#xff0c;监视回归、将事件标记为关键…

BIO、NIO 和 AIO

BIO、NIO 和 AIO的区别 https://www.cnblogs.com/barrywxx/p/8430790.html 三种 IO 的用法与原理 netty http://ifeve.com/netty5-user-guide/

#279. [SYZOI Round1] 滑稽♂树(树状数组套主席树)

#279. [SYZOI Round1] 滑稽♂树 子树上的问题&#xff0c;考虑dfsdfsdfs序&#xff0c;第kkk大&#xff0c;可以用主席树嘛&#xff0c;支持修改&#xff0c;那就树状数组上套主席树&#xff0c;参考P4175 [CTSC2008]网络管理&#xff08;树状数组套主席树&#xff09; #incl…

Net Core DocXCore 实现word模板导出

实际工作中&#xff0c;往往有这样的需求&#xff0c;需要导出word&#xff0c;还有各种各样的样式&#xff0c;于是有了word模板导出。实现以下几个需求&#xff1a;1、表单导出2、表格导出3、表单表格混合导出4、实际用例测试解决方案&#xff1a;实现是基于NET Core 2.1 &am…

主席树有关的一些题目(持续更新)

主席树 模板 P3919 【模板】可持久化线段树 1&#xff08;可持久化数组&#xff09; #include <bits/stdc.h>using namespace std;const int N 1e6 10;int root[N], n, m;int ls[N * 25], rs[N * 25], value[N * 25], tot;void build(int &rt, int l, int r) {r…

密钥分离,.Net程序猿不再背锅

引言互联网每隔一段时间就会爆出 【某程序猿在代码托管平台上传了公司机密配置信息&#xff0c;导致公司核心数据被黑客获取或修改】&#xff0c; 一茬又一茬背锅侠层出不穷。拯救自我软件工程理论早以加粗字体给出 经典原则&#xff1a;Never store production passwords or o…

Java ASM与Javassit

操作java字节码的工具有两个比较流行&#xff0c;一个是ASM&#xff0c;一个是Javassit 。 ASM &#xff1a;直接操作字节码指令&#xff0c;执行效率高&#xff0c;要是使用者掌握Java类字节码文件格式及指令&#xff0c;对使用者的要求比较高。 Javassit 提供了更高级的API…

Chrome让人失望,是时候转到Firefox或Edge?

我们曾报道过谷歌正在移除 Chrome 的关闭其他选项卡&#xff0c;这引起很多用户的讨论&#xff0c;很多用户表现出对 Chrome 浏览器的失望。而在国外&#xff0c;也有很多人认为谷歌这项修改使 Chrome 变得越来越臃肿。还有前段时间&#xff0c;谷歌在 Chrome 76 中隐藏了 “ht…

ICPC 徐州 H Yuuki and a problem (树状数组套主席树)

Yuuki and a problem 先不管第一问的修改操作&#xff0c;考虑如何达到第二问的查询操作&#xff0c; 题目要我们给出一个区间[l,r][l, r][l,r]中&#xff0c;不能通过权值得到的最小的数字是什么&#xff0c; 假设我们已经可以得到[1,x][1, x][1,x]之间的数了&#xff0c;且…

Java序列化、反序列化

目录 1&#xff1a;什么是序列化、反序列化&#xff1f; 2&#xff1a;序列化的用途&#xff1f; 3&#xff1a;序列化的n种方式 1&#xff1a;什么是序列化、反序列化&#xff1f; 把对象转换为字节序列的过程称为对象的序列化把字节序列转换为对象的过程中称为对象的反序列…

我认真写下9段如翔一般的代码,只为等你来品鉴

溪源 | 长沙.NET技术社区开篇我们总是很容易就能写出满足某个特定功能的代码&#xff0c;却很难写出优雅代码。又最欣赏那些优雅的代码&#xff0c;因为优雅代码更能体现一个开发者的积累。就像写一篇散文&#xff0c;有的就像初学者不得其门而入&#xff0c;遣词造句都非常困难…