牛客NOIPtg day5 B-demo的gcd

一句话题意：给定长度为n的序列，求任意两两之间gcd的积mod 998244353的值。

好像是莫比乌斯反演板子题？？？（反正noip估计不考这种毒瘤

考场上想到一个类似正解的思路 好像摊下来最多处理nlogn次就义无反顾地写了结果爆零

（你以为for while if continue 都不要时间的啊

后来一看发现自己的思路貌似毫无问题，就是实现丑了（逃

可以看到对于每一个质数p，对答案的贡献为p^k(k-1) 其中k=序列中包含质因数p的数个数

然后对于质数的幂次p^q，考场上的思路是每次计算质数之后将原序列中的该质数除掉，如此循环直到找不到该质因数为止（于是就T了还不如暴力qvq

其实可以考虑另一种很简单的思路 对于任意一个p^k，枚举他的倍数，对答案的贡献即为p^k'(k'-1)，与上述等价，复杂度为O（nlogn）（其实是ln

代码如下

#include <cstdio>
#include <iostream>
#define qvq register
#define int long long
using namespace std;
const int mod=;
const int maxn=;
int p[maxn],cnt,a[maxn];
bool np[maxn];

inline int read() {
    int x=,f=;
    char cr=getchar();
    while (cr>'' || cr<'') {
        if (cr=='-') f=-;
        cr=getchar();
    }
    while (cr>='' && cr<='') {
        x=(x<<)+(x<<)+cr-'';
        cr=getchar();
    }
    return x*f;
}

inline void euler() {
    for (qvq int i=;i<=;i++) {
        if (!np[i]) p[++cnt]=i;
        for (qvq int j=;j<=cnt && p[j]*i<=;j++) {
            np[p[j]*i]=;
            if (i%p[j]==) break;
        }
    }
}

inline int val(int k) {
    return k*(k-)/;
}

inline int power(int a,int b,int p) {
    if (b==) return ;
    if (b==) return a;
    if (b&) return a*power(a*a%p,b>>,p)%p;
    else return power(a*a%p,b>>,p);
}

int num[maxn],f[maxn];

signed main() {
    int n=read();
    int lim=-;
    for (qvq int i=;i<=n;i++) a[i]=read(),num[a[i]]++,lim=max(lim,a[i]);
    euler();
    for (int i=;i<=cnt;i++) {
        int temp=p[i];
        while (temp<=) {
            f[temp]=p[i];
            temp*=p[i];
        }
    }
    int ans=;
    for (int i=;i<=lim;i++) {
        if (!f[i]) continue;
        int res=;
        for (int j=i;j<=lim;j+=i) res+=num[j];
        ans*=power(f[i],val(res),mod),ans%=mod;
    }
    printf("%lld",ans);
}

我太菜了orz