HDU 5552 CDQ分治+NTT
来源:互联网 发布:数码绘画软件 编辑:程序博客网 时间:2024/06/07 10:51
题意:
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5552
统计n个点可以组成多少种不同的有环连通图,其中每条边可以染成m种颜色。
思路:
分别统计f(n)为n个点的连通图数目,g(n)为n个点能构成的图数目,h(n)为n个点能构成的树的数目,然后递推出公式,发现是个卷积形式,可以用FFT加速,这里会爆精度,改用NTT。
代码:
#include <bits/stdc++.h>using namespace std;typedef long long LL;const LL MOD = 152076289;const int MAXN = 4e4 + 10;const int G = 106;LL x1[MAXN], x2[MAXN];LL pow_mod(LL a, LL n) { LL res = 1; while (n) { if (n & 1) res = res * a % MOD; a = a * a % MOD; n >>= 1; } return res;}void change (LL *y, int len) { int i, j, k; for(i = 1, j = len / 2; i < len - 1; i++) { if (i < j) swap(y[i], y[j]); k = len / 2; while (j >= k) { j -= k; k /= 2; } if (j < k) j += k; }}void ntt (LL *y, int len, int on) { change (y, len); int id = 0; for(int h = 2; h <= len; h <<= 1) { id++; LL wn = pow_mod (G, (MOD - 1) / (1<<id)); for(int j = 0; j < len; j += h) { LL w = 1; for(int k = j; k < j + h / 2; k++) { LL u = y[k] % MOD; LL t = w * (y[k + h / 2] % MOD) % MOD; y[k] = (u + t) % MOD; y[k + h / 2] = ((u - t) % MOD + MOD) % MOD; w = w * wn % MOD; } } } if (on == -1) { for (int i = 1; i < len / 2; i++) swap (y[i], y[len - i]); LL inv = pow_mod(len, MOD - 2); for(int i = 0; i < len; i++) y[i] = y[i] % MOD * inv % MOD; }}LL f[MAXN], g[MAXN], h[MAXN];LL inv[MAXN], fac[MAXN];void init() { fac[0] = 1; inv[0] = 1; for (int i = 1; i <= 40000; i++) { fac[i] = fac[i - 1] * i % MOD; inv[i] = pow_mod(fac[i], MOD - 2); }}void solve(int l, int r) { if (l == r) { f[l] = (f[l] + g[l]) % MOD; return; } int m = (l + r) >> 1; solve(l, m); int len = 1; while (len <= r - l + 1) len <<= 1; for (int i = 0; i < len; i++) x1[i] = x2[i] = 0; for (int i = l; i <= m; i++) x1[i - l] = f[i] * inv[i - 1] % MOD; for (int i = 1; i <= r - l; i++) x2[i - 1] = g[i] * inv[i] % MOD; ntt(x1, len, 1); ntt(x2, len, 1); for (int i = 0; i < len; i++) x1[i] = x1[i] * x2[i] % MOD; ntt(x1, len, -1); for (int i = m + 1; i <= r; i++) { f[i] -= x1[i - l - 1] % MOD * fac[i - 1] % MOD; f[i] = (f[i] + MOD) % MOD; } solve(m + 1, r);}int main() { //freopen("in.txt", "r", stdin); init(); int T, cs = 0; scanf("%d", &T); while (T--) { int n, m; scanf("%d%d", &n, &m); for (int i = 1; i <= n; i++) g[i] = pow_mod(m + 1, i * (i - 1) / 2); memset(f, 0, sizeof(f)); solve(1, n); LL ans = f[n]; ans -= pow_mod(n, n - 2) * pow_mod(m, n - 1) % MOD; ans = (ans + MOD) % MOD; printf("Case #%d: %I64d\n", ++cs, ans); } return 0;}
阅读全文
0 0
- HDU 5552 (CDQ分治 NTT)
- HDU 5552 CDQ分治+NTT
- HDU 5896 (CDQ分治 NTT)
- HDU 5896 CDQ分治+NTT
- HDU 5322 Hope (CDQ分治+NTT)
- 【HDU】5279 YJC plays Minecraft【cdq分治+NTT】
- BNUOJ51279 组队活动(cdq分治&&NTT)
- 【cdq分治&NTT】BNUOJ51279组队活动
- HDU 4456 CDQ分治
- HDU-5324 cdq分治
- HDU - 1166 CDQ分治
- Hdu 5618 CDQ分治
- HDU 5322 NTT与分治
- hdu 5322 Hope 分治 NTT
- zoj 3874 Permutation Graph (cdq分治+NTT)
- BZOJ4836(CDQ分治+FFT(NTT会T))
- hdu 5126 stars cdq分治
- HDU 4456 Crowd (cdq分治)
- Java多线程编程学习chapter1
- caffe学习:数据层参数
- Traincascade Error: Bad argument (Can not get new positive sample. Themost possible reason is insuff
- (Y)小米miui系统调用系统裁剪图片功能报错 java.lang.RuntimeException: Failure delivering result ResultInfo{who=null, r
- 对于Java第一次面试的建议
- HDU 5552 CDQ分治+NTT
- Python语言常见异常错误类型
- 图像时域滤波(二)
- 论文都看不懂,你还搞什么人工智能?
- js获取IP地址 mac地址 主机名称
- 数据结构——第二讲、线性结构(2)
- LESS CSS 框架简介
- 趣味编程算法类文章征集列表
- EditText 取消多余的内边距