#include using namespace std; int calc(int a,int b,int mod){ if(a>b){ swap(a,b); } int fac[a+b+1],revfac[b+1],ans=0; long long temp[3]; fac[0]=1; for(int i=1;i<=a+b;++i){ temp[0]=fac[i-1]; temp[0]*=i; temp[0]%=mod; fac[i]=temp[0]; } for(int i=0;i<=b;++i){ temp[0]=fac[i]; temp[1]=mod-2; temp[2]=1; while(temp[1]){ if(temp[1]%2){ temp[2]*=temp[0]; temp[2]%=mod; } temp[0]*=temp[0]; temp[0]%=mod; temp[1]>>=1; } revfac[i]=temp[2]; } for(int i=0;i<=a;++i){ temp[0]=fac[a+b-i]; temp[0]*=revfac[a-i]; temp[0]%=mod; temp[0]*=revfac[b-i]; temp[0]%=mod; temp[0]*=revfac[i]; temp[0]%=mod; ans+=temp[0]; ans%=mod; } return ans; } int main(){ int a,b,mod=1e9+7; long long ans=mod; scanf("%d%d",&a,&b); ans+=calc(a,b,mod); ans-=calc(a-1,b,mod); ans*=(mod/2+1); ans%=mod; printf("%d\n",ans); }