UVa 11064 - Number Theory

标签：blog   io   os   ar   for   sp   数据   on   2014   

题目：求给顶一个数n，的所有的1 ≤ m ≤ n的m，使得gcd（m，n）≠ 1 且 gcd（m，n）≠ m。

分析：数论，素数筛法，欧拉函数。 

            设pi为n的第i个素数因，k1为第i个素数因子的个数，则有：

            1 ≤ m ≤ n，gcd（m，n）= 1 的m的个数为欧拉函数；

            欧拉函数：φ（n）= n *（1 - 1/p1）*（1 - 1/p2）*（1 - 1/p3）*…*（1 - 1/pt）；

            1 ≤ m ≤ n，gcd（m，n）= m 的m的个数为n的所有因数的个数；

            因数个数：f（n）= （k1+1）*（k2+1）*...*（kt+1）；

            这里利用筛法打表计算出150000内的素数，因为数据范围是20000000000内的，

            所以，不能被前150000内的素数整除的数，也一定是素数，并且每个数n中最多有一个；

            计算输出 n -  φ（n）- f（n）+ 1 即可{ gcd（1，n）计算了两次 }。

说明：UVa10299。

#include <iostream>  
#include <cstdlib>  
#include <cmath>
  
using namespace std;  
  
int fac[32],num[32];  
int prim[150000];  
int used[150000];  
  
int main()  
{  
    for (int i = 0 ; i < 150000 ; ++ i)  
        used[i] = 0;  
    int save = 0;  
    for (int i = 2 ; i < 150000 ; ++ i)  
        if (!used[i]) {  
            prim[save ++] = i;  
            for (int j = 2*i ; j < 150000 ; j += i)  
                used[j] = 1;  
        }  
      
    int n;  
    while (cin >> n && n) {  
        int count = 0,base = 0,m = n;  
        while (n > 1 && base < save) {  
            if (n%prim[base] == 0) {  
                fac[count] = prim[base];
                num[count] = 0;
                while (n%prim[base] == 0) {
                    n /= prim[base];  
                    num[count] ++;
				}
				count ++;
            }  
            base ++;
        }
        if (n > 1) {
			fac[count] = n;
			num[count] = 1;
			count ++;
		}
          
        long long ans = m;
        for (int i = 0 ; i < count ; ++ i)  
            ans = ans/fac[i]*(fac[i]-1);
        
        //gcd(m,n) = m
        int s = 1;
        for (int i = 0 ; i < count ; ++ i)
        	s = s*(num[i]+1);
        
        cout << m-ans-s+1 << endl;  
    }  
    return 0;  
}

UVa 11064 - Number Theory

标签：blog   io   os   ar   for   sp   数据   on   2014   

原文地址：http://blog.csdn.net/mobius_strip/article/details/40535315