很容易想到二维DP暴力吧,但是题目不允许,于是
开二维dp(i,j)表示堆到i层且r用了j个的方案数,
如果i*(i+1)/2<=g+j,dp(i,j)+=dp(i-1,j);表示这层如果不放r的话
如果j大于等于i,dp(i,j)+=dp(i-1,j-i);
那么压缩下空间应该不难,初始化dp(0)=1
#include<bits/stdc++.h>
using namespace std;
#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long
#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
#define pii pair<int,int>
#define fi first
#define se second
#define mk(x,y) make_pair(x,y)
const int mod=1e9+7;
const int maxn=2e5+100;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){
if(y==0) return x;
return gcd(y,x%y);
}
int r,g;
int dp[maxn];
int main(){
ios::sync_with_stdio(false);
cin>>r>>g;
if(r>g) swap(r,g);
int maxh=(int)sqrt(2.0*(r+g));
dp[0]=1;
int ans1=0,ans2=0,ret=0;
rep(i,1,maxh+1){
ret=0;
for(int j=r;j>=0;j--){
if(i*(i+1)/2>g+j) dp[j]=0;
if(j>=i) dp[j]=(dp[j]+dp[j-i])%mod;
///cout<<dp[j]<<" ";
ret=(ret+dp[j])%mod;
}
if(i>ans1&&ret) ans1=i,ans2=ret;
}
cout<<ans2<<" ";
return 0;
}