luoguP4389 付公主的背包 多项式exp
%%%dkw
话说这是个论文题来着...
考虑生成函数\(OGF\)
对于价值为\(v\)的物品,由于有\(10^5\)的件数,可以看做无限个
那么,其生成函数为\(x^0 + x^{v} + x^{2v} + ... = \frac{1}{1 - x^v}\)
我们所需的答案即\([x^n] \prod \frac{1}{1 - x^{v_i}}\)
只需考虑求出\(A = \prod \frac{1}{1 - x^{v_i}}\)
自然地想到取对数
\(In(A) = \sum In(\frac{1}{1 - x^{v_i}})\)
不难发现
\(In(\frac{1}{1 - x^v}) = - In(1 - x^v)\)
考虑用麦克劳林级数来模拟,那么
由于\(In^{(n)}(1 - x) = - \frac{1}{(1 - x)^n} * (n - 1)!\)
\(-In(1 - x^v) = \sum \frac{x^{vi}}{i}\)
于是,我们可以直接枚举倍数,在\(O(m \log m)\)的时间内完成计算
最后只要\(O(m \log m)\)的\(exp\)一下即可
#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define ri register int
#define rep(io, st, ed) for(ri io = st; io <= ed; io ++)
#define drep(io, ed, st) for(ri io = ed; io >= st; io --)
#define gc getchar
inline int read() {
int p = 0, w = 1; char c = gc();
while(c > '9' || c < '0') { if(c == '-') w = -1; c = gc(); }
while(c >= '0' && c <= '9') p = p * 10 + c - '0', c = gc();
return p * w;
}
const int sid = 500050;
const int mod = 998244353;
int n, m;
int V[sid], F[sid], inv[sid], rev[sid], ans[sid];
inline int Inc(int a, int b) { return (a + b >= mod) ? a + b - mod : a + b; }
inline int Dec(int a, int b) { return (a - b < 0) ? a - b + mod : a - b; }
inline int mul(int a, int b) { return 1ll * a * b % mod; }
inline int fp(int a, int k) {
int ret = 1;
for( ; k; k >>= 1, a = mul(a, a))
if(k & 1) ret = mul(ret, a);
return ret;
}
inline void init(int Maxn, int &n, int &lg) {
n = 1; lg = 0;
while(n < Maxn) n <<= 1, lg ++;
}
inline void NTT(int *a, int n, int opt) {
for(ri i = 0; i < n; i ++) if(i < rev[i]) swap(a[i], a[rev[i]]);
for(ri i = 1; i < n; i <<= 1)
for(ri j = 0, g = fp(3, (mod - 1) / (i << 1)); j < n; j += (i << 1))
for(ri k = j, G = 1; k < i + j; k ++, G = mul(G, g)) {
int x = a[k], y = mul(G, a[i + k]);
a[k] = (x + y >= mod) ? x + y - mod : x + y;
a[i + k] = (x - y < 0) ? x - y + mod : x - y;
}
if(opt == -1) {
int ivn = fp(n, mod - 2);
reverse(a + 1, a + n);
rep(i, 0, n) a[i] = mul(a[i], ivn);
}
}
int ia[sid], ib[sid];
inline void Inv(int *a, int *b, int n) {
if(n == 1) { b[0] = fp(a[0], mod - 2); return; }
Inv(a, b, n >> 1);
int N = 1, lg = 0; init(n + n, N, lg);
for(ri i = 0; i < N; i ++)
rev[i] = (rev[i >> 1] >> 1) | ((i & 1) << (lg - 1));
for(ri i = 0; i < N; i ++) ia[i] = ib[i] = 0;
for(ri i = 0; i < n; i ++) ia[i] = a[i], ib[i] = b[i];
NTT(ia, N, 1); NTT(ib, N, 1);
for(ri i = 0; i < N; i ++)
ia[i] = Dec((ib[i] << 1) % mod, mul(ia[i], mul(ib[i], ib[i])));
NTT(ia, N, -1);
for(ri i = 0; i < n; i ++) b[i] = ia[i];
}
inline void Inv_init(int n) {
inv[0] = inv[1] = 1;
rep(i, 2, n) inv[i] = mul(inv[mod % i], mod - mod / i);
}
inline void wf(int *a, int *b, int n) { for(ri i = 1; i < n; i ++) b[i - 1] = mul(a[i], i); }
inline void jf(int *a, int *b, int n) { for(ri i = 1; i < n; i ++) b[i] = mul(a[i - 1], inv[i]); }
int iv[sid], dx[sid];
inline void In(int *a, int *b, int n) {
for(ri i = 0; i < n + n; i ++) iv[i] = dx[i] = 0;
Inv(a, iv, n); wf(a, dx, n);
int N = 1, lg = 0; init(n + n, N, lg);
for(ri i = 0; i < N; i ++)
rev[i] = (rev[i >> 1] >> 1) | ((i & 1) << (lg - 1));
NTT(iv, N, 1); NTT(dx, N, 1);
for(ri i = 0; i < N; i ++) iv[i] = mul(iv[i], dx[i]);
NTT(iv, N, -1); jf(iv, b, n);
}
int inb[sid], fb[sid];
inline void Exp(int *a, int *b, int n) {
if(n == 1) { b[0] = 1; return; }
Exp(a, b, n >> 1);
for(ri i = 0; i < n + n; i ++) inb[i] = fb[i] = 0;
In(b, inb, n);
int N = 1, lg = 0; init(n + n, N, lg);
for(ri i = 0; i < N; i ++)
rev[i] = (rev[i >> 1] >> 1) | ((i & 1) << (lg - 1));
for(ri i = 0; i < n; i ++) fb[i] = Dec(a[i], inb[i]); fb[0] ++;
for(ri i = 0; i < n; i ++) inb[i] = b[i];
NTT(inb, N, 1); NTT(fb, N, 1);
for(ri i = 0; i < N; i ++) fb[i] = mul(fb[i], inb[i]);
NTT(fb, N, -1);
for(ri i = 0; i < n; i ++) b[i] = fb[i], b[i + n] = 0;
}
inline void calc() {
int N = 1, lg = 0;
init(m + 5, N, lg); Inv_init(N);
for(ri i = 1; i <= m; i ++)
for(ri j = i; j <= m; j += i)
F[j] = Inc(F[j], mul(V[i], inv[j / i]));
Exp(F, ans, N);
rep(i, 1, m) printf("%d\n", ans[i]);
}
int main() {
n = read(); m = read();
rep(i, 1, n) V[read()] ++;
calc();
return 0;
}
luoguP4389 付公主的背包 多项式exp的更多相关文章
- [luogu4389]付公主的背包(多项式exp)
完全背包方案计数问题的FFT优化.首先写成生成函数的形式:对重量为V的背包,它的生成函数为$\sum\limits_{i=0}^{+\infty}x^{Vi}=\frac{1}{1-x^{V}}$于是 ...
- LuoguP4389 付公主的背包【生成函数+多项式exp】
题目背景 付公主有一个可爱的背包qwq 题目描述 这个背包最多可以装10^5105大小的东西 付公主有n种商品,她要准备出摊了 每种商品体积为Vi,都有10^5105件 给定m,对于s\in [1,m ...
- 洛谷 4389 付公主的背包——多项式求ln、exp
题目:https://www.luogu.org/problemnew/show/P4389 关于泰勒展开: https://blog.csdn.net/SoHardToNamed/article/d ...
- luoguP4389 付公主的背包
luogu 显然这是个背包题 显然物品的数量是不用管的 所以考虑大小为\(v\)的物品可以装的体积用生成函数表示一下 \[ f(x)=\sum_{i=0}^{+\infty}x^{vi}=\frac{ ...
- [题解] LuoguP4389 付公主的背包
这个题太神辣- 暴力背包就能获得\(30\)分的好成绩...... \(60\)分不知道咋搞..... 所以直接看\(100\)分吧\(QwQ\) 用一点生成函数的套路,对于一个体积为\(v\)的物品 ...
- 洛谷 P4389 付公主的背包 解题报告
P4389 付公主的背包 题目背景 付公主有一个可爱的背包qwq 题目描述 这个背包最多可以装\(10^5\)大小的东西 付公主有\(n\)种商品,她要准备出摊了 每种商品体积为\(V_i\),都有\ ...
- Luogu4389 付公主的背包(生成函数+多项式exp)
显然构造出生成函数,对体积v的物品,生成函数为1+xv+x2v+……=1/(1-xv).将所有生成函数乘起来得到的多项式即为答案,设为F(x),即F(x)=1/∏(1-xvi).但这个多项式的项数是Σ ...
- 【Luogu4389】付公主的背包
题目 传送门 解法 答案显然是\(n\)个形如\(\sum_{i \geq 1} x^{vi}\)的多项式的卷积 然而直接NTT的时间复杂度是\(O(nm\log n)\) 我们可以把每个多项式求\( ...
- 洛谷P4389 付公主的背包--生成函数+多项式
题目链接戳这里 题目描述 有\(n\)件不同的商品,每件物品都有无限个,输出总体积为\([1,m]\)的方案数 思路 直接跑背包有\(30\) 考虑把每个物品的生成函数设出来,对于一件体积为\(v\) ...
随机推荐
- 【DS】排序算法之选择排序(Selection Sort)
一.算法思想 选择排序是一种简单直观的排序算法.它的工作原理如下: 1)将序列分成两部分,前半部分是已经排序的序列,后半部分是未排序的序列: 2)在未排序序列中找到最小(大)元素,放到已排序序列的末尾 ...
- [转载]Understanding the Bootstrap 3 Grid System
https://scotch.io/tutorials/understanding-the-bootstrap-3-grid-system With the 3rd version of the gr ...
- 基于canvas将图片转化成字符画
字符画大家一定非常熟悉了,那么如何把一张现有的图片转成字符画呢?HTML5让这个可能变成了现实,通过canvas,可以很轻松实现这个功能.其实原理很简单:扫描图片相应位置的像素点,再计算出其灰度值,根 ...
- javascript完美拖拽的实现
直接上代码: HTML代码: <!DOCTYPE HTML> <html lang="en-US"> <head> <meta chars ...
- 洛谷 P4838 P哥破解密码 题解
矩阵乘法 + 快速幂优化递推: 看到这个题目我们不难想到递推,题干中说3个连续的A出现在序列中是不合法的,所以可以分为三种情况: (1):序列前只有一个A,如:BA,BBA,BABA. (2):序列前 ...
- 转载 python多重继承C3算法
备注:O==object 2.python-C3算法解析: #C3 定义引用开始 C3 算法:MRO是一个有序列表L,在类被创建时就计算出来. L(Child(Base1,Base2)) = [ Ch ...
- 洛谷 P4389: 付公主的背包
题目传送门:洛谷 P4389. 题意简述: 有 \(n\) 个物品,每个物品都有无限多,第 \(i\) 个物品的体积为 \(v_i\)(\(v_i\le m\)). 问用这些物品恰好装满容量为 \(i ...
- blog迁移
blog迁移到: https://github.com/for-firecat/
- linux C守护进程编写
linux编程-守护进程编写 守护进程(Daemon)是运行在后台的一种特殊进程.它独立于控制终端并且周期性地执行某种任务或等待 处理某些发生的事件.守护进程是一种很有用的进程. Linux的大多数服 ...
- SOAP简单示例
看了网上的几个文章,SOAP的示例布局都不清晰,不能马上入手,特意写个例子与大家分享,同时记录备用. 当前环境:VS2013 + WPF private void Button_Click(objec ...