记录编号 |
286029 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[BZOJ 2820] YY的GCD |
最终得分 |
100 |
用户昵称 |
TenderRun |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
18.072 s |
提交时间 |
2016-07-26 09:05:11 |
内存使用 |
124.29 MiB |
显示代码纯文本
- #include <iostream>
- #include <cstring>
- #include <cstdio>
- using namespace std;
- const int maxn=10000010;
- int prime[maxn],cnt;
- int mu[maxn],sum[maxn];
- bool check[maxn];
-
- void Prepare(){
- mu[1]=1;
- for(int i=2;i<=10000000;i++){
- if(!check[i]){
- prime[++cnt]=i;
- mu[i]=-1;
- }
- for(int j=1;j<=cnt;j++){
- if(prime[j]*i>10000000)break;
- check[prime[j]*i]=true;
- if(i%prime[j]==0){
- mu[prime[j]*i]=0;
- break;
- }
- mu[prime[j]*i]=mu[i]*-1;
- }
- }
- for(int i=1;i<=10000000;i++)
- for(int j=1;j<=cnt&&prime[j]*i<=10000000;j++)
- sum[i*prime[j]]+=mu[i];
- for(int i=1;i<=10000000;i++)sum[i]+=sum[i-1];
- }
-
- int T;
- int a,b;
- long long C(int n,int m){
- if(n>m)swap(n,m);
- long long ret=0;int p;
- for(int i=1;i<=n;i=p+1){
- p=min(n/(n/i),m/(m/i));
- ret+=1ll*(sum[p]-sum[i-1])*(n/i)*(m/i);
- }
- return ret;
- }
-
- int main(){
- freopen("YYnoGCD.in","r",stdin);
- freopen("YYnoGCD.out","w",stdout);
- Prepare();
- scanf("%d",&T);
- while(T--){
- scanf("%d%d",&a,&b);
- printf("%lld\n",C(a,b));
- }
- return 0;
- }