【jzoj5358】【NOIP2017提高A组模拟9.12】【BBQ】
来源:互联网 发布:聊城气象数据 编辑:程序博客网 时间:2024/06/13 11:09
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;}
阅读全文
0 0
- JZOJ5358. 【NOIP2017提高A组模拟9.12】BBQ
- [JZOJ5358]【NOIP2017提高A组模拟9.12】BBQ
- 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
- 【NOIP2017提高A组模拟9.12】Shorten Diameter
- A【NOIP2017提高组模拟12.18】
- 【JZOJ4928】【NOIP2017提高组模拟12.18】A
- 【NOIP2017提高组模拟12.18】A
- 【JZOJ4928】【NOIP2017提高组模拟12.18】A
- 【NOIP2017提高A组模拟7.7】图
- 【NOIP2017提高A组模拟7.13】abcd
- 区间【NOIP2017提高A组模拟7.10】
- 浅谈CountDownLatch的用法
- [日推荐]『Brick4积木』乐高迷不可错过的小工具
- HDU1561 树形DFS+DP
- JS分页
- 支付宝
- 【jzoj5358】【NOIP2017提高A组模拟9.12】【BBQ】
- 线性化与串行化比较
- C++编程题1——输入字符串,输出数字最小的值
- K最近邻算法(介绍)
- MATLAB-轨道分叉
- spring父子容器
- python脚本制作dlib数据集
- Java泛型
- php5.5新特性之yield理解