Bzoj 2749: [HAOI2012]外星人 欧拉函数,数论,线性筛

2749: [HAOI2012]外星人

Time Limit: 3 Sec  Memory Limit: 128 MB
Submit: 568  Solved: 302
[Submit][Status][Discuss]

Description

Input

Output

输出test行，每行一个整数，表示答案。

Sample Input

1
2
2 2
3 1

Sample Output

3

HINT

Test<=50 Pi<=10^5,1<=Q1<=10^9

Source

很好的一道题。

就是要求能phi出多少个2。。。

奇数时要加一(就是刚开始的时候变成偶数时要用一次。)

 #include<bits/stdc++.h>
 using namespace std;
 #define MAXN 100000
 #define LL long long
 int prime[],phi[MAXN+],tot,f[MAXN+];
 bool vis[MAXN+];
 int read()
 {
     int s=,fh=;char ch=getchar();
     while(ch<''||ch>''){if(ch=='-')fh=-;ch=getchar();}
     while(ch>=''&&ch<=''){s=s*+(ch-'');ch=getchar();}
     return s*fh;
 }
 void getEular()
 {
     int i,j;
     phi[]=;tot=;
     for(i=;i<=MAXN;i++)
     {
         if(vis[i]==false)
         {
             prime[++tot]=i;
             phi[i]=i-;
         }
         for(j=;j<=tot&&prime[j]*i<=MAXN;j++)
         {
             vis[prime[j]*i]=true;
             if(i%prime[j]==)
             {
                 phi[prime[j]*i]=prime[j]*phi[i];
                 break;
             }
             phi[prime[j]*i]=phi[prime[j]]*phi[i];
         }
     }
 }
 int main()
 {
     freopen("alien.in","r",stdin);
     freopen("alien.out","w",stdout);
     int T,i,m,p,q,n;
     bool flag;
     LL ans;
     T=read();
     getEular();
     memset(f,,sizeof(f));//f[i]代表i需要phi多少次才能化为1.(也就是等于能phi多少个2.)
     f[]=-;
     for(i=;i<=MAXN;i++)f[i]=f[phi[i]]+;
     f[]++;f[]++;
     while(T--)
     {
         m=read();ans=;
         flag=false;//判断奇数时要加一.
         for(i=;i<=m;i++)
         {
             p=read();q=read();
             if(p==)flag=true;
             ans+=(LL)f[p]*q;
         }
         if(flag==false)ans++;
         printf("%lld\n",ans);
     }
     fclose(stdin);
     fclose(stdout);
     return ;
 }