标签：des   style   blog   http   io   ar   color   os   sp   

题目大意：给你一个m*n的方格，方格从（1，1）开始。每个点上有一棵树，然后让你从(0,0)点向下看，问你能看见几颗树。

解题思路：如果你的视线被后面的树和挡住的话以后在这条线上的树你是都看不见的啊。挡住的话就是这个小的方格内对角线的连线过顶点，如图：

建设A为(0,0)如果B能遮挡住C，那么B与C组成的矩形中nx, mx是不互质的。

所以x从1开始枚举n如果m中某一项与x互质那么就会有一个格子显示出来，所以这里需要找的是1-n中与枚举的x互质的个数。这里的话我们可以用到分解质因数，然后用状态压缩进行容斥，求出n中与x互质的个数。

PS：我代码里面写的n，m与这里解释的n，m是反过来的。

Visible Trees

2
1 1
2 3

1
5

#include <algorithm>
#include <iostream>
#include <stdlib.h>
#include <string.h>
#include <iomanip>
#include <stdio.h>
#include <string>
#include <queue>
#include <cmath>
#include <stack>
#include <map>
#include <set>
#define eps 1e-8
#define M 1000100
#define LL long long
//#define LL long long
#define INF 0x3f3f3f
#define PI 3.1415926535898
#define mod 1000000007


const int maxn = 100010;

using namespace std;

struct node
{
    int num[8];
    int ans;
} p[maxn];

int f[maxn];
int k[maxn];
int t;
void prime()
{
    t = 0;
    memset(f, 0, sizeof(f));
    for(int i = 2; i <= maxn-5; i++)
    {
        if(!f[i]) k[t++] = i;
        for(int j = 0; j < t; j++)
        {
            if(i*k[j] > maxn-5) break;
            f[i*k[j]] = true;
            if(i%k[j] == 0) break;
        }
    }
}


LL judge(int n, int x)
{
    LL cnt = 0;
    for(int i = 0; i < (1<<(p[x].ans)); i++)
    {
        int ss = 0;
        int sx = 1;
        for(int j = 0; j < p[x].ans; j++)
        {
            if((1<<j)&i)
            {
                ss++;
                sx *= p[x].num[j];
            }
        }
        if(ss%2) cnt += (n-n/sx);
        else cnt -= (n-n/sx);
    }
    return cnt;
}

int main()
{
    prime();
    for(int i = 2; i <= maxn-10; i++)
    {
        int x = i;
        p[i].ans = 0;
        for(int j = 0; j < t; j++)
        {
            if(x == 1) break;
            if(!f[x])
            {
                p[i].num[p[i].ans++] = x;
                break;
            }
            if(x%k[j]) continue;
            p[i].num[p[i].ans++] = k[j];
            while(x%k[j] == 0) x /= k[j];
        }
    }
    int n, m;
    int T;
    cin >>T;
    while(T--)
    {
        scanf("%d %d",&n, &m);
        LL sum = n;
        for(int i = 2; i <= m; i++) sum += judge(n, i);
        cout<<sum<<endl;
    }
    return 0;
}

HDU 2841 Visible Trees（数论）

标签：des   style   blog   http   io   ar   color   os   sp   

原文地址：http://blog.csdn.net/xu12110501127/article/details/41318289