【刷题】LOJ 556 「Antileaf's Round」咱们去烧菜吧

最新推荐文章于 2024-07-19 15:23:49 发布
最新推荐文章于 2024-07-19 15:23:49 发布
阅读量100 收藏
点赞数
文章标签： ui
原文链接：http://www.cnblogs.com/hongyj/p/10518744.html
版权

题目描述

烧菜

你有 \(m\) 种物品，第 \(i\) 种物品的大小为 \(a_i\) ​，数量为 \(b_i\)​（ \(b_i=0\) 表示有无限个）。

你还有 \(n\) 个背包，体积分别为 \(1\) 到 \(n\) ，现在你很想知道用这些物品填满某个背包的方案数。

为了满足你的好奇心，你决定把填满每个背包的方案数都算一遍。

因为你其实只是闲得无聊，所以你只想知道方案数对 \(998244353\)（ \(7\times 17\times 2^{23}+1\)，一个质数）取模后的值。

输入格式

第一行两个非负整数，分别表示 \(n,m\)。

以下 \(m\) 行，每行两个非负整数，分别表示 \(a_i,b_i\) ​。

输出格式

输出 \(n\) 个非负整数表示答案。

样例

样例输入

5 3
1 0
1 1
3 2

样例输出

2
2
3
4
4

样例解释

拼出 \(1\) ~ \(5\) 的方案分别如下：

\(\{1_1\},\{1_2\}\)

\(\{1_1,1_1\},\{1_1,1_2\}\)

\(\{1_1,1_1,1_1\},\{1_1,1_1,1_2\},\{3\}\)

\(\{1_1,1_1,1_1,1_1\},\{1_1,1_1,1_1,1_2\},\{1_1,3\},\{1_2,3\}\)

\(\{1_1,1_1,1_1,1_1,1_1\},\{1_1,1_1,1_1,1_1,1_2\},\{1_1,1_1,3\},\{1_1,1_2,3\}\)

数据范围与提示

\(0<n,m \leq 10^5,0 \leq a_i​ \leq 110000,0 \leq b_i \leq 10^6\)。

题解

数量无限的物品的生成函数为：\(\sum_{i=0}^{+\infty}x^{ia_i}=\frac{1}{1-x^{a_i}}\)

数量有限的物品的生成函数为：\(\sum_{i=0}^{b_i}x^{ia_i}=\frac{1-x^{(b_i+1)a_i}}{1-x^{a_i}}\)

将所有物品的生成函数乘起来就是答案的生成函数，即（假设全部都有限）：

\(\begin{aligned} ans &= \prod_{i=1}^m\frac{1-x^{(b_i+1)a_i}}{1-x^{a_i}} \\ &= \exp{\sum_{i=1}^m (\ln (1-x^{(b_i+1)a_i})-\ln (1-x^{a_i}))} \end{aligned}\)

然后我们知道：\(-\ln(1-x)=\sum_{i=1}^{+\infty}\frac{x^i}{i}\) ，所以

\(\ln (1-x^{(b_i+1)a_i})=\sum_{i=1}^{+\infty}-\frac{1}{i}x^{(b_i+1)a_i}\)

\(\ln (1-x^{a_i})=\sum_{i=1}^{+\infty}-\frac{1}{i}x^{a_i}\)

那么对于物品，可以 \(O(n \ln{n} )\) 地把它们的贡献加到OGF上，然后做多项式exp就好了

