/*
一、容斥原理应用:
1.简单的排列问题由0到9的数字组成排列,要求第一个数大于1,最后一个数小于8,一共有多少种排列?
我们可以来计算它的逆问题,即第一个元素<=1或者最后一个元素>=8的情况。我们设第一个元素<=1时有X组排列,最后一个元素>=8时有Y组排列。
那么通过容斥原理来解决就可以写成:
|X|+|Y|-|X∩Y|然后被总的排列数10!减,就是最终的答案了。
2.序列问题
长度为n的由数字0,1,2组成的序列,要求每个数字至少出现1次,这样的序列有多少种?
3.求指定区间内与n互素的数的个数:给出整数n和r,求区间[1,r]中与n互素的数的个数.
去解决它的逆问题,求不与n互素的数的个数。考虑n的所有素因子pi(i=1…k)在[1,r]中有多少数能被pi整除呢?
它就是:[r/pi]我们可以用2^k的算法求出所有的pi组合,然后计算每种组合的pi乘积,通过容斥原理来对结果进行加减处理。
第一步:
for(i=2;i*i<=n;i++) //对n进行素因子分解
if(n&&n%i==0)
{
prime[m++]=i;
while(n&&n%i==0)
n/=i;
}
if(n>1)
prime[m++]=n;
第二步状压版:
LL prime[maxn];
LL solve(LL num,int m)
{
LL ans=0,mult,i,j,bits;
for(i=1;i<(LL)(1<<m);i++) //用二进制来1,0来表示第几个素因子是否被用到,如m=3,三个因子是2,3,5,则i=3时二进制是011,表示第2、3个因子被用到
{
mult=1,bits=0;
for(j=0;j<m;j++)
if(i&((LL)(1<<j)))//判断第几个因子目前被用到
++bits,mult*=prime[j];
if(bits&1)//容斥原理,奇加偶减
ans+=num/mult;
else
ans-=num/mult;
}
return ans;
}
DFS版
void dfs(int id,int flag,LL num,int ncnt,LL now,int st)//参数暴int注意
{
LL i;
LL temp;
if(id==ncnt+1)
return ;
for(i=st;i<ncnt;++i)
{
temp=now*prime[i];
ans+=flag*(num/temp);
dfs(id+1,(-1)*flag,num,ncnt,temp,i+1);
}
}
hdu 2204
#include
#include
#include
#define ll __int64
using namespace std;
ll prime[20]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61};
ll ans;
ll n;
void dfs(ll st,ll now,ll tt)
{
ll mult;
ll i;
for(i=st;i<19;i++)
{
mult = (ll)(pow(n,1.0/(now*prime[i])));
if(tt%2)
ans+=mult;
else
ans-=mult;
if(mult == 1)
return ;
dfs(i+1,now*prime[i],tt+1);
}
}
int main()
{
while(~scanf("%I64d",&n))
{
ans=0;
dfs(0,1,1);
printf("%I64d\n",ans);
}
return 0;
}
*/
二:欧拉函数
/*
欧拉函数 版本一
*/
#include<cstdio>
#define inf 4000000
int phi[inf];
void getphi()
{
int i,j;
phi[1]=1;
for(i=2;i<=inf;i++)
phi[i]=i;
for(i=2;i<=inf;i++)
if(phi[i]==i)
for(j=i;j<=inf;j+=i)
phi[j]=phi[j]/i*(i-1);
}
int main()
{
getphi();
int a,b,i,sum;
while(scanf("%d%d",&a,&b)!=EOF)
{
__int64 sum=0;
for(i=a;i<=b;i++)
sum+=phi[i];
printf("%I64d\n",sum);
}
return 0;
}
int eular(LL n) //欧拉函数 版本2
{
int i, ans = n;
for(i = 2; i * i <= n; i++)
{
if(n%i == 0)
{
ans -= ans/i;
while(n%i == 0)
n /= i;
}
}
if(n > 1) ans -= ans/n;
return ans;
}
三:筛选素数
/*
筛选素数
*/
#include<cstdio>
#include<cstring>
int b[10050];
int main()
{
int i,j,n;
memset(b,0,sizeof(b));
for(i=2;i<=10000;i++)
if(!b[i])//b[]==0表示为素数
for(j=2*i;j<=10000;j+=i)
b[j]=1;
while(scanf("%d",&n)!=EOF)
{
for(i=n/2;i<n;i++)
if(b[i]==0 && b[n-i]==0)
{
printf("%d %d\n",n-i,i);
break;
}
}
return 0;
}
/*
线性筛选
*/
#define N 100000+10
int pr[N],p[N/10],lp;
void gp(){
for(int i=2;i<N;i++){
if(!pr[i])p[lp++]=pr[i]=i;
for(int j=0;j<lp && i*p[j]<N;j++){
pr[i*p[j]]=p[j];
if(i%p[j]==0)break;
}
}
}
数论研究
最新推荐文章于 2020-04-13 19:29:56 发布