HDU 6333.Harvest of Apples

Problem

有若干次询问，每次询问给出\(n, m ( 1\le n, m \le 10^5)\)，询问\[\sum_{i=0}^{m}C_n^i\]

Solution

我们定义\[T_n^m=\sum_{i=0}^mC_n^i\]显然有\[T_n^m=2T_{n-1}^m-C_{n-1}^m\]\[T_n^m=T_n^{m-1}+C_n^m\]然后就可以在已知一个\(T_n^m\)的情况下推知其他的\(T_n^m\)了，为了降低复杂度可以使用莫队算法或者分块打表。

更一般的形如式子\[\sum_{i=0}^m{f(i)*g(n-i)}\]这种卷积可以使用分块FFT求解。

（赛场上找到了递推式没有想到用莫队真是愚蠢啊，看来过往对莫队算法的理解还是不够）

Code

/*
 Author: LargeDumpling
 Email: LargeDumpling@qq.com
 Edit History:
    2018-08-01  File created.
*/

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int MAXN=100050;
const long long mod=1000000007;
const int BLOCK=316;
struct jz
{
    int id;
    long long n,m,ans;
    bool operator<(const jz &X)const { return n/BLOCK==X.n/BLOCK?m<X.m:n<X.n; }
}Q[MAXN];
long long fac[MAXN],ifac[MAXN];
long long pOw(long long a,long long n)
{
    long long ans;
    for(ans=1;n;n>>=1,a=a*a%mod)
        if(n&1)
            ans=ans*a%mod;
    return ans;
}
bool cmp(jz A,jz B) { return A.id<B.id; }
void pre_calc()
{
    ifac[0]=fac[0]=1;
    for(int i=1;i<MAXN;i++)
        fac[i]=fac[i-1]*i%mod;
    ifac[MAXN-1]=pOw(fac[MAXN-1],mod-2);
    for(int i=MAXN-2;i;i--)
        ifac[i]=ifac[i+1]*(i+1)%mod;
    return;
}
long long C(long long n,long long m) { return fac[n]*ifac[m]%mod*ifac[n-m]%mod; }
int main()
{
    int T_T;
    long long curn=1,curm=0,curans=1,idu=pOw(2,mod-2);
    pre_calc();
    scanf("%d",&T_T);
    for(int i=1;i<=T_T;i++)
    {
        scanf("%lld%lld",&Q[i].n,&Q[i].m);
        Q[i].id=i;
    }
    sort(Q+1,Q+T_T+1);
    for(int i=1;i<=T_T;i++)
    {
        while(curn<Q[i].n)
        {
            curans=(curans*2LL-C(curn,curm)+mod)%mod;
            curn++;
        }
        while(Q[i].n<curn)
        {
            curans=(curans+C(curn-1,curm))%mod*idu%mod;
            curn--;
        }
        while(curm<Q[i].m)
        {
            curans=(curans+C(curn,curm+1))%mod;
            curm++;
        }
        while(Q[i].m<curm)
        {
            curans=(curans-C(curn,curm)+mod)%mod;
            curm--;
        }
        Q[i].ans=curans;
    }
    sort(Q+1,Q+T_T+1,cmp);
    for(int i=1;i<=T_T;i++)
        printf("%lld\n",Q[i].ans);
    fclose(stdin);
    fclose(stdout);
    return 0;
}

转载于:https://www.cnblogs.com/LargeDumpling/p/9403481.html