#include<bits/stdc++.h>
#define ui unsigned int
#define ll long long
#define db double
#define ld long double
#define ull unsigned long long
#define ft first
#define sd second
#define pb(a) push_back(a)
#define PII std::pair<int,int>
#define PLL std::pair<ll,ll>
#define mp(a,b) std::make_pair(a,b)
#define ITR(a,b) for(auto a:b)
#define REP(a,b,c) for(register int a=(b),a##end=(c);a<=a##end;++a)
#define DEP(a,b,c) for(register int a=(b),a##end=(c);a>=a##end;--a)
const int MAXN=1<<21,Mod=998244353;
int M[MAXN];
ll OGF[MAXN],ans[MAXN];
template<typename T> inline void read(T &x)
{
    T data=0,w=1;
    char ch=0;
    while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
    if(ch=='-')w=-1,ch=getchar();
    while(ch>='0'&&ch<='9')data=((T)data<<3)+((T)data<<1)+(ch^'0'),ch=getchar();
    x=data*w;
}
template<typename T> inline void write(T x,char ch='\0')
{
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+'0');
    if(ch!='\0')putchar(ch);
}
template<typename T> inline bool chkmin(T &x,T y){return y<x?(x=y,true):false;}
template<typename T> inline bool chkmax(T &x,T y){return y>x?(x=y,true):false;}
template<typename T> inline T min(T x,T y){return x<y?x:y;}
template<typename T> inline T max(T x,T y){return x>y?x:y;}
namespace poly{
    const int Mod=998244353;
    int rev[MAXN];
    ll b[MAXN],c[MAXN],tmp[MAXN],inv2,inv3,invm,pw[MAXN],ipw[MAXN],iv[MAXN];
    inline ll qexp(ll a,ll b)
    {
        ll res=1;
        while(b)
        {
            if(b&1)res=res*a%Mod;
            a=a*a%Mod;
            b>>=1;
        }
        return res;
    }
    inline void init(int n)
    {
        inv2=qexp(2,Mod-2);inv3=qexp(3,Mod-2);
        for(register int l=1;l<=(n<<2);l<<=1)
            pw[l]=qexp(3,(Mod-1)/l),ipw[l]=qexp(inv3,(Mod-1)/l);
        REP(i,1,n+n+1)iv[i]=qexp(i,Mod-2);
    }
    inline void NTT(ll *A,int n,int tp)
    {
        REP(i,0,n-1)if(i<rev[i])std::swap(A[i],A[rev[i]]);
        for(register int l=2;l<=n;l<<=1)
        {
            ll wn=tp>0?pw[l]:ipw[l];
            for(register int i=0;i<n;i+=l)
            {
                ll w=1;
                for(register int j=0;j<(l>>1);++j)
                {
                    ll A1=A[i+j],A2=1ll*A[i+j+(l>>1)]*w%Mod;
                    A[i+j]=(A1+A2)%Mod,A[i+j+(l>>1)]=(A1-A2+Mod)%Mod;
                    w=1ll*w*wn%Mod;
                }
            }
        }
    }
    inline void Get_Mul(int n,ll *A,ll *B,ll *C)
    {
        int m,cnt;
        for(m=1,cnt=0;m<n;m<<=1,++cnt);
        REP(i,0,m-1)rev[i]=(rev[i>>1]>>1)|((i&1)<<(cnt-1));
        NTT(A,m,1);NTT(B,m,1);
        REP(i,0,m-1)A[i]=1ll*A[i]*B[i]%Mod;
        NTT(A,m,-1);
        invm=qexp(m,Mod-2);
        REP(i,0,n-1)C[i]=1ll*A[i]*invm%Mod;
    }
    inline void Get_Inv(int len,ll *A,ll *B)
    {
        if(len==1)
        {
            B[0]=qexp(A[0],Mod-2);
            return ;
        }
        Get_Inv((len+1)>>1,A,B);
        int m,cnt;
        for(m=1,cnt=0;m<(len<<1);m<<=1,++cnt);
        REP(i,0,m-1)rev[i]=(rev[i>>1]>>1)|((i&1)<<(cnt-1));
        REP(i,0,len-1)tmp[i]=A[i];
        REP(i,len,m-1)tmp[i]=0;
        NTT(tmp,m,1);NTT(B,m,1);
        REP(i,0,m-1)B[i]=(2ll*B[i]%Mod-1ll*B[i]*B[i]%Mod*tmp[i]%Mod+Mod)%Mod;
        NTT(B,m,-1);
        invm=qexp(m,Mod-2);
        REP(i,0,len-1)B[i]=1ll*B[i]*invm%Mod;
        REP(i,len,m-1)B[i]=0;
    }
    inline void Get_Sqr(int len,ll *A,ll *B)
    {
        if(len==1)
        {
            B[0]=1;
            return ;
        }
        Get_Sqr((len+1)>>1,A,B);
        memset(b,0,sizeof(b));Get_Inv(len,B,b);
        int m,cnt;
        for(m=1,cnt=0;m<(len<<1);m<<=1,++cnt);
        REP(i,0,m-1)rev[i]=(rev[i>>1]>>1)|((i&1)<<(cnt-1));
        REP(i,0,len-1)tmp[i]=A[i];
        REP(i,len,m-1)tmp[i]=0;
        NTT(B,m,1);NTT(tmp,m,1);NTT(b,m,1);
        REP(i,0,m-1)B[i]=1ll*(tmp[i]+1ll*B[i]*B[i]%Mod)%Mod*b[i]%Mod*inv2%Mod;
        NTT(B,m,-1);
        invm=qexp(m,Mod-2);
        REP(i,0,len-1)B[i]=1ll*B[i]*invm%Mod;
        REP(i,len,m-1)B[i]=0;
    }
    inline void Get_Div(int n1,ll *A,int n2,ll *B,ll *C,ll *D)
    {
        memset(c,0,sizeof(c));
        REP(i,0,n2-1)c[i]=B[n2-1-i];
        Get_Inv(n1-n2+1,c,b);
        int m,cnt;
        for(m=1,cnt=0;m<((n1-n2+1)<<1);m<<=1,++cnt);
        REP(i,0,m-1)rev[i]=(rev[i>>1]>>1)|((i&1)<<(cnt-1));
        REP(i,0,n1-n2)tmp[i]=A[n1-1-i];
        REP(i,n1-n2+1,m-1)tmp[i]=0;
        NTT(tmp,m,1);NTT(b,m,1);
        REP(i,0,m-1)tmp[i]=1ll*tmp[i]*b[i]%Mod;
        NTT(tmp,m,-1);
        invm=qexp(m,Mod-2);
        memset(c,0,sizeof(c));
        REP(i,0,n1-n2)C[n1-n2-i]=c[n1-n2-i]=1ll*tmp[i]*invm%Mod;
        for(m=1,cnt=0;m<n1;m<<=1,++cnt);
        REP(i,0,m-1)rev[i]=(rev[i>>1]>>1)|((i&1)<<(cnt-1));
        NTT(B,m,1);NTT(c,m,1);
        REP(i,0,m-1)B[i]=1ll*B[i]*c[i]%Mod;
        NTT(B,m,-1);
        invm=qexp(m,Mod-2);
        REP(i,0,n2-2)D[i]=(A[i]-1ll*B[i]*invm%Mod+Mod)%Mod;
    }
    inline void Get_Der(int n,ll *A,ll *B)
    {
        REP(i,0,n-2)B[i]=1ll*A[i+1]*(i+1)%Mod;B[n-1]=0;
    }
    inline void Get_Int(int n,ll *A,ll *B)
    {
        REP(i,1,n-1)B[i]=1ll*A[i-1]*iv[i]%Mod;B[0]=0;
    }
    inline void Get_Ln(int n,ll *A,ll *B)
    {
        memset(b,0,sizeof(b));
        Get_Inv(n,A,b);Get_Der(n,A,c);
        int m,cnt;
        for(m=1,cnt=0;m<(n<<1);m<<=1,++cnt);
        REP(i,0,m-1)rev[i]=(rev[i>>1]>>1)|((i&1)<<(cnt-1));
        NTT(b,m,1);NTT(c,m,1);
        REP(i,0,m-1)b[i]=1ll*b[i]*c[i]%Mod;
        NTT(b,m,-1);
        invm=qexp(m,Mod-2);
        REP(i,0,n-1)b[i]=1ll*b[i]*invm%Mod;
        Get_Int(n,b,B);
    }
    inline void Get_Exp(int len,ll *A,ll *B)
    {
        if(len==1)
        {
            B[0]=1;
            return ;
        }
        Get_Exp((len+1)>>1,A,B);Get_Ln(len,B,c);
        int m,cnt;
        for(m=1,cnt=0;m<(len<<1);m<<=1,++cnt);
        REP(i,0,m-1)rev[i]=(rev[i>>1]>>1)|((i&1)<<(cnt-1));
        REP(i,0,len-1)tmp[i]=(A[i]-c[i]+(i==0)+Mod)%Mod;
        REP(i,len,m-1)tmp[i]=0;
        NTT(tmp,m,1);NTT(B,m,1);
        REP(i,0,m-1)B[i]=1ll*tmp[i]*B[i]%Mod;
        NTT(B,m,-1);
        invm=qexp(m,Mod-2);
        REP(i,0,len-1)B[i]=1ll*B[i]*invm%Mod;
        REP(i,len,m-1)B[i]=0;
    }
}
int main()
{
    int n,m;read(n);read(m);poly::init(n+1);
    REP(i,1,m)
    {
        int a,b;read(a);read(b);
        if(a<=n)M[a]++;
        if(1ll*(b+1)*a<=n&&b)M[(b+1)*a]--;
    }
    REP(i,1,n)if(M[i])REP(j,1,n/i)OGF[i*j]=((OGF[i*j]+1ll*M[i]*poly::iv[j]%Mod)%Mod+Mod)%Mod;
    poly::Get_Exp(n+1,OGF,ans);
    REP(i,1,n)write(ans[i],'\n');
    return 0;
}

