#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7,inv2=(mod+1>>1);
void add(int &x,const int y){
	if((x+=y)>=mod)x-=mod; 
}
int mul(int x,int y){
	return 1ll*x*y%mod;
}
int n,m;
#define N 3005
int f[N],g[N][N],C[N*2][N*2],ipw[N*2],fac[N*2];
int slv(int p,int q){
	int s=0,ct=m+q;
	for(int i=0;i<=min(p,q);++i){
		int s1=mul(mul(C[p][i],C[q][i]),fac[i]);
		int ct1=ct-2*i;
		s1=mul(s1,mul(C[ct1][n-p],fac[n-p]));
		if((ct1-n+p)&1)continue;
		s1=mul(s1,f[ct1-n+p>>1]);
		if(i&1)add(s,mod-s1);
		else add(s,s1);
	}
	return mul(s,ipw[q]);
}
int main(){
	ios::sync_with_stdio(0);
	cin>>n>>m,ipw[0]=fac[0]=1;
	for(int i=1;i<=(N-5<<1);++i)fac[i]=mul(i,fac[i-1]);
	for(int i=0;i<=(N-5<<1);++i){
		C[i][0]=1;if(i)ipw[i]=mul(inv2,ipw[i-1]);
		for(int j=1;j<=i;++j)C[i][j]=C[i-1][j-1],add(C[i][j],C[i-1][j]);
	}
	f[0]=1;int ans=0;
	for(int i=1;i<=m;++i)f[i]=mul(f[i-1],C[i*2][2]);
	for(int i=0;i<=n;++i)
		for(int j=0;j<=m;++j)
			if(n+i==m+j)add(ans,mul(C[n][i],mul(C[m][j],slv(i,j))));
	cout<<ans;
	return 0;
}