BZOJ4833: [Lydsy1704月赛]最小公倍佩尔数

可以看一下这篇东西

---

推一下f(n),e(n)的柿子，有
$(1+\sqrt 2)^{n+1}$
$=(1+\sqrt 2)^n(1+\sqrt 2)$
$=(e(n)+f(n)\sqrt 2)(1+\sqrt 2)$
$=e(n)+2f(n)+(f(n)+e(n))\sqrt 2$
于是$e(n+1)=e(n)+2f(n),f(n+1)=f(n)+e(n)$
即$f(n+1)=2f(n)+f(n-1)$

然后类比斐波那契的性质：$(\ fib[a],fib[b]\ )=fib[(a,b)]$
长得像$f(n+1)=af(n)+bf(n-1)$ 的东西很多都有这个性质？
反正这题有= =

---

那么做个lcm，gcd的容斥，令全集$S=1,2,3....n$有
$g(n)=lcm(f(i),i\in S)$

$$g(n)=\prod_{T\subseteq S,T\neq \phi}f(gcd(T))^{(-1)^{|T|-1}}$$

---

我们定义一个$h(n)$，满足$f(i)=\prod_{d|i}h(d)$，那么$h(i)=\prod_{d|i}f(d)^{μ(i/d)}$

那么

$g(n)=\prod_{T\subseteq S,T\neq \phi}(\ \prod_{d|gcd(T)}\ h(d))^{(-1)^{|T|-1}}$

$g(n)=\prod_dh(d)^{\sum_{T\subseteq S,T\neq \phi,d|gcd(T)}(-1)^{|T|-1}}$

然后观察$h(d)$的指数$\sum_{T\subseteq S,T\neq \phi,d|gcd(T)}(-1)^{|T|-1}$ 这个东西
因为$T\neq \phi$，于是当$d\in S$时值为1，否则为0

于是有$g(n)=\prod_{i=1}^nh(i)$

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int maxn = 1e6+10;

int p[maxn],pri,mu[maxn];
bool v[maxn];
void pre()
{
    mu[1]=1;
    for(int i=2;i<maxn;i++)
    {
        if(!v[i]) p[++pri]=i,mu[i]=-1;
        for(int j=1,k=p[j]*i;k<maxn;j++,k=p[j]*i)
        {
            v[k]=true;
            if(i%p[j]==0) break;
            mu[k]=-mu[i];
        }
    }
}

int mod,n;
int pw(int x,int k)
{
    int re=1;
    for(;k;k>>=1,x=(ll)x*x%mod) if(k&1)
        re=(ll)re*x%mod;
    return re;
}
int inv(int x){ return pw(x,mod-2); }
int f[maxn],invf[maxn],g[maxn];

int main()
{
    //freopen("tmp.in","r",stdin);
    //freopen("tmp.out","w",stdout);

    // fi=2fi-1+fi-2
    pre();

    int T; scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&mod);
        f[0]=0,f[1]=1; for(int i=2;i<=n;i++) f[i]=(2ll*f[i-1]+f[i-2])%mod;
        for(int i=1;i<=n;i++) invf[i]=inv(f[i]);

        for(int i=1;i<=n;i++) g[i]=1;
        for(int i=1;i<=n;i++)
        {
            for(int j=1;i*j<=n;j++) if(mu[j])
                g[i*j]=(ll)g[i*j]*(mu[j]==1?f[i]:invf[i])%mod;
        }
        int ans=0;
        for(int i=1,t=1;i<=n;i++)
        {
            t=(ll)t*g[i]%mod;
            ans=(ans+(ll)t*i%mod)%mod;
        }
        ans=(ans+mod)%mod;
        printf("%d\n",ans);
    }


    return 0;
}