[NTT][斯特林数]BZOJ 5093: 图的价值
来源:互联网 发布:恋夜秀场破解软件酷安 编辑:程序博客网 时间:2024/05/17 14:19
Description
一个带标号的图的价值定义为每个点度数的
给定
Solution
推式子:
NTT就好了。
#include <bits/stdc++.h>using namespace std;const int G = 3;const int N = 808080;const int MOD = 998244353;typedef long long ll;inline char get(void) { static char buf[100000], *S = buf, *T = buf; if (S == T) { T = (S = buf) + fread(buf, 1, 100000, stdin); if (S == T) return EOF; } return *S++;}template<typename T>inline void read(T &x) { static char c; x = 0; int sgn = 0; for (c = get(); c < '0' || c > '9'; c = get()) if (c == '-') sgn = 1; for (; c >= '0' && c <= '9'; c = get()) x = x * 10 + c - '0'; if (sgn) x = -x;}int g, ig, num, L;int n, p, m, ans;int w[2][N];int R[N];int inv[N], fac[N], ifac[N];int Sp[N], A[N], B[N], C[N];inline int Pow(int a, int b) { int c = 1; while (b) { if (b & 1) c = (ll)c * a % MOD; b >>= 1; a = (ll)a * a % MOD; } return c;}inline int Inv(int x) { return Pow(x, MOD - 2);}void Prep(int n) { g = Pow(G, (MOD - 1) / n); ig = Inv(g); num = n; w[0][0] = w[1][0] = 1; for (int i = 1; i <= n; i++) { w[0][i] = (ll)w[0][i - 1] * ig % MOD; w[1][i] = (ll)w[1][i - 1] * g % MOD; }}inline void FFT(int *a, int n, int r) { static int x, y, INV; for (int i = 0; i < n; i++) if (R[i] > i) swap(a[i], a[R[i]]); for (int i = 1; i < n; i <<= 1) for (int j = 0; j < n; j += (i << 1)) for (int k = 0; k < i; k++) { x = a[j + k]; y = (ll)a[j + k + i] * w[r][num / (i << 1) * k] % MOD; a[j + k] = (x + y) % MOD; a[j + k + i] = (x - y + MOD) % MOD; } if (!r) { INV = Inv(n); for (int i = 0; i < n; i++) a[i] = (ll)a[i] * INV % MOD; }}inline void Add(int &x, int a) { x = (x + a) % MOD;}int main(void) { freopen("1.in", "r", stdin); read(n); read(p); for (m = 1; m <= p; m <<= 1) ++L; m <<= 1; Prep(m); for (int i = 1; i < m; i++) R[i] = (R[i >> 1] >> 1) | ((i & 1) << L); inv[1] = 1; for (int i = 2; i < m; i++) inv[i] = (ll)(MOD - MOD / i) * inv[MOD % i] % MOD; ifac[0] = fac[0] = 1; for (int i = 1; i < m; i++) { ifac[i] = (ll)ifac[i - 1] * inv[i] % MOD; fac[i] = (ll)fac[i - 1] * i % MOD; } for (int i = 0; i <= p; i++) { A[i] = (ll)Pow(i, p) * ifac[i] % MOD; if (i & 1) B[i] = MOD - ifac[i]; else B[i] = ifac[i]; } C[0] = 1; for (int i = 1; i <= p; i++) C[i] = (ll)C[i - 1] * (n - i) % MOD * inv[i] % MOD; FFT(A, m, 1); FFT(B, m, 1); for (int i = 0; i < m; i++) Sp[i] = (ll)A[i] * B[i] % MOD; FFT(Sp, m, 0); for (int i = 0; i <= p; i++) Add(ans, (ll)Sp[i] * fac[i] % MOD * C[i] % MOD * Pow(2, n - i - 1) % MOD); ans = (ll)ans * n % MOD * Pow(2, (ll)(n - 1) * (n - 2) / 2 % (MOD - 1)) % MOD; cout << ans << endl; return 0;}
阅读全文
2 0
- [NTT][斯特林数]BZOJ 5093: 图的价值
- bzoj 5093: 图的价值 NTT+数学+第二类斯特林数
- 【XSY1301】原题的价值 第二类斯特林数 NTT
- 【BZOJ】3992 [SDOI2015]序列统计 【离散对数下的NTT】
- BZOJ 3992 [SDOI2015]序列统计 NTT
- BZOJ 3456 城市规划 NTT 生成函数计数 ***
- bzoj 3992: [SDOI2015]序列统计 NTT
- bzoj 4555: [Tjoi2016&Heoi2016]求和 NTT
- BZOJ 3992 [SDOI2015]序列统计 NTT
- bzoj 4555: [Tjoi2016&Heoi2016]求和 (NTT)
- bzoj 3456: 城市规划 NTT+多项式求逆
- bzoj 2179 FFT快速傅立叶(NTT)
- NTT
- NTT
- 【 bzoj 3992 】 [SDOI2015]序列统计 - NTT 生成函数
- bzoj 3992: [SDOI2015]序列统计 (NTT+快速幂+DP)
- bzoj 3456: 城市规划 (NTT+多项式求逆)
- [生成函数][NTT][多项式求逆]BZOJ 3456: 城市规划
- 数据库_索引总结
- Redis初识
- JAVA核心技术卷2 第四章Socket
- 冒泡排序、快速排序
- LeetCode-028 Implement strStr()
- [NTT][斯特林数]BZOJ 5093: 图的价值
- 哈尔滨理工大学第七届程序设计竞赛决赛(网络赛-高年级组)A 所有情况的和
- debain9的亮度调整失败
- 鼠标事件的使用
- 闭包的理解
- 批量给pdf添加目录(最完整详细方法)
- LeetCode14. Longest Common Prefix
- [DP]Atcoder CODE FESTIVAL 2017 Final (Parallel) D
- Lambda 表达式