社区讨论

求帮助卡常 90 pts TLE on 9

P2257YY的GCD参与者 8已保存回复 20

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
20 条
当前快照
1 份
快照标识符
@mhjntewv
此快照首次捕获于
2025/11/04 05:36
4 个月前
此快照最后确认于
2025/11/04 06:33
4 个月前
查看原帖
CPP
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=1e7+5;
int pri[N],u[N],S[N],pum=0;
bool prime[N];
void init(int n){
	u[1]=1;
	for(int i=1;i<=n;i++) u[i]=1;
	for(int i=2;i<=n;i++){
		if(prime[i]==0)
			u[i]=-1,pri[++pum]=i;
		for(int j=1;j<=pum&&i*pri[j]<=n;j++){
			prime[i*pri[j]]=1;
			u[i*pri[j]]=u[pri[j]]*u[i];
			if(i%pri[j]==0)
				u[i*pri[j]]=0;
		}
	}
	for(int d=1;d<=n;d++)
		if(u[d]!=0)
			for(int i=1;i<=pum&&d*pri[i]<=n;i++)
				S[d*pri[i]]+=u[d];
	for(int i=1;i<=n;i++)
		pri[i]=pri[i-1]+S[i];
}
long long work(int n,int m){
	long long ans=0; 
	for(int l=1,r;l<=min(n,m);l=r+1){
		r=min(n/(n/l),m/(m/l));
		ans+=(pri[r]-pri[l-1])*(n/r)*(m/r);
	}
	return ans;
}
signed main() {
	ios::sync_with_stdio(false);
	cin.tie(0);cout.tie(0);
	init(1e7);
	int t;
	cin>>t;
	while(t--){
		int n,m;
		cin>>n>>m;
		cout<<work(n,m)<<'\n';
	} 
  	return 0;
}


回复

20 条回复,欢迎继续交流。

正在加载回复...