提交时间:2023-08-16 15:40:46

运行 ID: 98820

#include<bits/stdc++.h> #define int long long using namespace std; const int mod=1e9+7; int qpow(int a,int b) { int ans=1ll; while(b) { if(b&1) ans=(ans*a)%mod; a=(a*a)%mod; b>>=1; } return ans; } const int N=1e5+7; int jie[N],inv[N],s1[N],s2[N]; int C(int n,int k) { return ((jie[n]*inv[k])%mod*inv[n-k])%mod; } const int kl=1000; struct query { int n,k,id; friend bool operator<(query i,query j) { if(i.n/kl==j.n/kl) return i.k<j.k; return i.n<j.n; } }q[N]; int n,t,k; int ans[N],nowans; signed main() { srand(114524); ios::sync_with_stdio(0); cin.tie(0); cout.tie(0); jie[1]=jie[0]=1; for(int i=2;i<N;i++) jie[i]=(jie[i-1]*i)%mod; inv[N-1]=qpow(jie[N-1],mod-2); for(int i=N-2;i>=0;i--) inv[i]=(inv[i+1]*(i+1))%mod; cin>>t; for(int i=1;i<=t;i++) { cin>>q[i].n>>q[i].k; q[i].id=i; } sort(q+1,q+t+1); n=1;k=0; nowans=1; for(int i=1;i<=t;i++) { if(k>q[i].k&&k-q[i].k>N/2) { k=0; nowans=1; } while(n<q[i].n) { nowans=nowans*2-C(n,k); n++; nowans%=mod; } while(k<q[i].k) { k++; nowans=nowans+C(n,k); nowans%=mod; } while(n>q[i].n) { n--; nowans=(nowans+C(n,k))*inv[2]; nowans%=mod; } while(k>q[i].k) { nowans=nowans-C(n,k); nowans%=mod; k--; } ans[q[i].id]=(nowans%mod+mod)%mod; } for(int i=1;i<=t;i++) cout<<ans[i]<<'\n'; }