第八届山东省赛 C sdut 3895 fireworks（组合数+逆元）

题意：

一个烟花（位置为x）每秒爆炸并往两边（x-1,x+1）分裂成两部分，不断循环。给你n个位置，每个位置上有c个烟花，问t秒后，w位置上有多少个烟花。

解题思路：

容易推出对于一个烟花来说，t秒后周围的烟花数量分部就是杨辉三角的t+1行，不过每两个数之间要多一个0.

n的上限1e5，不能打出杨辉三角需要o（n^2）的时间显然是不行的。

所以我们只能直接求组合数，求的时候除法取模需要求一下逆元。但如果每除一次都算一次逆元，仍然会超时。我们可以直接求出1-m的乘积的逆元，1-m累乘取模，再求逆元，求出来的逆元与1-m累乘求的逆元虽然可能不一样，但是都是有效的逆元。然后我们预处理一下1-50000的阶乘的值即可。

代码：

#include <bits/stdc++.h>
#define LL long long 
using namespace std;
const LL  maxn=1e5+5;
const int mod=1e9+7;
LL a[maxn+5];
long long pw(int a,int b)
{
    long long sum=1,base=a;
    while(b)
    {
        if(b&1)sum=(sum*base)%mod;
        base=(base*base)%mod;
        b>>=1;
    }
    return sum;
}


void init()
{
    LL i, j;
    a[0]=1;
    for(i=1; i<=maxn-3; i++)
    {
        a[i]=(a[i-1]*i)%mod;

    }
    
}
long long C(int m,int n)
{
    long long sum=a[n];
    sum=sum*pw(a[n-m],mod-2)%mod;
    sum=sum*pw(a[m],mod-2)%mod;
    return sum;
}

int main()
{
    int n, t, w, i, j, x;
    LL c;
    init();
    while(~scanf("%d%d%d", &n, &t, &w))
    {
        LL ans=0;
        for(i=1; i<=n; i++)
        {
            scanf("%d%lld", &x, &c);
            if(abs(x-w)%2==t%2 && abs(x-w)<=t)
            {
                ans=(ans+(c*C(abs(t/2-(abs(x-w))/2), t))%mod)%mod;                

          //      printf("%d %d %lld\n", abs(t/2- abs(x-w)/2), t, C(abs(t/2-(abs(x-w))/2), t));
            }
        }
        printf("%lld\n", ans%mod);
    }
}