2020杭电暑假多校第一场

第五题： Fibonacci Sum

我们写出斐波那契的通项公式，然后令a=1+sqrt(5)/2, b=1-sqrt(5)/2，因为5是1e9+9的二次剩余。用x来替代，那么我们a就可以变成(1+x)*inv2，同理b变成(1-x)*inv2。写出替换之后我们二项式展开然后就可以发现当我们r和c固定的时候，C后面就是一串等比数列，所以用等比数列求和公式和欧拉降幂就可以得到答案了。

#include<bits/stdc++.h>
//#define int long long
typedef long long ll;
using namespace std;
const int N=1e5+10,inv2=500000005,p=383008016,mod=1e9+9;
const int a=691504013;
const int b=308495997;
ll n,c,k,inv[N],f[N],mia[N],mib[N];
inline ll qmi(ll a,ll b){
   
    ll ret=1;
    for (;b;b>>=1,a=a*a%mod)
        if (b&1)
            ret=ret*a%mod;
    return (ret+mod)%mod;
}

inline ll C(int n,int m) {
   
    return f[n]*inv[m]%mod*inv[n-m]%mod;
}
ll ph(ll x){
   
    ll res=x,a=x;
    for(ll i=2;i*i<=x;i++){
   
        if(a%i==0){
   
            res=res/i*(i-1);
            while(a%i==0) a/=i;
        }
    }
    if(a>1) res=res/a*(a-1);
    return res;
}
inline ll Inv(ll x){
   
    return qmi(x,mod-2);
}