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),如果结构体中有数组的话还是需要对数组单独进行初始化处理的。
*/
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值