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
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 ;
}