显示代码纯文本
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=1e5+10,p=23333333;
typedef long long ll;
int n,dp[2][N],*x=dp[0],*y=dp[1],size,g[N],f[320][N];
int inc(int x,int y){x+=y;return x>=p?x-p:x;}
int dec(int x,int y){x-=y;return x<0?x+p:x;}
int mul(int x,int y){return (ll)x*y%p;}
int main()
{
freopen("get_bag.in","r",stdin);
freopen("get_bag.out","w",stdout);
scanf("%d",&n);
size=sqrt(n);
f[0][0]=1;
for (int i=1;i<=size;i++)
for (int j=size+1;j<=n;j++)
f[i][j]=inc(f[i][j-i],f[i-1][j-size-1]);
for (int i=0;i<=size;i++)
for (int j=0;j<=n;j++)
g[j]=inc(g[j],f[i][j]);
x[0]=1;
for (int i=1;i<=size;i++,swap(x,y))
for (int j=0;j<=n;j++){
y[j]=x[j];
if (j>=i) y[j]=inc(x[j],y[j-i]);
if (j>=i*(i+1)) y[j]=dec(y[j],x[j-i*(i+1)]);
}
int ans=0;
for (int i=0;i<=n;i++) ans=inc(ans,mul(x[i],g[n-i]));
printf("%d\n",ans);
return 0;
}