容斥原理:所有情况相加-重复的情况。
hdu1796
推荐指数:****
这是最最基础的模板
题意:给n和m,m个数,在小于n的数中,有多少为m个数的倍数。
思路:m个数的倍数-重复计入的部分,主要是理解容斥原理。
//二进制
#include<cstdio>
#include<iostream>
using namespace std;
int gcd(int a,int b)
{
return b==0?a:gcd(b,a%b);
}
int arr[25];
int main()
{
int n,m;
while(scanf("%d %d",&n,&m)!=EOF)
{
int num=0;
int ans=0;
for(int i=1;i<=m;i++)
{
int ai;
scanf("%d",&ai);
if(ai!=0)//数据好坑,还可能为0
arr[num++]=ai;
}
//num种情况,对应二进制num位(0~num-1)
for(int i=1;i<(1<<num);i++)
{
int lcm=1,sum=0;
for(int j=0;j<num;j++)
{
if(i&(1<<j))
{
sum++;//位上为一的个数
lcm=lcm/gcd(lcm,arr[j])*arr[j];//最小公倍数
}
}
if(sum%2==0)ans-=(n-1)/lcm;//偶减
else ans+=(n-1)/lcm;//奇加
}
printf("%d\n",ans);
}
return 0;
}
//dfs
#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
int arr[25];
int n,m,num,ans;
int gcd(int a,int b)
{
return b==0?a:gcd(b,a%b);
}
void dfs(int a,int lcm,int id)
{
lcm=arr[a]/gcd(lcm,arr[a])*lcm;
if(id%2==0)
ans-=(n-1)/lcm;
else
ans+=(n-1)/lcm;
//深搜,eg:A->AB->ABC,B->BC,C
for(int i=a+1;i<num;i++)
dfs(i,lcm,id+1);
}
int main()
{
while(scanf("%d %d",&n,&m)!=EOF)
{
num=0,ans=0;
for(int i=1;i<=m;i++)
{
int ai;//注意0这种情况
cin>>ai;
if(ai!=0)
arr[num++]=ai;
}
for(int i=0;i<num;i++)
dfs(i,arr[i],1);
cout<<ans<<endl;
}
return 0;
}
hdu 4135
推荐指数***
题意:给定区间[a,b],n,求在[a,b]中有多少个和n互质的数;
思路:[a,b]总数-与n不互质的数=ans;求出n的质因子,利用容斥原理,求出质因子在[a,b]的倍数个数-重复部分;
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
ll num;
ll arr[20];
void init(ll n)//分解质因子+二进制化
{
num=0;
for(ll i=2;i*i<=n;i++)
if(n%i==0)
{
while(n%i==0)n/=i;
arr[num++]=i;
}
if(n>1)arr[num++]=n;
}
ll solve(ll n)//
{
ll ans=0;
for(ll i=1;i<(1ll<<num);i++)
{
ll l=1,s=0;
for(ll j=0;j<num;j++)
if(i>>j&1)
{
s++;//该为为1
l*=arr[j];//已是质数,不用gcd
}
if(s&1)ans+=n/l;
else ans-=n/l;
}
return ans;
}
int main()
{
int t;
ll a,b,n;
scanf("%d",&t);
int cas=0;
while(t--)
{
scanf("%lld%lld%lld",&a,&b,&n);
init(n);
printf("Case #%d: %lld\n",++cas,b-a+1-solve(b)+solve(a-1));
}
return 0;
}
hdu6286
推荐指数***
题意:给两个区间[a,b] [c,d],求有多少个(x,y)存在(a<=x<=b,c<=y<=d),使xy=k*2018(k是整数);
思路:2018的因子只有1 2 1009 2018;只需判断在[a,b]区间其因子个数乘对应另一因子在[c,d]的个数,再减去重复计算的。
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cmath>
using namespace std;
typedef long long ll;
//在当前[l,r]区间x因子的个数
ll fac(int l,int r,int x)
{
return r/x-(l-1)/x;
}
int main()
{
int l1,r1,l2,r2;
while(~scanf("%d%d%d%d",&l1,&r1,&l2,&r2))
{
ll ans=0;
ans+=(fac(l1,r1,2018)*(r2-l2+1));
ans+=(fac(l2,r2,2018)*(r1-l1+1));
ans-=(fac(l1,r1,2018)*fac(l2,r2,2018));
ans+=((fac(l1,r1,1009)-fac(l1,r1,2018))*(fac(l2,r2,2)-fac(l2,r2,2018)));
ans+=((fac(l2,r2,1009)-fac(l2,r2,2018))*(fac(l1,r1,2)-fac(l1,r1,2018)));
printf("%lld\n",ans);
}
return 0;
}
hdu2204
推荐指数:*****
题意:给你一个数n,求[1,n]之间有多少数可以表示为mk,eg:8=23;64=82=26;可见求出质数的次方包含了合数的情况。
下班了下班了 明天写
上线了理解借鉴
#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
typedef long long ll;
int prime[30];
int main()
{
//求质数:10^18约等于2^64,求100以内即可
int num=0;
prime[0]=2;
for(int i=3;i<100;i++)
{
int flag=1;
for(int j=2;j<i;j++)
{
if(i%j==0)
{
flag=0;
break;
}
}
if(flag)prime[++num]=i;
}
//2,3,5次方相当于容斥定理的三个集合;
ll n;
while(cin>>n)
{
ll ans=0;
ll temp;
for(int i=0;;i++)
{
temp=pow(n,1.0/prime[i]);
if(temp<2)break;
ans+=temp-1;
//eg:10,10^(1/2)=3;所以1,2,3的平方都小于10,减去1;(1单独考虑)
for(int j=i+1;;j++)
{
temp=pow(n,1.0/prime[i]/prime[j]);
if(temp<2)break;
ans-=temp-1;
for(int k=j+1;;k++)
{
temp=pow(n,1.0/prime[i]/prime[j]/prime[k]);
if(temp<2)break;
ans+=temp-1;
}
}
}//终于懂了。第一轮:开2次、2*3次、2*3*5次;第二轮:开3次、3*5次、3*5*7次;所以素数存到>64;
cout<<ans+1<<'\n';
}
return 0;
}