Problem: fireworks

#include <cstdio>#include <cstring>#include <algorithm>#include <cmath>#include <cstdlib>#include <cctype>#define Clear(a, b) memset(a, b, sizeof(a))#define In() freopen("in.txt", "r", stdin)#define Out() freopen("out.txt", "w", stdout)#define For(l, h) for(int i = l; i < h; i++)typedef long long ll;using namespace std;const int mod = 1e9 + 7;const ll INF = 1e18;const int maxn = 1e5+5;//-----------------------------------------------------------------------//ll f[maxn] = {1};void fac(int n) {    for(int i = 1; i <= n; i++) {        f[i] = f[i-1] * i;        f[i] %= mod;    }}ll qpow(ll x, ll n) {    ll ret = 1;    while(n) {        if(n&1) ret = ret*x%mod;        x = x*x%mod;        n >>= 1;    }    return ret;}ll C(int m, int n) {    ll a = f[n], b = qpow(f[m], mod-2), c = qpow(f[n-m], mod-2);//printf("%lld %lld %lld\n", a, b, c);    return a*b%mod*c%mod;}int main() {    int n, t, w;    while(scanf("%d%d%d", &n, &t, &w) == 3) {        fac(t);        ll ans = 0;        for(int i = 0; i < n; i++) {            int p, k;            ll s;            scanf("%d%d", &p, &k);            int x = t-abs(w-p); //printf("x = %d\n", x);            if(x < 0 || x&1) s = 0;            else {                x >>= 1;                s = C(x, t);            }            ans += s*k;            ans %= mod;        }        printf("%lld\n", ans);    }    return 0;}/*在纸上写一写发现是个杨辉三角形，杨辉三角形的性质：1.第n行有n项。2.第n行的数字之和为2的n-1次。3.第n行的m个数为C(n-1, m-1)。这题转化为求C(n, m);]，由于n过大，直接计算阶乘超时，所以先打表计算出n的阶乘。对组合数C(n, m) = n!/m!*(n-m!)，由于求模运算对除法不可分配，故采用求逆元的方法来变除为乘。对于正整数和，如果有 a*b % m ≡ 1//(ax ≡ 1(mod m))，那么把这个同余方程中x的最小正整数解叫做模的逆元。感觉有点像倒数是不是。如果m是素数那么由费马小定理：a^p-1%p≡1, a^p-2*a%p≡1,即a^p-2为a的逆元。*/

/*顺便查了一下memset函数的使用：1.memset函数按字节对内存块进行初始化，只能将int数组初始化为0与-1。2.无论多大的数，只有后八位一个字节有效，例如memset(-1)与memset(511),他们的二进制码分别为‘11111111 11111111 11111111 11111111’与‘00000000 00000000 00000001 11111111’最后8位都为1，于是赋值时数组的每个数被初始化为‘11111111 11111111 11111111 11111111’，就是-1。3.举个反例，想把int数组置为1，就会出现00000001000000010000000100000001，结果等于16843009。4.另外用memset初始化结构体：memset(&st, 0, sizeo(structsample_struct));初始化结构体数组：memset(TEST,0,sizeof(structsample_struct)*10），如果结构体中有数组的话还是需要对数组单独进行初始化处理的。*/