容斥原理——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(){//质因子分解,埃氏筛
    for(int i=2;i<=100000;i++){
        if(fac[i].size())continue;
        fac[i].push_back(i);
        for(int j=2;j*i<=100000;j++)
            fac[i*j].push_back(i);
    }
}
//枚举的m,在m的质因子里的位置pos,当前的因数num,步长
void dfs(int m,ll pos,ll num,int cnt){
    //cout<<fac[m].size()<<‘\n‘;
    if(num>n)return;
    if(pos==fac[m].size())return;
    if(cnt%2)ans+=n/num;
    else ans-=n/num;
    for(int i=pos+1;i<fac[m].size();i++)
        dfs(m,i,num*fac[m][i],cnt+1);
}
int main(){
    init();
    int m,t;
    cin>>t;while(t--){
        cin>>m>>n;
        ans=n;//x==1时
        for(int i=2;i<=m;i++)dfs(i,-1,1,1);
        cout<<ans<<endl;
    }
}

原文地址:https://www.cnblogs.com/zsben991126/p/10861997.html

时间: 2024-10-09 21:13:02

容斥原理——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容斥原理

人在(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&

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的集合,那么容斥原理就变成了: 这个