Bzoj3309 DZY Loves Math

Time Limit: 20 Sec  Memory Limit: 512 MB
Submit: 992  Solved: 589

Description

对于正整数n，定义f(n)为n所含质因子的最大幂指数。例如f(1960)=f(2^3 * 5^1 * 7^2)=3, f(10007)=1, f(1)=0。
给定正整数a,b，求sigma(sigma(f(gcd(i,j)))) (i=1..a, j=1..b)。

Input

第一行一个数T，表示询问数。
接下来T行，每行两个数a,b，表示一个询问。

Output

对于每一个询问，输出一行一个非负整数作为回答。

Sample Input

4
7558588 9653114
6514903 4451211
7425644 1189442
6335198 4957

Sample Output

35793453939901
14225956593420
4332838845846
15400094813

HINT

【数据规模】

T<=10000

1<=a,b<=10^7

Source

莫比乌斯反演 脑洞题

题解传送门http://blog.csdn.net/sdfzyhx/article/details/72854335

 1 /*by SilverN*/
 2 #include<iostream>
 3 #include<algorithm>
 4 #include<cstdio>
 5 #include<cmath>
 6 #include<cstring>
 7 #define LL long long
 8 using namespace std;
 9 const int mxn=10000010;
10 int read(){
11     int x=0,f=1;char ch=getchar();
12     while(ch<‘0‘ || ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();}
13     while(ch>=‘0‘ && ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();}
14     return x*f;
15 }
16 int pri[mxn],cnt;
17 int mu[mxn],last[mxn],w[mxn],g[mxn];
18 bool vis[mxn];
19 LL smm[mxn];
20 void init(){
21     mu[1]=1;
22     for(int i=2;i<mxn;i++){
23         if(!vis[i]){
24             pri[++cnt]=i;last[i]=1;w[i]=1;
25             mu[i]=-1;g[i]=1;
26         }
27         for(int j=1;j<=cnt && pri[j]*i<mxn;j++){
28             int tmp=pri[j]*i;
29             vis[tmp]=1;
30             if(i%pri[j]==0){
31                 last[tmp]=last[i];
32                 w[tmp]=w[i]+1;
33                 mu[tmp]=0;
34                 if(last[tmp]==1) g[tmp]=1;
35                 else    if(w[last[tmp]]==w[tmp])g[tmp]=-g[last[tmp]];
36                         else g[tmp]=0;
37                 break;
38             }
39             last[tmp]=i;
40             mu[tmp]=-mu[i];
41             w[tmp]=1;
42             g[tmp]=(w[i]==1)?(-g[i]):0;
43         }
44     }
45     for(int i=1;i<mxn;i++)smm[i]=smm[i-1]+g[i];
46     return;
47 }
48 int a,b;
49 void solve(){
50     if(a>b)swap(a,b);
51     LL ans=0;
52     for(int i=1,pos;i<=a;i=pos+1){
53         pos=min(a/(a/i),b/(b/i));
54         ans+=(smm[pos]-smm[i-1])*(LL)(a/i)*(b/i);
55     }
56     printf("%lld\n",ans);
57     return;
58 }
59 int main(){
60 //    freopen("in.txt","r",stdin);
61     init();
62     int T=read();
63     while(T--){
64         a=read();b=read();
65         solve();
66     }
67     return 0;
68 }