【jzoj5358】【NOIP2017提高A组模拟9.12】【BBQ】
2017-09-16 11:07
399 查看
description
solution
可以发现组合数的意义就是在二维平面上一个点只能往上往右走走到另一个点的方案数,这个问题可以用递推来解决,相当于在(-a[i],-b[i])上都加一,求f(a[i],b[i])的和,减去i,j相等的情况再除以二即可。
code
#include<cstdio> #include<cmath> #include<cstring> #include<algorithm> #define LF double #define LL long long #define ULL unsigned int #define fo(i,j,k) for(int i=j;i<=k;i++) #define fd(i,j,k) for(int i=j;i>=k;i--) #define fr(i,j) for(int i=begin[j];i;i=next[i]) using namespace std; int const mn=2*1e5+3,mm=8*1e3,mo=1e9+7; int n,a[mn],b[mn],fact[mm+3],ni[mm+3],f[4000+5][4000+5]; int c(int x,int y){ return 1ll*fact[x]*ni[y]%mo*ni[x-y]%mo; } int Pow(LL x,int y){ LL z=1; while(y){ if(y&1)z=z*x%mo; x=x*x%mo; y>>=1; } return z; } int main(){ //freopen("a.in","r",stdin); //freopen("a.out","w",stdout); freopen("d.in","r",stdin); freopen("d.out","w",stdout); scanf("%d",&n); int bit=2*1e3+1; fo(i,1,n)scanf("%d%d",&a[i],&b[i]),f[bit-a[i]][bit-b[i]]++; fact[0]=1; fo(i,1,mm)fact[i]=1ll*fact[i-1]*i%mo; ni[mm]=Pow(fact[mm],mo-2); fd(i,mm-1,0)ni[i]=1ll*ni[i+1]*(i+1)%mo; int ans=0; fo(i,1,bit*2)fo(j,1,bit*2)f[i][j]=(1ll*f[i][j]+f[i-1][j]+f[i][j-1])%mo; fo(i,1,n)ans=(ans+f[bit+a[i]][bit+b[i]])%mo; fo(i,1,n)ans=(ans-c(a[i]*2+b[i]*2,a[i]*2))%mo; printf("%d",(1ll*ans*Pow(2,mo-2)%mo+mo)%mo); return 0; }
相关文章推荐
- JZOJ5358. 【NOIP2017提高A组模拟9.12】BBQ
- [JZOJ5358]【NOIP2017提高A组模拟9.12】BBQ
- JZOJ5358. 【NOIP2017提高A组模拟9.12】BBQ
- 【JZOJ 5358】【NOIP2017提高A组模拟9.12】BBQ
- 5358. 【NOIP2017提高A组模拟9.12】BBQ
- 【jzoj5360】【NOIP2017提高A组模拟9.12】【Shorten Diameter】
- 【jzoj5359】【NOIP2017提高A组模拟9.12】【Arrays and Palindrome】
- JZOJ5360. 【NOIP2017提高A组模拟9.12】Shorten Diameter
- JZOJ5358【NOIP2017提高A组模拟9.12】BBQ
- [JZOJ5391]【NOIP2017提高A组模拟9.26】卡常题
- 【jzoj5338】【NOIP2017提高A组模拟8.25】【影子】【点分治】
- 【jzoj5346】【NOIP2017提高A组模拟9.5】【NYG的背包】【贪心】
- JZOJ4919. 【NOIP2017提高组模拟12.10】神炎皇
- 【JZOJ4920】【NOIP2017提高组模拟12.10】降雷皇
- 【JZOJ4928】【NOIP2017提高组模拟12.18】A
- JZOJ 5393. 【NOIP2017提高A组模拟10.5】Snake vs Block
- jzoj5336 【NOIP2017提高A组模拟8.24】提米树 (dfs序dp,奇异姿势dp)
- [JZOJ5395]【NOIP2017提高A组模拟10.6】Count
- JZOJ 5398. 【NOIP2017提高A组模拟10.7】Adore
- jzoj5400. 【NOIP2017提高A组模拟10.7】Repulsed