转载于:https://www.cnblogs.com/hongyj/p/10518744.html

weixin_30376083
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
【LibreOJ #556.「Antileaf's Round」咱们去烧菜吧】【生成函数+牛顿迭代】
beginend
08-13 228
题意 有mmm种物品，第iii种物品的大小是aia_iai​，数量是bib_ibi​，其中bi=0b_i=0bi​=0表示数量无限。问对于大小在[1,n][1,n][1,n]之间的每个背包，把背包装满的方案。 n,m≤105n,m\le10^5n,m≤105 分析 如果物品无限的话生成函数就是11−xai\frac{1}{1-x^{a_i}}1−xai​1​ 否则为1−xai(bi+1)1−xai...
LOJ 556: 「Antileaf's Round」咱们去烧菜吧
weixin_30781107的博客
04-27 134
题目传送门：LOJ #556。 题意简述： 有 \(m\) 种物品，第 \(i\) 种物品的体积为 \(a_i\)，数量为 \(b_i\)。 对于 \(1\le i\le n\)，求出使用这些物品恰好填满容积为 \(i\) 的背包的方案数。 两种方案不同当且仅当存在一种物品在两种方案中被选取的数量不同。 \(1\le n,m\le 10^5\)。 题解： 显然可以 \(\mathcal{...
[loj556][Antileaf's Round]咱们去烧菜吧
zhouyuheng2003的博客
02-28 254
前言 感觉挺套路的？ 反正是付公主的背包的严格加强版 题目大意 给出nnn个商品，第iii个商品的体积为viv_ivi​，有bi个 现在给出一个mmm，对背包大小s∈[1,m]s\in[1,m]s∈[1,m]求填满方案数 数据范围 n≤100000,m≤100000n\le100000,m\le100000n≤100000,m≤100000 做法 此处引用[付公主的背包]一题的思路直接作修改了 考...
LOJ #556. 「Antileaf's Round」咱们去烧菜吧
weixin_34239169的博客
12-25 88
好久没更博了 咕咕咕 现在多项式板子的常数巨大...周末好好卡波常吧.... LOJ #556 题意 给定$ m$种物品的出现次数$B_i$以及大小$A_i$ 求装满大小为$[1..n]$的背包的方案数各是多少 数据范围全是$ 10^5$ $ Solution$ 转化成生成函数求解 即是要求 $Ans=\prod\limits_{i=1}^m \sum\limit...
【LOJ556】咱们去烧菜吧（多项式exp）
zxyoi_dreamer的博客（不定期诈尸）
02-10 217
传送门 和这道题是一样的，就不写题解了：付公主的背包 代码： #include<bits/stdc++.h> #define ll long long #define re register #define cs const namespace IO{ inline char gc(){ static cs int Rlen=1<<22|1;static cha...
Loj X:创建逻辑电路-开源
07-04
Loj X 是一个开源项目，专注于提供一个平台，让人们能够创建和设计逻辑电路，从而实现基本的计算功能，甚至构建小型的完整处理器。这个工具集成了各种元件，如逻辑门、存储器和程序读取器，使得用户可以构建出复杂的...
log4j配置，jar包
06-01
总结来说，Log4j是Java开发中不可或缺的工具，通过配置`log4j.properties`文件，我们可以精确地控制日志的输出，便于调试和监控系统运行状态。而`log4j-1.2.11.jar`则包含了实现这些功能所需的类库。
up_717085_google-api-php-client_71loj.rar
12-22
标题中的"up_717085_google-api-php-client_71loj.rar"表明这是一个关于Google API PHP客户端的压缩文件，其中包含了用于与Google API交互的PHP代码示例。这个工具通常用于开发人员在PHP环境中访问和利用Google的...
mybaits 架包
01-30
MyBatis 是一款著名的Java持久层框架，它简化了数据库操作与对象之间的映射工作，提供了灵活的SQL映射和事务管理。MyBatis-3.3.0 是该框架的一个稳定版本，它包含了对SQL的动态构造能力，以及对XML配置和注解方式的...
[多项式exp] LOJ#556. 「Antileaf's Round」咱们去烧菜吧
CE玩家
03-06 537
学一发分治FFT求多项式exp 感觉这种打法应该比牛顿迭代好打多了 #include &lt;cstdio&gt; #include &lt;iostream&gt; #include &lt;algorithm&gt; using namespace std; const int N=1000010,P=998244353; inline void add(int &amp;x,int...
LOJ#556. 「Antileaf's Round」咱们去烧菜吧（牛顿迭代）
DZYO的博客
07-16 543
传送门 题解： 考前练练手感。 G=∏FiG=∏FiG = \prod F_i Fi={1−xai(bi+1)1−xai11−xaibi&gt;0bi=0Fi={1−xai(bi+1)1−xaibi&gt;011−xaibi=0F_i = \begin{cases} \frac{1-x^{a_i(b_i+1)}}{1-x^{a_i}} & b_i>0\\ \frac{1}{1-x^{a...
洛谷P4389 付公主的背包/loj #556. 「Antileaf's Round」咱们去烧菜吧
OZY的博客
12-23 266
学了几天，总算是把loj这个题做出来了。。 随便开了一套题。。然后第一题就不会实在是太没有面子了 有挺多前置知识的 我的学习资料 两个题的题意基本都是一样的 只不过第一个比较简单，可以视为完全没有限制 第二个的话则是有限制的 很显然，第二个是第一个的加强版。。 当然两题都可以用第二个的做法来做 其实 从推导的难度来说，第二个做法简单一些 从知识的要求来说，第一个做法简单一些 首先，当然是先构造生成...
loj556 咱们去烧菜吧 polyExp+调和级数枚举
BPM136(BPMThor)
08-13 192
咋一看就是要算∏i=1n(1−xai(bi+1))∏i=1n(1−xai)\frac{\prod_{i=1}^n (1-x^{a_i(b_i+1)})}{\prod_{i=1}^n (1-x^{a_i})}∏i=1n​(1−xai​)∏i=1n​(1−xai​(bi​+1))​ 大力分治NTT？完蛋发现好像是O(nmlognlogm)O(nmlognlogm)O(nmlognlogm)的 看到一堆...
[多项式ln][多项式exp][背包DP][生成函数] LOJ #556. 咱们去烧菜吧
Vectorxj的博客
03-06 862
SolutionSolutionSolution 就是个背包DP咯。 答案长这个样子ans=[xn]∏i=1m1−xaibi+ai1−xaians=[xn]∏i=1m1−xaibi+ai1−xaians=[x^n]\prod_{i=1}^m{1 - x^{a_ib_i+a_i}\over 1-x^{a_i}}考虑取个lnln\lnlnans=∑i=1m(ln(1−xaibi+ai)−ln(1−...
#4155. 咱们去烧菜吧
Johnny817的博客
04-09 202
题意 内存限制：512 MiB 时间限制：1000 ms 你有 mmm 种物品，第 iii 种物品的大小为 aia_iai​，数量为 bib_ibi​（bi=0b_i=0bi​=0 表示有无限个）。 你还有 nnn 个背包，体积分别为 111 到 nnn，现在你很想知道用这些物品填满某个背包的方案数。 为了满足你的好奇心，你决定把填满每个背包的方案数都算一遍。 因为你其实只是闲得无聊，所以你只想知...
20200630 多项式难题训练
Master.Yi的博客
07-02 2418
LOJ#6247. 九个太阳（单位根反演） 给定 n≤1015,k≤220n\le10^{15},k\le2^{20}n≤1015,k≤220，kkk 是 2 的幂，求： ∑k∣i,0≤i≤n(ni)\sum_{k|i,0\le i\le n}\binom nik∣i,0≤i≤n∑​(in​) 单位根反演 相当于就是求 f(x)=∑i=0n(ni)xi=(x+1)nf(x)=\sum_{i=0}^n \binom nix^i=(x+1)^nf(x)=∑i=0n​(in​)xi=(x+1)n 的 kkk 的倍
[启发式分裂] BZOJ5200 [NWERC2017]Factor-Free Tree
CE玩家
03-07 823
当找到一个根的时候会把区间分成两个部分，然后递归做 问题就是怎么找根了 可以用bzoj4059那种思路，大概叫启发式分裂吧 对每个数，可以预处理一下在它左边第一个与他不互质的数和在它右边第一个与它不互质的数 然后就 O(nlogn)O(nlog⁡n)\mathcal{O}(n\log n) 搞就可以了 #include &lt;cstdio&gt; #include &lt;iost...
Windows 12 网页版
最新发布
weixin_56800842的博客
07-19 98
Windows 12 网页版是一个在线开源项目，由一个初三学生用前端开发三件套 HTML+CSS+JavaScript 实现的一个网页版的 Windows 12，灵感来源于Windows 12 概念版 (Powered by PowerPoint)，直接打开网页即可在线体验Windows12了，希望让用户在网络上预先体验 Windows 12，不过内容可能与 Windows 12 正式版本不一致，精美的ui设计，流畅丰富的动画，各种高级的功能(相较于网页版)。
有源汇有上下界的最小流（LOJ117）
06-11
具体地，对于一条容量为$c$，下界为$l$，上界为$u$的边$(u,v)$，我们在新图中从$u$向$v$连一条容量为$u-l$的边和一条容量为$u-l$，费用为$c$的边，再从源点向$u$连一条容量为$l$的边，费用为$0$，最后从$v$向汇点连...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值