社区讨论
求帮助卡常 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 条回复,欢迎继续交流。
正在加载回复...