11.2 NOIP模拟赛 (afternoon)

欧拉函数(phi)

【题目描述】

已知N，求phi(N)

【输出说明】

正整数N

【输入说明】

输出phi(N)

【样例输入】

8

【样例输出】

4

【数据范围】

对于20%的数据,N<=10^5
对于60%的数据,N<=10^9
对于100%的数据,N<=10^18

思路:

裸欧拉函数

代码:

#include<iostream>
#include<cstdio>
#define ll long long
using namespace std;
ll n,ans;
int main()
{
    freopen("phi.in","r",stdin);
    freopen("phi.out","w",stdout);
    cin>>n;ans=n;
    for(ll i=2;i*i<=n;i++){
        if(n%i==0){
            while(n%i==0)n/=i;
            ans=ans/i*(i-1);
        }
    }
    if(n>1)ans=ans/n*(n-1);
    cout<<ans<<endl;
    fclose(stdin);fclose(stdout);
    return 0;
}

思路:

公式变形 F[n]=p1^(a1-1)*(p1-1) + p2^(a2-1)*(p2-1) + ....
Dfs 每个 p a 
这里只筛了10^7 的素数 如果有质因子 比这个大
可以特判掉 判断比较大的数是不是素数 这里用的 Miller Rabin 或者 Fermat
这不是重点 说不定sqrt暴力判常熟写的小一点也能过 没试过233 
然后 利用欧拉函数的性质 F[x]=x-1 x是素数 剪剪枝就好了 大约两秒左右的时间 

代码;

#include<cstdio>
#include<ctime>
#include<iostream>
#include<cstdlib>
#include<algorithm>
#define ll long long
#define maxn 100010
#define MAXN 10000000
using namespace std;
ll n,k,prime[MAXN/10],num,ans[maxn];
bool f[MAXN+10];
inline void prim(){
    for(int i=2;i<=MAXN;i++){
        if(f[i]==0)prime[++num]=i;
        for(int j=1;j<=num&&i*prime[j]<=MAXN;j++){
            f[i*prime[j]]=1;
            if(i%prime[j]==0)break;
        }
    }
}
inline ll mul(ll x,ll y,ll z){
    ll ans=0;
    while(y){
        if(y&1)
            ans+=x,ans%=z,y--;
        x<<=1;x%=z;y>>=1;
    }
    return ans;
}
inline ll Qc(ll x,ll y,ll z){
    ll ans=1;
    while(y){
        if(y&1)ans=mul(ans,x,z);
        y>>=1;x=mul(x,x,z);
    }
    return ans;
}
inline bool Is_prime(ll n)
{  
    if(n<2)return 0;
    if(n==2)return 1;
    if(n%2==0)return 0;
    ll m=n-1,j=0;
    while(m%2==0) {
          j++;
          m >>=1;
    }
    for(int i=1;i<=5;i++){
        ll a=rand()%(n-1)+1;
        ll x=Qc(a,m,n);
        ll y; 
        for(int k=1;k<=j;k++){
            y=mul(x,x,n);
            if(y==1&&x!=1&&x!=n-1)return 0;
            x=y;
        }
        if(x!=1)return 0;
    }
    return 1;
}
inline void dfs(ll x,ll y,int z){
    if(x+1>prime[num]&&Is_prime(x+1)) ans[++ans[0]]=y*(x+1);
    for(int i=z;i>=1;i--)
        if(x%(prime[i]-1)==0){
            ll a=x/(prime[i]-1),b=y,c=1;
            while(a%c==0){
                b*=prime[i];dfs(a/c,b,i-1);c*=prime[i];
            }
        }
    if(x==1){
        ans[++ans[0]]=y;return;
    }
}
int main()
{
//    freopen("arc.in","r",stdin);
//    freopen("arc.out","w",stdout);
    cin>>n>>k;srand(time(0));
    prim();dfs(n,1,num);
    sort(ans+1,ans+1+ans[0]);
    for(int i=1;i<=k;i++)
        cout<<ans[i]<<" ";
    return 0;
}

思路:

暴力筛法只有60分 
满分做法听得一脸懵逼
相信大家60分做法都会那就只贴std了

代码:

/*
    Problem:
        sigma(phi(i)),i<=N
    Algorithm:
        线性筛法+miu的高效计算 
    Note:
        O(N^(2/3)) 
*/
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define rep(i,x,y) for(i=x;i<=y;i++)

typedef long long int64;
const int MAXN=20000000;
int n,nn,i,j,mul;int64 ans;
bool pd[MAXN+1];int p[MAXN/5],M[MAXN+1];

int get(int n)
{
    if(n<=nn)return M[n];
    int i,j,ans=0;
    for(i=2;i<=n;i=j+1)
    {
        j=n/(n/i);
        ans+=get(n/i)*(j-i+1);   
    }
    return 1-ans;
}

int main()
{
    freopen("sum.in","r",stdin);
    freopen("sum.out","w",stdout);

    scanf("%d",&n);
    M[1]=1;nn=min(n,MAXN);
    for(i=2;i<=nn;i++)
    {
        if(!pd[i]){p[++p[0]]=i;M[i]=-1;}
        for(j=1;j<=p[0]&&(mul=p[j]*i)<=nn;j++)
        {
            pd[mul]=1;
            if(i%p[j]==0){M[mul]=0;break;}
            M[mul]=-M[i];
        }
    }
    rep(i,1,nn)M[i]+=M[i-1];

    for(i=1;i<=n;i=j+1)
    {
        j=n/(n/i);
        ans+=1LL*(get(j)-get(i-1))*(n/i)*(n/i+1)/2;
    }
    printf("%I64d\n",ans);
    return 0;
}