zoukankan      html  css  js  c++  java
  • 数论模版-欧拉函数、莫比乌斯函数和素数

    /*Author:WNJXYK*/
    #include<cstdio>
    #include<iostream>
    #include<cstring>
    #include<string>
    #include<algorithm>
    #include<queue>
    #include<set>
    #include<map>
    using namespace std;
    
    #define LL long long
    #define Inf 2147483647
    #define InfL 10000000000LL
    
    inline void swap(int &x,int &y){int tmp=x;x=y;y=tmp;}
    inline void swap(LL &x,LL &y){LL tmp=x;x=y;y=tmp;}
    inline int remin(int a,int b){if (a<b) return a;return b;}
    inline int remax(int a,int b){if (a>b) return a;return b;}
    inline LL remin(LL a,LL b){if (a<b) return a;return b;}
    inline LL remax(LL a,LL b){if (a>b) return a;return b;}
    
    const int Maxn=10000000;
    
    
    LL prime[Maxn/2+10];
    bool valid[Maxn+10];
    int primes;
    inline void getPrime(){
    	memset(valid,true,sizeof(valid));
    	for (int i=2;i<=Maxn;i++){
    		if (valid[i])prime[++primes]=i;
    		for (int j=1;j<=primes && prime[j]*i<=Maxn;j++){
    			valid[prime[j]*i]=false;
    			if (i%prime[j]==0) break;
    		}
    	}
    }
    
    
    LL miu[Maxn+10];
    inline void getMiu(){
    	for (int i=1;i<=Maxn;i++){
    		int target=(i==1?1:0);
    		int delta=target-miu[i];
    		miu[i]=delta;
    		for (int j=i+i;j<=Maxn;j+=i) miu[j]+=delta;
    	}
    }
    
    LL phi[Maxn+10];
    LL minDiv[Maxn+10];
    LL sum[Maxn+10];
    inline void getPhi(){
    	for (int i=1;i<=Maxn;i++) minDiv[i]=i;
    	for (int i=2;i<=i*i;i++)
    		if (minDiv[i]=i)
    			for (int j=i*i;j<=Maxn;j+=i)
    				minDiv[j]=i;
    	phi[1]=1;
    	for (int i=2;i<=Maxn;i++){
    		phi[i]=phi[i/minDiv[i]];
    		if ((i/minDiv[i])%minDiv[i]==0){
    			phi[i]*=minDiv[i];
    		}else{
    			phi[i]*=minDiv[i]-1;
    		}
    	}
    }
    
    /*
    int primes,prime[Maxn],mu[Maxn],g[Maxn],sum[Maxn];
    void getPrime_Miu_G_S(){
        static bool vis[Maxn];
        memset(vis,0,sizeof(vis));
        mu[1]=1;primes=0;
        for(int i=2;i<Maxn;i++){
            if(!vis[i])
                prime[primes++]=i,mu[i]=-1,g[i]=1;
            for(int j=0;j<primes&&i*prime[j]<Maxn;j++){
                vis[i*prime[j]]=1;
                if(i%prime[j])
                    mu[i*prime[j]]=-mu[i],g[i*prime[j]]=mu[i]-g[i];
                else{ 
    <span style="white-space:pre">				</span>mu[i*prime[j]]=0;g[i*prime[j]]=mu[i];break;}
            }
        }
        sum[0]=0;
        for(int i=1;i<Maxn;i++)
            sum[i]=sum[i-1]+g[i];
    }*/
    
    int main(){
    
    	return 0;
    }
    

  • 相关阅读:
    JavaScript--事件绑定及深入(26)
    JavaScript--事件对象(25)
    JavaScript--事件入门(24)
    JavaScript--动态加载脚本和样式(23)
    JavaScript--DOM元素尺寸和位置(22)
    JavaScript--DOM操作表格及样式(21)
    JavaScript--DOM进阶(20)
    jQuery-瀑布流的取巧做法分享
    Less开发指南(二)- 基本语法
    less开发指南(一)- 小牛试刀
  • 原文地址:https://www.cnblogs.com/WNJXYK/p/4063923.html
Copyright © 2011-2022 走看看