Hdu2841容斥原理

人在(0,0)点,问n*m的矩阵上的点有多少可以与人直接可见,其实就是矩阵上点与(0,0)点的形成的斜率种类数。

#pragma comment(linker,"/STACK:102400000,102400000")
#define _CRT_SECURE_NO_WARNINGS
#include<stdio.h>
#include<iostream>
#include<string>
#include<queue>
#include<stack>
#include<list>
#include<stdlib.h>
#include<algorithm>
#include<vector>
#include<map>
#include<cstring>
#include<set>
using namespace std;
typedef long long LL;

LL cnt, ans, m;
LL q[20];
void gao(LL x, LL pre, LL flag, LL key)
{
    if (x == cnt) {
        if (flag == 0) ans -= key / pre;
        else ans += key / pre;
        return;
    }
    gao(x + 1, pre, flag, key);
    gao(x + 1, pre*q[x], flag ^ 1, key);
}

LL ask(LL n) {
    cnt = 0;
    LL t = n;
    ans = 0;
    for (LL i = 2; i*i <= t; i++) {
        if (t%i)  continue;
        while (t%i == 0) t /= i;
        q[cnt++] = i;
    }
    if (t > 1) q[cnt++] = t;
    gao(0, 1, 1, m);
    return ans;

}

int main()
{
    LL T, n;
    cin >> T;
    while (T--) {
        cin >> n >> m;
        if (m < n) swap(n, m);
        LL sum = 0;
        for (LL i = 1; i <= n; i++) {
            sum += ask(i);
        }
        cout << sum << endl;
    }
    return 0;
}
时间: 2024-10-09 21:12:56

Hdu2841容斥原理的相关文章

HDU2841 Visible Trees (容斥原理)

主题链接:http://acm.hdu.edu.cn/showproblem.php?pid=2841 题意: 一个人在(0,0)点,然后前面有一个m*n的格子 ,每一个格子的节点上有一棵树.问这个人站在原地能看到多少棵树 假设两棵树在一条直线上那么仅仅能看到最前面的一棵树. 分析: 假设一个数的坐标为(a,b).那么坐标为(a*k,b*k)的都不能看见.假设a,b有公因子c那么我们肯定仅仅能看到(a/c,b/c). 因此我们得出结论,能看到的树的横纵坐标一定互质. 那么我们就能够把问题转化为,

HDU2841 Visible Trees(容斥原理)

题目..大概就是有个m*n个点的矩形从(1,1)到(m,n),问从(0,0)出发直线看过去最多能看到几个点. 如果(0,0)->(x,y)和(0,0)->(x',y')两个向量平行,那后面的那个点就看不到了. 因此给出一个点(x,y),判断它能否被看到,就是是否能找到一个大于1的k,使k|x且k|y. 这样,问题就能转变为有几个点的x.y找不到公约数,即有几对x.y,满足x和y互质. 可以通过枚举x,看有几个y与其互质累加.这样问题就又变成,区间有几个数与某个数互质,经典的容斥问题HDU413

容斥原理——hdu2841

记得要开ll /* 莫比乌斯反演模板题,也可以直接算phi来做 容斥的解法 求x[1..m],在[1,n]中和其互质的数的个数即可 那么就是n-和x不互质的数个数即可 */ #include<bits/stdc++.h> #include<vector> using namespace std; #define maxn 100005 #define ll long long vector<int>fac[maxn]; ll n,ans; void init(){//质

UVA 11014 - Make a Crystal(容斥原理)

UVA 11014 - Make a Crystal 题目链接 题意:给定一个NxNxN的正方体,求出最多能选几个整数点.使得随意两点PQ不会使PQO共线. 思路:利用容斥原理,设f(k)为点(x, y, z)三点都为k的倍数的点的个数(要扣掉一个原点O).那么全部点就是f(1),之后要去除掉共线的,就是扣掉f(2), f(3), f(5)..f(n).n为素数.由于这些素数中包括了合数的情况,而且这些点必定与f(1)除去这些点以外的点共线,所以扣掉.可是扣掉后会扣掉一些反复的.比方f(6)在f

BZOJ3198 SDOI2013 spring HASH+容斥原理

题意:给定6个长度为n的数列,求有多少个数对(i,j)((i,j)≡(j,i))使得i和j位置恰好有K个数相同,其中0≤K≤6 题解: 设fi=至少有K个数相同的位置对的数量,用2^6枚举每一种可能,根据容斥原理,答案就是\[\sum\limits_{i = K}^N {{f_i}C_i^K{{\left( { - 1} \right)}^{i - K}}} \] 至于多乘一个组合数,是因为当前枚举到有x个数相同,一个位置对有i个相同的数,那么累计的时候就会算成$C_x^i$,因此实际上这个位置

hdu 4790 Just Random 容斥原理+数学

Just Random Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 1723    Accepted Submission(s): 483 Problem Description Coach Pang and Uncle Yang both love numbers. Every morning they play a game w

[BZOJ 1042] 硬币购物 容斥原理

题意 有四种货币, 它们的价值分别是 c[0], c[1], c[2], c[3] . n 次询问, 每次给定 d[0], d[1], d[2], d[3], s, 问凑出 s , 且第 i 种货币不超过 c[i] 个的方案数. c[i], d[i], s <= 100000 , n <= 1000 . 分析 设第 i 种货币取了 x[i] 个. 问题转化为求不定方程 c[0]x[0] + c[1]x[1] + c[2]x[2] + c[3]x[3] = s 的非负整数解的个数. 且满足 4

HDU - 4135 Co-prime(容斥原理)

Question 参考 题意找出[a,b]中与n互质的数的个数分析通常我们求1-n中与n互质的数的个数都是用欧拉函数.但如果n比较大或者是求1-m中与n互质的数的个数等等问题,要想时间效率高的话还是用容斥原理.先对n分解质因数,分别记录每个质因数, 那么所求区间内与某个质因数不互质的个数就是 m/r(i),假设r(i)是r的某个质因子 假设只有三个质因子, 总的不互质的个数应该为p1+p2+p3-p1*p2-p1*p3-p2*p3+p1*p2*p3. pi代表m/r(i),即与某个质因子不互质的

容斥原理

对容斥原理的描述 容斥原理是一种重要的组合数学方法,可以让你求解任意大小的集合,或者计算复合事件的概率. 描述 容斥原理可以描述如下: 要计算几个集合并集的大小,我们要先将所有单个集合的大小计算出来,然后减去所有两个集合相交的部分,再加回所有三个集合相交的部分,再减去所有四个集合相交的部分,依此类推,一直计算到所有集合相交的部分. 关于集合的原理公式 上述描述的公式形式可以表示如下:                  它可以写得更简洁一些,我们将B作为所有Ai的集合,那么容斥原理就变成了: 这个