求助
  • 板块P1891 疯狂 LCM
  • 楼主suzhikz
  • 当前回复0
  • 已保存回复0
  • 发布时间2024/9/19 20:25
  • 上次更新2024/9/19 20:36:13
查看原帖
求助
800499
suzhikz楼主2024/9/19 20:25

这份代码

#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];

就爆炸了怎么回事呢

2024/9/19 20:25
加载中...