UPC 问题 D: GCDLCM （推理 + 处理技巧 + 米勒拉宾素数检测）

题目说 给A ， B 两个数，找x，x和 y1 的gcd 是 A， x 和y2的 lcm是 B   

问满足条件的x有几个

首先 A是 x 的因数， B是x的倍数，那就是 x=k*A,  B=q*x       B=A*k*q         B/A=k*q    设B/A = t 

把问题转化成了 求 k的个数  就是求 t的因子的个数 ，，因子个数？  很自然联想到 唯一分解定理。

 t的范围是 1e18的， 打 个 1e18的素数表？？？ 

其实 打个 1e6的素数表就行了，

假设 1e6后面还有素数，，那一定比 1e6 要大， 3个1e6 相乘 就是 1e18了

那就是  1e6 表 的 后面 ，最多会有两个超大素数。也可能没有

1e6 表的后面 的情况 

1，   前面除尽了，还剩1

2， 最后剩一个超大素数   p1（米勒拉宾素数检测 判断 这个数是不是素数）（待会说）

3，剩两个超大素数，但是他们相等  p1*p1

4， 剩两个超大素数，不相等，，p1*p2

 

 

怎么算因数的个数呢？

那就把 2的个数，3的个数，5的个数。。。。分别+1   再相乘

解释：  假设剩下两个2 一个3（可以不选2，，可以选1个2，，可以选2个2 ）（2的每一种情况对应  没有3，，有一个3。

取3 就这两种情况）  所以 每个素因子的个数  ，分别+1，   再相乘

假设剩的 是 1，不用管了

剩的是 p1  再乘2

剩p1*p1  就乘3

剩p1*p2 乘4

 

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1000006;


int phi[maxn]={0},prime[maxn];
int tot=0;
void Euler_on(int N){
    phi[1] = 1;
    for(int i = 2; i < N; i ++){
        if(!phi[i]){
            phi[i] = i-1;
            prime[tot++] = i;
        }
        for(int j = 0; j < tot && (1LL)*i*prime[j] < N; j ++){
            if(i % prime[j]) phi[i * prime[j]] = phi[i] * (prime[j]-1);
            else{
                phi[i * prime[j] ] = phi[i] * prime[j];
                break;
            }
        }
    }
}



ll add_mod(ll a,ll b,ll mod){
    ll ans=0;
    while(b){
        if(b&1)
            ans=(ans+a)%mod;
        a=a*2%mod;
        b>>=1;
    }
    return ans;
}

ll pow_mod(ll a,ll n,ll mod){
    if(n>1){
        ll tmp=pow_mod(a,n>>1,mod)%mod;
        tmp=add_mod(tmp,tmp,mod);
        if(n&1) tmp=add_mod(tmp,a,mod);
        return tmp;
    }
    return a;
}

bool Miller_Rabbin(ll n,ll a){
    ll d=n-1,s=0,i;
    while(!(d&1)){
        d>>=1;
        s++;
    }
    ll t=pow_mod(a,d,n);
    if(t==1 || t==-1)
        return 1;
    for(i=0;i<s;i++){
        if(t==n-1)
            return 1;
        t=add_mod(t,t,n);
    }
    return 0;
}

bool is_prime(ll n){
    ll i,tab[4]={3,4,7,11};
    for(i=0;i<4;i++){
        if(n==tab[i])
            return 1;
        if(!n%tab[i])
            return 0;
        if(n>tab[i] && !Miller_Rabbin(n,tab[i]))
            return 0;
    }
    return 1;
}

int cou[maxn];
int main(){

     Euler_on(1000000);

    ll A,B;
    cin>>A>>B;

    if(B%A!=0) cout<<0<<endl;
    else{
        int c=0;
        ll t=B/A;

        for(int i=1;i<tot;i++)
        {
            if(t%prime[i]==0)
                c++;
            while(t%prime[i]==0){
                t/=prime[i];
                cou[c]++;
            }
        }


        ll k=1;
            for(int i=1;i<=c;i++)
                k*=(cou[i]+1);


        if(t==1){
            cout<<k<<endl;

        }
        else if( is_prime(t) ){
            cout<<k*2<<endl;
        }
        else
        {
            ll tmp=sqrt(t);
            if( tmp*tmp==t )
            {
                cout<<k*3<<endl;
            }
            else
                cout<<k*4<<endl;
        }
    }
    return 0;
}

米勒拉宾素数检测的板子

咱也不知道，咱也不敢问，咱也不敢说啥，用就A了！

ll add_mod(ll a,ll b,ll mod){
    ll ans=0;
    while(b){
        if(b&1)
            ans=(ans+a)%mod;
        a=a*2%mod;
        b>>=1;
    }
    return ans;
}

ll pow_mod(ll a,ll n,ll mod){
    if(n>1){
        ll tmp=pow_mod(a,n>>1,mod)%mod;
        tmp=add_mod(tmp,tmp,mod);
        if(n&1) tmp=add_mod(tmp,a,mod);
        return tmp;
    }
    return a;
}

bool Miller_Rabbin(ll n,ll a){
    ll d=n-1,s=0,i;
    while(!(d&1)){
        d>>=1;
        s++;
    }
    ll t=pow_mod(a,d,n);
    if(t==1 || t==-1)
        return 1;
    for(i=0;i<s;i++){
        if(t==n-1)
            return 1;
        t=add_mod(t,t,n);
    }
    return 0;
}

bool is_prime(ll n){
    ll i,tab[4]={3,4,7,11};
    for(i=0;i<4;i++){
        if(n==tab[i])
            return 1;
        if(!n%tab[i])
            return 0;
        if(n>tab[i] && !Miller_Rabbin(n,tab[i]))
            return 0;
    }
    return 1;
}

每次只能等大佬讲完才有勇气做，，

自己太菜了