GCDLCM 【米勒_拉宾素数检验 (判断大素数)】

GCDLCM

题目链接（点击）

题目描述

In FZU ACM team, BroterJ and Silchen are good friends, and they often play some interesting games.
One day they play a game about GCD and LCM. firstly BrotherJ writes an integer A and Silchen writes an integer B on the paper. Then BrotherJ gives Silchen an integer X. Silchen will win if he can find two integers Y1 and Y2 that satisfy the following conditions:
• GCD(X, Y1) = A
• LCM(X, Y2) = B
• Fuction GCD(X, Y ) means greatest common divisor between X and Y .
• Fuction LCM(X, Y ) means lowest common multiple between X and Y .
BrotherJ loves Silchen so much that he wants Silchen to win the game. Now he wants to calculate how many number of X he can give to Silchen.

输入

Input is given from Standard Input in the following format:
A B
Constraints
1 ≤ A, B ≤ 1018
Both A and B are integers.

输出

Print one integer denotes the number of X.

样例输入

3 12

样例输出

题意：

给出A和B要求找到x、y1和y2满足条件的x个数：（y1和y2任意取）

gcd (x,y1) = A

gcd (x,y2) = B

例如当A=3 B=12时

当x=3的时候满足：gcd（3,3）= 3并且 lcm（3，12）

当x=6的时候满足：gcd（6,3）= 3并且 lcm（6，12）

当x=6的时候满足：gcd（12,3）= 3并且 lcm（12，3）

只有上面三种x符合条件所以结果是 3

思路：

可以列出表达式或者根据上面规律可以推出：

满足条件的 x 一定是A的倍数并且是B的因子也就是求B/A的因子个数（前提是A是B的倍数）

之前解决这类问题都是素数打表+唯一分解定理解决可是这个题的数范围太大打表基本不可能实现

就引入下面代码中的方法

米勒_拉宾素数检验：

具体推导过程和原理这个博客写的很好：

感谢：https://blog.csdn.net/qq_40564464/article/details/81774129

AC代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int mod=1e9+7;
LL qmul_mod(LL m,LL q,LL mod) ///快速乘 保证高精度 避免快速幂相乘爆LL
{LL ans=0;while(q){if(q%2){ans=(ans+m)%mod;}m=(m*2)%mod;q/=2;}return ans;
}
LL qpow_mod(LL m,LL q,LL mod) ///快速幂+取模
{LL ans=1;while(q){if(q%2){ans=qmul_mod(ans,m,mod);}m=qmul_mod(m,m,mod);q/=2;}return ans;
}
///2 7 61
bool Miller_Rabbin(LL x) ///米勒拉宾素数检验算法 博客里面给的是用2、3、5、7、11正确性很高 但我
{                            ///改了那个之后会Wa没办法 可用这个就过了 顿时…if(x==2||x==7||x==61) return true; //肯定是素数if(x%2==0||x%7==0||x%61==0) return false; //素数倍数if(qpow_mod(2,x-1,x)==1&&qpow_mod(7,x-1,x)==1&&qpow_mod(61,x-1,x)==1) ///均满足费马小return true;                      ///定理说明一定是素数（从推荐的那个博客中有讲原因）elsereturn false;
}
int main()
{LL A,B;scanf("%lld%lld",&A,&B);if(B%A==0){B/=A;LL ans=1;for(int i=2;i<=1000000;i++){ ///这种方法不需要素数打表 只要跑一边for循环就可以把小于LL Count=0;                  ///1e6的因子除去 还挺好while(B%i==0){Count++;B/=i;}ans=ans*(Count+1)%mod;}if(B>1){ ///下面可能有点难理解 if(Miller_Rabbin(B)){  ///被上面的循环筛过一次之后 剩下的如果B>1 剩下的肯定是由素ans=ans*2%mod;    ///数构成的数 可能是一个大素数 也可能是两个素数}                  ///而不可能是三个素数构成 因为三个>1e6的数相乘肯定会超1e18 else{LL t=sqrt(B);if(t*t==B) ///是两个相同的素数构成的ans=ans*3%mod; else ///两个不同的素数构成的ans=ans*4%mod;}}printf("%lld\n",ans);}elseprintf("0\n");return 0;
}