#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int P=1e9+7,M=1e6+5;
ll exgcd(ll a,ll b,ll &x,ll &y){
if(!b){x=1;y=0;return a;}
ll t=x-(a/b)*y;
return exgcd(b,a%b,y,t);
}
inline ll gcd(ll x,ll y){
while(y)y^= x^= y ^=x %=y;
return x;
}
//快速幂
inline int fast(int b,int n){
ll a=b,ans=1;
while(n){
if(n&1)ans=(ans*a)%P;
a=1ll*a*a%P;
n>>=1;
}
}
//快速乘
inline int Pre(ll x,ll n,ll mod){
int res=0;
x%=p;
while(n){
if(n&1)res=x+res%mod;
x=(x+x)%p;
n>>=1;
}return res;
}
//我们要保证每个数只被它最小的质因子标记一次。
int check[M],prime[M],num,ans[M];
inline void init(int n){
memset(check,0,sizeof(check));
for(int i=2;i<=n;i++){
if(!check[i])prime[++num]=i;
for(int j=1;j<=num;j++){
if(prime[j]*i>n)break;
//记录最小质因子
check[i*prime[j]]=prime[j];
if(i%prime[j]==0)break;
}
}
}
//1e9因数分解
void solve_1e9_div(int n){
while(check[n]){
ans[++num]=check[n];
n/=check[n];
}ans[++num]=n;
}
//对于1e18的质数判定
//利用费马小定理
//对于互质的整数x和质数p,
//x^p-1=1(mod p)
//但是它有一定的局限性 如341
//首先判掉2和偶数。
//我们可以考虑这个式子x^2≡1 (mod n)
//x^2-1 ≡0 (mod n)
//(x-1) (x+1) = kn
//由于0≤x <n如果n是一个质数,显然只能(x-1)|n或者(x+1)|n,也就是x=1,n-1。
//所以如果有一个x不是这两个数的解,则n一定不是质数。
//对于10^18以内的数,x取前9个质数就可以准确无误的判断了。
bool is_inv(int x,int p){
ll a=x,u=p-1,res=1;
while(u){
ll nw=res;
ll nx=nw*nw%p;
if(nx==1&&nw!=1&&nw!=p-1)return 0;
//x^2-1 ≡0 (mod n)
//(x-1) (x+1) = kn
//x=1,n-1;
if(u&1)res=res*a%p;
a=a*a%p;
u>>=1;
}
return res==1;//x^(p-1)%p==1
}
bool miller_rabin(ll n){
if(n==2)return 1;
if(n<2||n%2==0)return 0;
//非常大的优化
for(int i=1;prime[i]<n&&i<=20;i++){
if(!is_inv(prime[i],n))return 0;
}return 1;
}
//所以如果我们在(mod m)的意义下,随机生成一个序列。
//当序列长度为√m的时候,就有两个数相同。(生日悖论)
//如果我们要分解n,对于n的最小质因子,假设为p,那么在√p步之内,很可能找到两个数在(mod p)下同余。
//如果我们找到了两个数x_1和x_2,他们在模p下同余,却在模n下不同余,那么我们求 gcd(n,x_1-x_2),就找到一个质因子。
map<ll,int>res;
map<ll,int>::iterator it;
ll pollard_rho(ll n,ll c){
ll x,y,d,i=1,k=2;
x=rand()%n+1;
y=x;
while(1){
i++;
x=(Pre(x,x,n)+c)%n;
//x[i+1]=x[i]*x[i]+c;来充当随机数组
d=gcd(y-x,n);
if(1<d&&d<n)return d;//找到了
if(y==x)return n;//没有用
if(i==k)y=x,k<<=1;//找到第2^k项为 y
}
}
void div_(ll n,ll c){
if(n==1)return;
if(miller_rabin(n)){
res[n]++;
return;
}
ll p=n;
while(p==n)p=pollard_rho(p,c--);
div_(p,c);
div_(n/p,c);
}
void solve_1e18_div(ll n){
init(1e6);
srand((int)time(0));
res.clear();
div_(n,(int)rand()+1e9);
if(res.empty())printf("%lld\n",n);
for(it= res.begin(); it!= res.end();){
printf( "%lld^%d",it->first,it->second);
if((++it) != res.end())printf( " * " );
}
}
简单数学(关于素数判断和因数分解)
最新推荐文章于 2023-05-06 20:53:42 发布