#include <cstdio>
#include <iostream>
using namespace std;
const int M = 3000005;
const int MOD = 1e9+7;
#define int long long
int read()
{
	int x=0,f=1;char c;
	while((c=getchar())<'0' || c>'9') {if(c=='-') f=-1;}
	while(c>='0' && c<='9') {x=(x<<3)+(x<<1)+(c^48);c=getchar();}
	return x*f;
}
int n,m,x,y,ans,fac[M],inv[M];
void init(int n)
{
	inv[0]=inv[1]=fac[0]=1;
	for(int i=1;i<=n;i++) fac[i]=fac[i-1]*i%MOD;
	for(int i=2;i<=n;i++) inv[i]=inv[MOD%i]*(MOD-MOD/i)%MOD;
	for(int i=2;i<=n;i++) inv[i]=inv[i-1]*inv[i]%MOD;
}
int C(int n,int m)
{
	if(m<0 || n<m) return 0;
	return fac[n]*inv[m]%MOD*inv[n-m]%MOD;
}
void flip1(int &x,int &y) {swap(x,y);x--;y++;}
void flip2(int &x,int &y) {swap(x,y);x+=m+2;y-=m+2;}
int cal(int x,int y) {return C(x+y,y);}
signed main()
{
	init(3e6);
	n=read();m=read();x=n+m+1;y=n;
	ans=cal(x,y);
	while(x>=0 && y>=0)
	{
		flip1(x,y);ans-=cal(x,y);
		flip2(x,y);ans+=cal(x,y);
		ans=(ans%MOD+MOD)%MOD;
	}
	x=n+m+1;y=n;
	while(x>=0 && y>=0)
	{
		flip2(x,y);ans-=cal(x,y);
		flip1(x,y);ans+=cal(x,y);
		ans=(ans%MOD+MOD)%MOD;
	}
	printf("%lld\n",ans);
}