BZOJ 3288: Mato矩阵

最新推荐文章于 2020-01-13 20:23:43 发布

aisen1985

最新推荐文章于 2020-01-13 20:23:43 发布

阅读量107

点赞数

文章标签： c/c++

原文链接：http://www.cnblogs.com/beiyuoi/p/6200470.html

版权

Description

一个 $n*n$ 行列式,$(i,j)=gcd(i,j)$

Sol

线性筛.

这道题神奇的筛出来 $phi$ ...

打表可以发现,一个数会被他所有的因子减掉因子的 $phi$ ...

然后我就不会证明了...

Code

/**************************************************************
    Problem: 3288
    User: BeiYu
    Language: C++
    Result: Accepted
    Time:432 ms
    Memory:9100 kb
****************************************************************/
 
#include <bits/stdc++.h>
using namespace std;
 
typedef long long LL;
const int N = 1e6+50;
const LL p = 1e9+7;
 
int n;LL ans;
int pr[N],cp;
int phi[N];
 
void Pre() {
    phi[1]=1,ans=1;
    for(int i=2;i<=n;i++) {
        if(!phi[i]) phi[i]=i-1,pr[++cp]=i;
        for(int j=1;j<=cp && (LL)pr[j]*i<=n;j++) {
            if(i%pr[j]) {
                phi[i*pr[j]]=phi[i]*(pr[j]-1);
            }else {
                phi[i*pr[j]]=phi[i]*pr[j];
            }
        }
    }
//  for(int i=1;i<=n;i++) cout<<phi[i]<<" ";cout<<endl;
    for(int i=1;i<=n;i++) ans=(ans*phi[i])%p;
}
 
int main() {
    cin>>n;
    Pre();
    cout<<ans<<endl;
    return 0;
}

转载于:https://www.cnblogs.com/beiyuoi/p/6200470.html