#include<bits/stdc++.h>
using namespace std;

const int mod=1e9+7;
int pls(int x,int y){return (x=(x+y)%mod)<0?x+mod:x;}
int mul(int x,int y){return pls(1ll*x*y%mod,0);}
int power(int x,int y=mod-2){
    int c=1;
    while(y){if(y&1){c=mul(c,x);}x=mul(x,x),y>>=1;}
    return c;
}

int main(){
    int t,Max=0;scanf("%d",&t);
    vector< pair<int,int> > a(t);
    for(pair<int,int>& i:a){
        scanf("%d%d",&i.first,&i.second);
        Max=max(Max,max(i.first,i.second));
    }Max<<=1;
    vector<int> fac(Max+1);
    fac[0]=1;
    for(int i=1;i<=Max;i++) fac[i]=mul(fac[i-1],i);
    #define C(N,M) mul(fac[N],mul(power(fac[M]),power(fac[N-(M)])))
    for(int i=0,n,m;i<t;i++){
        n=a[i].first,m=a[i].second;
        printf("%d\n",mul(mul(n,m)-1,C(n+m-2,n-1)));
    }
    return 0;
}