这份代码
#include<bits/stdc++.h>
#define ll long long
#define reg register
#define db double
#define il inline
using namespace std;
int t,n;
ll ans[1000005];
ll phi[1000005],prime[1000005],tot;
void init(){
phi[1]=1;
for(reg int i=2;i<=1000000;i++){
if(phi[i]==0){
phi[i]=i-1;
prime[++tot]=i;
}
for(int j=1;j<=tot&&i*prime[j]<=1000000;j++){
// if(i==16129)cout<<j*prime[]
if(i%prime[j]==0){
// if(prime[j]*i>1000000)cout<<prime[j]*i<<endl;
phi[prime[j]*i]=phi[i]*prime[j];
break;
}else{
phi[prime[j]*i]=phi[i]*phi[prime[j]];
}
}
}
for(int i=1;i<=1000000;i++){
for(int j=i;j<=1000000;j+=i){
if(j==2){
// cout<<i<<' '<<phi[i]<<endl;
}
if(i==1)ans[j]++;
else ans[j]+=phi[i]*i/2;
}
}
}
int main(){
init();
cin>>t;
while(t--){
cin>>n;
cout<<ans[n]*n<<endl;
}
return 0;
}
将24行改为
phi[prime[j]*i]=phi[i]*prime[j]-phi[i];
就爆炸了怎么回事呢