数论研究

/*

一、容斥原理应用：

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;  
        }  
    }  
}