HDU1286 找新朋友(容斥原理+欧拉函数)

Problem Description

新年快到了，“猪头帮协会”准备搞一个聚会，已经知道现有会员N人，把会员从1到N编号，其中会长的号码是N号，凡是和会长是老朋友的，那么该会员的号码肯定和N有大于1的公约数，否则都是新朋友，现在会长想知道究竟有几个新朋友？请你编程序帮会长计算出来。

Input

第一行是测试数据的组数CN（Case number，1

Output

对于每一个N，输出一行新朋友的人数，这样共有CN行输出。

Sample Input

2
25608
24027

Sample Output

7680
16016

思路

求1~n中与n互质的数的个数，欧拉函数和容斥模板题

代码1(欧拉函数)

#include <cstdio>
#include <cstring>
#include <cctype>
#include <stdlib.h>
#include <string>
#include <map>
#include <iostream>
#include <stack>
#include <cmath>
#include <queue>
#include <vector>
#include <algorithm>
using namespace std;
typedef long long ll;
#define inf 1000000
#define mem(a,b) memset(a,b,sizeof(a))
int euler_phi(int n)
{
    int res = n;
    int m = (int)sqrt(n);
    for(int i = 2; i <= m; i++)
        if(n % i == 0)
        {
            res = res / i * (i-1);
            while(n % i == 0) n /= i;
        }
    if(n > 1) res = res / n * (n-1);
    return res;
}
int main()
{
    int t,n;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        printf("%d\n",euler_phi(n));
    }
    return 0;
}

代码2(容斥)

#include <cstdio>
#include <cstring>
#include <cctype>
#include <stdlib.h>
#include <string>
#include <map>
#include <iostream>
#include <stack>
#include <cmath>
#include <queue>
#include <vector>
#include <algorithm>
using namespace std;
typedef long long ll;
#define inf 1000000
#define mem(a,b) memset(a,b,sizeof(a))
int a[50],b[1010];//a保存n的质因子,a[0]表示质因子个数
void div(int n)//分解质因数
{
    int j=0;
    for(int i=2; i*i<=n; i++)
        if(n%i==0)
        {
            while(n%i==0)
                n/=i;
            a[++j]=i;
        }
    if(n>1) a[++j]=n;
    a[0]=j;
}
int get_cnt(int mid)//1--mid之间与n互质的数有多少个
{
    int g=0,sum=mid,t;
    b[++g]=1;
    for(int i=1; i<=a[0]; i++)
    {
        t=g;
        for(int j=1; j<=g; j++)
        {
            b[++t]=b[j]*a[i]*-1;
            sum+=mid/b[t];
        }
        g=t;
    }
    return sum;
}

int main()
{
    int t,n;
    scanf("%d",&t);
    while(t--)
    {
        mem(a,0);mem(b,0);
        scanf("%d",&n);
        div(n);
        printf("%d\n",get_cnt(n));
    }
    return 0;
}