hdu 2841（容斥原理+欧拉函数）

题意：
给出一个矩阵问在点(0,0)的人最多能看得到多少点。不难发现，gcd(x,y)=1的时候就是满足的点，其实可以容斥直接秒了，但是我写的时候觉得1e7有点虚，就先用欧拉函数算了正方形的地方，剩下的再用容斥写，确实跑的比裸的容斥要快，就是求在1~n内和 2 ~m互质的有多少个，正难则反，质因数分解后用容斥算出不互质的，减去即可，其实欧拉函数、莫比乌斯等函数都是用容斥证明出来的
所以答案就是ans=1+∑phi[i] i=(2~min(n,m))+(min(n,m)+1 ~max(m,n)）的地方容斥

#include<iostream>
#include<cstdio>
#include<cstring>

using namespace std;
typedef long long ll;
const int maxn=1e5+10;
bool v[maxn];
int phi[maxn],cnt=0,prime[maxn],n,m,num,p[10],sum,i;
ll ans;

void init()
{
    v[1]=v[0]=1;phi[1]=1;
    for(int i=2;i<=maxn-9;++i)
    {
        if(!v[i])phi[i]=i-1,prime[++cnt]=i;
        for(int j=1;j<=cnt&&prime[j]*i<=maxn-9;++j)
        {
            v[prime[j]*i]=1;
            if(i%prime[j]==0)
            {
                phi[prime[j]*i]=phi[i]*prime[j];
                break; 
            } 
            else phi[prime[j]*i]=phi[prime[j]]*phi[i];
        } 
    } 
}
void dfs(int pos,int lcm,int id)
{
    if(id==0)
    {
       if(i&1)num+=(m/lcm);
       else num-=(m/lcm); 
       return;
    }
    if(pos>sum)return;
    if(lcm*p[pos]<=m)
        dfs(pos+1,lcm*p[pos],id-1);
    dfs(pos+1,lcm,id);
    return;
}
int getans(int k)
{
    for(int i=2;i*i<=k;++i)
    {
       if(k%i==0){ 
          p[++sum]=i;   
          while(k%i==0)k/=i; 
       }
    } 
    if(k>1)p[++sum]=k;
    for(i=1;i<=sum;++i)
        dfs(1,1,i);
    return m-num;
}
int main()
{
    init(); 
    int t;
    scanf("%d",&t);
    while(t--)
    { 
        scanf("%d%d",&m,&n);
        ans=1;
        if(m>n)swap(n,m);
        for(int i=2;i<=m;++i)
            ans+=2ll*phi[i];
        for(int i=m+1;i<=n;++i)
        { 
            num=sum=0;
            ans+=getans(i);
        }
        printf("%lld\n",ans);
    }
    return 0; 
}