Bzoj 2813 奇妙的Fibonacci （线性筛）

题目大意：

题目链接

(来源)

解题思路：

在做这道题之前，首先要知道$F_{gcd(i,j)}=gcd(F_i,F_j)(j>i)(0)$

证明如下：

先证明$F_n=F_k\ast F_{n-k+1}+F_{n-k}(1)$

   $F_n=F_{n-1}+F_{n-2}$(斐波那契数列递推式)

   $F_{n-1}=F_{n-2}+F_{n-3}$

   $F_{n-2}=F_{n-3}+F_{n-4}$

 所以$F_n=3\ast F_{n-3}+F_{n-4}$

   $F_n=F_4\ast F_{n-3}+F_{n-4}$

归纳法$F_n=F_k\ast F_{n-k+1}+F_{n-k}$

接下来就是利用数学归纳法证明即可

设$j\%i=k$，所以$j=ti+k$

由$(1)$得$F_j=F_{ti}\ast F_{j-ti+1}+F_{j-ti}$

   $F_j=F_{ti}\ast F_{j-ti+1}+F_k(2)$

 又有$j=ti+k(3)$

可以发现$(2)(3)$的式子是极其相似的，而这正是 辗转相除法的样子 ，所以只要按照类似的方法递推下去，就能证明$(0)$

所以这道题就可以转换为求给定的每个j的因子 有多少个

设$d[i]$代表$i$的约数个数，$num[i]$代表$i$的最小质因子个数

设$f[i]$代表$i$的约数平方, $g[i]$代表$i$的最小质因子$p$(最高次数为$k$)$p^0+...+p^{2k}$的和

因为都是积性函数，所以可以线性筛出来

void init() {
    d[1] = f[1] = g[1] = 1;
    for (ll i = 2; i < maxn; i++) {
        if (!vis[i]) p[++cnt] = i, d[i] = 2, num[i] = 1, f[i] = g[i] = i * i + 1;
        for (ll j = 1; j <= cnt && i * p[j] < maxn; j++) {
            if (i * p[j] >= maxn) break;
            vis[p[j] * i] = true;
            if (i % p[j] == 0) {//p[j]是i的最小质因子，所以p[j]也是i*p[j]的因子
                num[i * p[j]] = (num[i] + 1);//更新num
                d[i * p[j]] = d[i] / num[i * p[j]] * (num[i * p[j]] + 1);//利用num来更新d
                g[i * p[j]] = g[i] * p[j] * p[j] + 1;//更新g
                f[i * p[j]] = f[i] / g[i] * g[i * p[j]];//利用g来更新f
                break;
            }
            else { //i与p[j]互质
                d[i * p[j]] = d[i] * 2;
                num[i * p[j]] = 1;
                f[i * p[j]] = f[i] * f[p[j]];
                g[i * p[j]] = p[j] * p[j] + 1;
            }
        }   
    }
}

到这其实就差不多了，不过最后要注意的一定是
当询问是奇数时，要加上$F_2$的结果

AC代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll maxn = 1e7 + 1e6;
const ll mod = 1e9 + 7;
bool vis[maxn];
ll d[maxn], num[maxn]; //d[i]代表i的约数个数，num[i]代表i的最小质因子的次方
int cnt;
ll p[maxn], f[maxn], g[maxn];//cnt代表质数的个数，p[]存储质数，f[i]约数的平方和，g[i]代表最小质因子的p^0+...+p^2k和（最高次数为k）
void init() {
    d[1] = f[1] = g[1] = 1;
    for (ll i = 2; i < maxn; i++) {
        if (!vis[i]) p[++cnt] = i, d[i] = 2, num[i] = 1, f[i] = g[i] = i * i + 1;
        for (ll j = 1; j <= cnt && i * p[j] < maxn; j++) {
            if (i * p[j] >= maxn) break;
            vis[p[j] * i] = true;
            if (i % p[j] == 0) {
                num[i * p[j]] = (num[i] + 1);
                d[i * p[j]] = d[i] / num[i * p[j]] * (num[i * p[j]] + 1);
                g[i * p[j]] = g[i] * p[j] * p[j] + 1;
                f[i * p[j]] = f[i] / g[i] * g[i * p[j]];
                break;
            }
            else {
                d[i * p[j]] = d[i] * 2;
                num[i * p[j]] = 1;
                f[i * p[j]] = f[i] * f[p[j]];
                g[i * p[j]] = p[j] * p[j] + 1;
            }
        }   
    }
}
ll Q, F, A, B, C;
ll ans1 = 0, ans2 = 0;
int main() {
    init();
    cin >> Q >> F >> A >> B >> C;
    while (Q--) {
        ans1 = (ans1 + d[F]) % mod;
        ans2 = (ans2 + f[F]) % mod;
        F = (F * A % C + B) % C + 1;
        if (F & 1) ans1 = (ans1 + 1) % mod, ans2 = (ans2 + 4) % mod;
    }
    cout << ans1 << endl << ans2 << endl;
    return 0;
}