题意
Mato同学最近正在研究一种矩阵,这种矩阵有n行n列第i行第j列的数为gcd(i,j)。
例如n=5时,矩阵如下:
1 1 1 1 1
1 2 1 2 1
1 1 3 1 1
1 2 1 4 1
1 1 1 1 5
Mato想知道这个矩阵的行列式的值,你能求出来吗?
答案模10^9+7
n<=1000000
分析
一开始想到可以先根据高斯消元推出来把该矩阵消成下三角矩阵后对角线上的元素是啥,结果推了半天发现推不出来。
看了题解,发现原来打表就可以发现(i,i)的元素就是
φ(i)
。。。
那么就直接求
φ(1)到φ(n)
的乘积即可。
代码
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const int N=1000005;
const int MOD=1000000007;
int n,tot,phi[N],prime[N];
bool not_prime[N];
void get_prime(int n)
{
phi[1]=1;
for (int i=2;i<=n;i++)
{
if (!not_prime[i]) prime[++tot]=i,phi[i]=i-1;
for (int j=1;j<=tot&&i*prime[j]<=n;j++)
{
not_prime[i*prime[j]]=1;
if (i%prime[j]==0)
{
phi[i*prime[j]]=phi[i]*prime[j];
break;
}
phi[i*prime[j]]=phi[i]*(prime[j]-1);
}
}
}
int main()
{
scanf("%d",&n);
get_prime(n);
int ans=1;
for (int i=1;i<=n;i++) ans=(LL)ans*phi[i]%MOD;
printf("%d",ans);
return 0;
}