[Luogu P3704] [BZOJ 4816] [SDOI2017]数字表格

洛谷传送门

BZOJ传送门

题目描述

Doris刚刚学习了fibonacci数列。用$f[i]$表示数列的第$i$项，那么

$f[0]=0,f[1]=1$,

$f[n]=f[n-1]+f[n-2],n\ge 2$

Doris用老师的超级计算机生成了一个$n\times m$的表格，

第$i$行第$j$列的格子中的数是$f[\gcd (i,j)]$，其中$\gcd (i,j)$表示$i,j$的最大公约数。

Doris的表格中共有$n\times m$个数，她想知道这些数的乘积是多少。

答案对$10^9+7$取模。

输入输出格式

输入格式：

有多组测试数据。

第一个一个数$T$，表示数据组数。

接下来$T$行，每行两个数$n,m$

输出格式：

输出$T$行，第$i$行的数是第$i$组数据的结果

输入输出样例

输入样例#1：

3
2 3
4 5
6 7

输出样例#1：

1
6
960

说明

对$10\%$的数据，$1\le n,m\le 100$

对$30\%$的数据，$1\le n,m$

另外存在$30\%$的数据，$T\le 3$

对$100\%$的数据，$T\le 1000,1\le n,m\le 10^6$

解题分析

$$\begin{gathered} \prod _{i=1}^n\prod _{j=1}^{m}F[gcd(i,j)] \\ =\prod _{d=1}^{n}F[d{]}^{{\sum }_{i=1}^n{\sum }_{j=1}^m[gcd(i,j)=d]} \\ =\prod _{d=1}^{n}F[d{]}^{{\sum }_{p=1}^{\lfloor \frac{n}{d}\rfloor }\mu (p)\lfloor \frac{n}{dp}\rfloor \lfloor \frac{m}{dp}\rfloor } \end{gathered}$$

套路地设$T=dp$
$$\begin{gathered} \prod _{d=1}^{n}F[d{]}^{{\sum }_{p=1}^{\lfloor \frac{n}{d}\rfloor }\mu (p)\lfloor \frac{n}{dp}\rfloor \lfloor \frac{m}{dp}\rfloor } \\ =\prod _{T=1}^n\prod _{d|T}F[d{]}^{\mu (\frac{T}{d})\lfloor \frac{n}{T}\rfloor \lfloor \frac{m}{T}\rfloor } \\ =\prod _{T=1}^n\prod _{d|T}(F[d{]}^{\mu (\frac{T}{d})}{)}^{\lfloor \frac{n}{T}\rfloor \lfloor \frac{m}{T}\rfloor } \end{gathered}$$
然后就可以筛出$\mu (i)$后$O(Nln(N))$预处理， $O(T\sqrt{N}log(N))$回答询问了。

代码如下：

#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <cctype>
#include <algorithm>
#define R register
#define IN inline
#define W while
#define gc getchar()
#define MOD 1000000007ll
#define ll long long 
#define MX 1000500
template <class T>
IN void in(T &x)
{
    x = 0; R char c = gc;
    for (; !isdigit(c); c = gc);
    for (;  isdigit(c); c = gc)
    x = (x << 1) + (x << 3) + c - 48;
}
template <class T> IN T max(T a, T b) {return a > b ? a : b;}
template <class T> IN T min(T a, T b) {return a < b ? a : b;}
IN void add(R int ad, int &tar) {tar += ad; if (tar > MOD) tar -= MOD;}
int pcnt;
int miu[MX], pri[MX];
ll sum[MX], F[MX], inv[MX];
bool npr[MX];
IN ll fpow(R ll base, R ll tim)
{
    ll ret = 1;
    W (tim)
    {
        if (tim & 1) ret = ret * base % MOD;
        base = base * base % MOD, tim >>= 1;
    }
    return ret;
}
void get()
{
    miu[1] = 1; R int tar;
    for (R int i = 2; i <= 1e6; ++i)
    {
        if (!npr[i]) pri[++pcnt] = i, miu[i] = -1;
        for (R int j = 1; j <= pcnt; ++j)
        {
            tar = i * pri[j];
            if (tar > 1e6) break;
            npr[tar] = true;
            if (!(i % pri[j])) {miu[tar] = 0; break;}
            miu[tar] = -miu[i];
        }
    }
    sum[0] = F[1] = 1;
    for (R int i = 2; i <= 1e6; ++i) F[i] = (F[i - 1] + F[i - 2]) % MOD;
    for (R int i = 1; i <= 1e6; ++i) inv[i] = fpow(F[i], MOD - 2), sum[i] = 1;
    for (R int i = 1; i <= 1e6; ++i)
    for (R int j = 1, tar; (tar = j * i) <= 1e6; ++j)
    {
        if (miu[j] == 1) sum[tar] = sum[tar] * F[i] % MOD;
        else if (miu[j] == -1) sum[tar] = sum[tar] * inv[i] % MOD;
    }
    for (R int i = 2; i <= 1e6; ++i) sum[i] = sum[i - 1] * sum[i] % MOD;
}
int main(void)
{
    ll ans;
    int T, n, m; in(T); get();
    W (T--)
    {
        in(n), in(m); ans = 1;
        if (n > m) std::swap(n, m);
        for (R int lef = 1, rig; lef <= n; lef = rig + 1)
        {
            rig = std::min(n / (n / lef), m / (m / lef));
            ans = ans * fpow(sum[rig] * fpow(sum[lef - 1], MOD - 2) % MOD, 1ll * (n / lef) * (m / lef)) % MOD;
        }
        printf("%lld\n", ans);
    }
}