Poj.2480.Longge's problem

J - A/B

Time Limit: 1000 MS Memory Limit: 32768 KB

64-bit integer IO format: %I64d , %I64u Java class name: Main

[Submit] [Status]

Description

要求(A/B)%9973，但由于A很大，我们只给出n(n=A%9973)(我们给定的A必能被B整除，且gcd(B,9973) = 1)。

Input

数据的第一行是一个T，表示有T组数据。
每组数据有两个数n(0 <= n < 9973)和B(1 <= B <= 10^9)。

Output

对应每组数据输出(A/B)%9973。

Sample Input

2
1000 53
87 123456789

Sample Output

7922
6060

考点：欧拉函数

#include<cmath>
#include<cstdio>

long long phi(long long n){
    long long ans=n;
    for(long long i=2;i<=sqrt(n);i++){
        if(n%i==0){
            ans=ans/i*(i-1);
            while(n%i==0)n/=i;
        }
    }
    if(n>1)ans=ans/n*(n-1);
    return ans;
}

int main()
{
    long long n;
    while(~scanf("%lld",&n)){
        long long ans=0,i;
        for(i=1;i<sqrt(n);i++){
    if(n%i==0)
        ans+=i*phi(n/i)+n/i*phi(i);}
    if(i*i==n)ans+=i*phi(i);
        printf("%lld\n",ans);
    }
}