第五题: 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);
}