这个题太神辣~


暴力背包就能获得\(30\)分的好成绩......

\(60\)分不知道咋搞.....

所以直接看\(100\)分吧\(QwQ\)

用一点生成函数的套路,对于一个体积为\(v\)的物品,我们构造一个序列\(f_n = [v \mid n]\ (n \ge 0)\)

其生成函数\(F(x) = \sum\limits_{i=0}^{\infty} [v \mid i]x^i = \sum\limits_{i=0}^{\infty} x^{vi} = \frac{1}{1-x^v}\)

那么答案就是这\(n\)个\(F\)的乘积。

考虑怎么把这个东西乘起来,大力乘是\(O(nm \log m)\),其封闭形式\(\frac{1}{1-x^v}\)好像也没有什么好优化的......

一个套路的想法是把每个\(F\),先\(\ln\)一下,然后再加起来,最后\(\exp\)回去。

问题变成了怎么求\(\ln F(x)\),不放设\(\ln F(x) = G(x)\)。

套路的两边对\(x\)求导

\[\frac{F'(x)}{F(x)} = G'(x)
\]

注意到\(F(x) = \frac{1}{1-x^v}\)

\[\begin{aligned}
(1 - x^v)\sum\limits_{i=0}^{\infty} vix^{vi-1} &= G'(x) \\
\sum\limits_{i=0}^{\infty} vix^{vi-1} - \sum\limits_{i=0}^{\infty} vix^{v(i+1)-1} &= G'(x)
\end{aligned}
\]

注意到\(i=0\)时候两个\(\Sigma\)都是\(0\),所以

\[\begin{aligned}
\sum\limits_{i=1}^{\infty} vix^{vi-1} - \sum\limits_{i=1}^{\infty} v(i-1) x^{vi-1} &= G'(x) \\
\sum\limits_{i=1}^{\infty} vx^{vi-1} &= G'(x)
\end{aligned}
\]

两边同时积分回去

\[G(x) = \int\sum\limits_{i=1}^{\infty} vx^{vi-1}dx = \sum\limits_{i=1}^{\infty}v \cdot \frac{1}{vi}x^{vi} = \sum\limits_{i=1}^{\infty} \frac{1}{i}x^{vi}
\]

可以对每个\(v\),\(O(\frac{m}{v})\)的搞出\(G\),然后\(\exp\)回去就好了。

\(v\)有重复的话开个桶可以保证复杂度,预处理逆元能快一点。

\(Code:\)

#include <bits/stdc++.h>
using namespace std;
const int N=3e5+10,P=998244353,gen=3,igen=(P+1)/gen;
int add(int x,int y){return x+y>=P?x+y-P:x+y;}
int sub(int x,int y){return x-y<0?x-y+P:x-y;}
int fpow(int x,int y){
int ret=1; for(x%=P;y;y>>=1,x=1ll*x*x%P)
if(y&1) ret=1ll*x*ret%P;
return ret;
}
namespace Poly{
int rev[N];
void init(int n){
for(int i=0;i<n;i++)
rev[i]=rev[i>>1]>>1|((i&1)?n>>1:0);
}
void ntt(int *f,int n,int flg){
for(int i=0;i<n;i++)
if(rev[i]<i) swap(f[i],f[rev[i]]);
for(int len=2,k=1;len<=n;len<<=1,k<<=1){
int wn=fpow(flg==1?gen:igen,(P-1)/len);
for(int i=0;i<n;i+=len)
for(int j=i,w=1;j<i+k;j++,w=1ll*w*wn%P){
int tmp=1ll*w*f[j+k]%P;
f[j+k]=sub(f[j],tmp),f[j]=add(f[j],tmp);
}
}
if(flg==-1){
int inv=fpow(n,P-2);
for(int i=0;i<n;i++) f[i]=1ll*f[i]*inv%P;
}
}
#define glim(n) int limit=1;while(limit<=(n))limit<<=1;init(limit)
void getinv(int *f,int n,int *G){
if(n==1){G[0]=fpow(f[0],P-2);return;}
getinv(f,(n+1)>>1,G); glim(2*n);
static int F[N];
for(int i=0;i<limit;i++)
i>=n?F[i]=G[i]=0:F[i]=f[i],G[i]=G[i];
ntt(F,limit,1),ntt(G,limit,1);
for(int i=0;i<limit;i++) G[i]=1ll*G[i]*sub(2,1ll*F[i]*G[i]%P)%P;
ntt(G,limit,-1);
for(int i=n;i<limit;i++) G[i]=0;
}
void dao(int *f,int n,int *G){
static int F[N]; for(int i=0;i<=n;i++) F[i]=f[i];
for(int i=1;i<=n;i++) G[i-1]=1ll*F[i]*i%P; G[n]=0;
}
void jifen(int *f,int n,int *G){
static int F[N]; for(int i=0;i<=n;i++) F[i]=f[i];
for(int i=0;i<=n;i++) G[i+1]=1ll*F[i]*fpow(i+1,P-2)%P; G[0]=0;
}
void getln(int *f,int n,int *G){
static int F[N],iF[N]; for(int i=0;i<n;i++) F[i]=f[i];
getinv(F,n,iF),dao(F,n-1,F); glim(n*2);
for(int i=0;i<limit;i++) F[i]=i>=n?0:F[i],iF[i]=i>=n?0:iF[i];
ntt(F,limit,1),ntt(iF,limit,1);
for(int i=0;i<limit;i++) G[i]=1ll*F[i]*iF[i]%P;
ntt(G,limit,-1),jifen(G,n-1,G);
for(int i=n;i<limit;i++) G[i]=0;
}
void getexp(int *f,int n,int *G){
if(n==1){G[0]=1;return;}
getexp(f,(n+1)>>1,G); glim(n*2);
static int F[N],lnG[N];
for(int i=0;i<limit;i++) i>=n?F[i]=G[i]=0:F[i]=f[i],G[i]=G[i];
getln(G,n,lnG),ntt(F,limit,1),ntt(G,limit,1),ntt(lnG,limit,1);
for(int i=0;i<limit;i++) G[i]=1ll*G[i]*add(sub(1,lnG[i]),F[i])%P;
ntt(G,limit,-1);
for(int i=n;i<limit;i++) G[i]=0;
}
}
int F[N],ans[N],inv[N],cnt[N];
int main(){
int n,m; scanf("%d%d",&n,&m);
inv[1]=1; for(int i=2;i<=m;i++)
inv[i]=1ll*(P-P/i)*inv[P%i]%P;
for(int _=1;_<=n;_++){
int v;scanf("%d",&v);
if(v<=m) cnt[v]=add(cnt[v],1);
}
for(int v=1;v<=m;v++) if(cnt[v])
for(int i=1;i*v<=m;i++) F[i*v]=add(F[i*v],1ll*cnt[v]*inv[i]%P);
Poly::getexp(F,m+1,ans);
for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
return 0;
}

[题解] LuoguP4389 付公主的背包的更多相关文章

  1. LuoguP4389 付公主的背包【生成函数+多项式exp】

    题目背景 付公主有一个可爱的背包qwq 题目描述 这个背包最多可以装10^5105大小的东西 付公主有n种商品,她要准备出摊了 每种商品体积为Vi,都有10^5105件 给定m,对于s\in [1,m ...

  2. luoguP4389 付公主的背包

    luogu 显然这是个背包题 显然物品的数量是不用管的 所以考虑大小为\(v\)的物品可以装的体积用生成函数表示一下 \[ f(x)=\sum_{i=0}^{+\infty}x^{vi}=\frac{ ...

  3. luoguP4389 付公主的背包 多项式exp

    %%%dkw 话说这是个论文题来着... 考虑生成函数\(OGF\) 对于价值为\(v\)的物品,由于有\(10^5\)的件数,可以看做无限个 那么,其生成函数为\(x^0 + x^{v} + x^{ ...

  4. 洛谷 P4389 付公主的背包 解题报告

    P4389 付公主的背包 题目背景 付公主有一个可爱的背包qwq 题目描述 这个背包最多可以装\(10^5\)大小的东西 付公主有\(n\)种商品,她要准备出摊了 每种商品体积为\(V_i\),都有\ ...

  5. 洛谷 P4389: 付公主的背包

    题目传送门:洛谷 P4389. 题意简述: 有 \(n\) 个物品,每个物品都有无限多,第 \(i\) 个物品的体积为 \(v_i\)(\(v_i\le m\)). 问用这些物品恰好装满容量为 \(i ...

  6. 洛谷P4389 付公主的背包--生成函数+多项式

    题目链接戳这里 题目描述 有\(n\)件不同的商品,每件物品都有无限个,输出总体积为\([1,m]\)的方案数 思路 直接跑背包有\(30\) 考虑把每个物品的生成函数设出来,对于一件体积为\(v\) ...

  7. luogu P4389 付公主的背包

    传送门 神仙题鸭!orz dkw 暴力就是完全背包 而完全背包可以和生成函数扯上关系,记第i种物品质量为\(a_i\),那么这种物品的生成函数\(G(i)=\sum_{j=0}^{\infty}x^{ ...

  8. luogu4389 付公主的背包

    题目链接:洛谷 题目大意:现在有$n$个物品,每种物品体积为$v_i$,对任意$s\in [1,m]$,求背包恰好装$s$体积的方案数(完全背包问题). 数据范围:$n,m\leq 10^5$ 这道题 ...

  9. P3489 付公主的背包

    题意:n<=1e5,m<=1e5,跑n个物品1到m容量的完全背包. 考虑暴力的做法就是把一些1/(1+x^a)的多项式乘起来即可. 考虑优化,取一下ln,转化为加法,然后exp回去就好了.

随机推荐

  1. Windows7 wampServer3.0.6 Mutillidae2.7.12

    在Mac上访问虚拟机中的mutillidae,报403: By default, Mutillidae only allow access from localhost ***: Parallels ...

  2. Navicat导入json文件到数据库

    最近做小程序商城系统,物流管理这块需要存储国际和中国的物流地址. 所以,天哪,地址那么的多!!!! www.baidu.com,搜索大佬是如何做的.有的是通过访问阿里云快递物流api接口获取数据存入数 ...

  3. vue-i18n突然所有的都解析不出来了

    因为浏览器cookie缓存的语言标志与代码中不一致

  4. logback日志

    一.什么是日志框架? 是一套能够实现日志输出的工具包 能够描述系统运行状态的所有时间都可以算作日志 用户下线,接口超时,数据崩溃 二.日志框架的能力 1.定制输出目标(文件,回滚策略,数据库,网络的第 ...

  5. CPU、内存、硬盘之间的关系

    要完完全全地讲清楚cpu.内存.硬盘之间的关系,博客的篇幅是不够的.这里简单的介绍以下它们之间的关系,抛砖引玉. 1.CPU即中央处理器,是英语“Central Processing Unit”的缩写 ...

  6. webpack 命令 Module build failed (from ./node_modules/babel-loader/lib/index.js) 错误问题解决方案

    在项目中运行的时候出现报错,错误为Module build failed (from ./node_modules/babel-loader/lib/index.js) 解决方案: 控制台输入  np ...

  7. element设置headers添加token

    <template>   <div>     <el-upload       action="http://localhost:3000/picture&qu ...

  8. python里的def 方法中->代表什么意思?

    功能注释 函数注释是关于用户定义函数使用的类型的完全可选元数据信息(请参阅PEP 3107和 PEP 484了解更多信息). 注释__annotations__ 作为字典存储在函数的属性中,对函数的任 ...

  9. JS写一个列表跑马灯效果--基于touchslide.js

    先放上效果图: 类似于这样的,在列表中循环添加背景样式的跑马灯效果. 准备引入JS插件: <script type="text/javascript" src="x ...

  10. 数学软件实训1-MATLAB程序设计及应用初步

    数学软件实训任务一 一 题目:MATLAB程序设计及应用初步 二 目的:掌握MATLAB程序设计的基本方法,会利用MATLAB程序设计思想编程处理一些简单问题. 三 要求: 1 掌握控制流的基本语法结 ...