1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85
|
#include<bits/stdc++.h> using namespace std; #define debug(x) cerr << #x << ": " << (x) << endl #define rep(i,a,b) for(int i=(a);i<=(b);i++) #define dwn(i,a,b) for(int i=(a);i>=(b);i--) #define SZ(a) ((int) (a).size()) #define pb push_back #define all(x) (x).begin(), (x).end() #define x first #define y second #define int long long using pii = pair<int, int>; using ll = long long; inline void read(int &x){ int s=0; x=1; char ch=getchar(); while(ch<'0' || ch>'9') {if(ch=='-')x=-1;ch=getchar();} while(ch>='0' && ch<='9') s=(s<<3)+(s<<1)+ch-'0',ch=getchar(); x*=s; }
const int N=2e5+5, mod=1e9+7;
int add(int x, int y){ return (x+y)%mod; }
int mul(int x, int y){ return 1LL*x*y%mod; }
int fpow(int x, int p){ int res=1; for(; p; p>>=1, x=1LL*x*x%mod) if(p&1) res=1LL*res*x%mod; return res%mod; }
int inv(int x){ return fpow(x, mod-2); }
int fac[N]; int invfac[N];
void getFac(int n){ fac[0]=invfac[0]=1; for(int i=1; i<=n; i++) fac[i]=1LL*fac[i-1]*i%mod, invfac[i]=1LL*invfac[i-1]*inv(i)%mod; }
int C(int a, int b){ if(b<0 || b>a) return 0; return 1LL*fac[a]*invfac[b]%mod*invfac[a-b]%mod; }
int cal(int x, int y){ return C(x+y, y); }
void solve(){ getFac(N-1); int n, m; cin>>n>>m; int a, b; cin>>a>>b; int res=0; rep(j, b+1, m){ res=add(res, mul(cal(n-a-1, j-1), cal(a-1, m-j))); } cout<<res<<"\n"; }
signed main(){ solve(); return 0; }
|