當(dāng)時(shí)比賽時(shí)靈機(jī)一動(dòng) 把gigi?1喂給了OEIS
然后就找到了 233 就是這個(gè)咯
然后就水過(guò)去了
題解？題解我還沒(méi)看 先挖個(gè)坑

UPD：跟這個(gè)題是一毛一樣的咯

#include<cstdio> #include<cstdlib> #include<algorithm> using namespace std; typedef long long ll;const int N=1000005; int P; inline ll Pow(ll a,int b){ll ret=1;for (;b;b>>=1,a=a*a%P)if (b&1)ret=ret*a%P;return ret; }int n; ll f[N],g[N];#define read(x) scanf("%d",&(x))int main(){int Q;freopen("t.in","r",stdin);freopen("t.out","w",stdout);read(Q);while (Q--){read(n); read(P);f[0]=0; f[1]=g[1]=1;for (int i=2;i<=n;i++) f[i]=((f[i-1]<<1)%P+f[i-2])%P,g[i]=f[i];for (int i=1;i<=n;i++){ll inv=Pow(g[i],P-2);for (int j=i+i;j<=n;j+=i)g[j]=g[j]*inv%P;}ll lcm=1,ans=0;for (int i=1;i<=n;i++)lcm=lcm*g[i]%P,ans+=lcm*i%P;printf("%lld\n",ans%P);}return 0; }

總結(jié)

以上是生活随笔為你收集整理的[数论 反演] BZOJ 4833 最小公倍佩尔数的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。