很容易想到二维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;
}