Hdu 2841 Visible Trees 容斥原理

第一次写容斥原理，虽然有点受挫还是学到了不少

题意：给出一个m*n的矩阵，从（1，1）开始，一个人站在（0,0）位置，问人可以看到矩阵里的几棵树，要求两棵树和人在同一直线上的时候只能看到离自己最近的一棵

思路：首先 我们发现，对于一个点（x,y）只要x与y存在最大公约数g,则可以知道在(x/g,y/g)出有一棵树挡起了点(x,y)，所以(x,y)是看不到的，因此我们要判断一个点是否能看到，就看它是否存在最大公约数不为1的数，若不存在则可以看到，相反则看不到

而满足最大公约数为1的数就是x 和y 互质，所以我们只要找到矩阵内互质的数之和即可,我们可以固定n,然后从[1,m]找有互质的数

问题转化为，对于一个数x，从[1,m]里找与x互质的个数，在求和

两个数字互质，其实就是它们没有公共的质因子，反过来两个数字有公共的质因子则一定不互质，那么我们可以求反面，x与1到m里面多少个数字不互质，然后用m减去即可

然后就会用到容斥原理：先找到有多少个数和x有1个公共的质因子，然后加上；再找到有多少个数与x有2个公共的质因子，然后减去；再找到有多少个数有多少个数与x有3个公共的质因子，然后加上……最后得到的个数，就是有多少个数与x不互质

别人的解释：因为容斥原理一个最基本的准则就是——要计算几个集合并集的大小，我们要先将所有单个集合的大小计算出来，然后减去所有两个集合相交的部分，再加回所有三个集合相交的部分，再减去所有四个集合相交的部分，依此类推，一直计算到所有集合相交的部分。(奇数加，偶数减)

代码：

#include <iostream>
#include <stdio.h>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <vector>
#include <map>
#include <queue>
using namespace std;
#define ll long long
const int M=100010;
vector<int>p[M];
int prime[M];
void init()
{
    for(int i=0; i<M; i++)//1既不是质数，也不是合数。
        p[i].clear();
    memset(prime,0,sizeof(prime));
    for(int i=2; i<M; i++)
    {
        if(!prime[i])//是质数（素数）
        {
            p[i].push_back(i);
            for(int j=2*i; j<M; j+=i)
            {
                prime[j]=1;//不是质因子
                p[j].push_back(i);
            }
        }
    }
}
int cal(int n,int m)
{
    int size=(int)p[n].size();
    int ss=1<<size;
    int sum=0;
    for(int s=1;s<ss;s++)
    {
       int k=0,num=0,phi=1;
       for(int i=0;i<size;i++)
       {
           if(s&(1<<i))
           {
               num++;
               phi*=p[n][i];
           }
       }
       if(num&1)//根据容斥原理
       sum+=m/phi;
       else
       sum-=m/phi;
    }
    return m-sum;
}
int main()
{
    init();
    int t,h,w;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&h,&w);
        if(w>h)swap(w,h);
        long long ans=h;//固定最大的那个值
        for(int i=2;i<=w;i++)
        ans+=cal(i,h);
        printf("%I64d\n",ans);
    }
    return 0;
}