2021-05-25acwing分解质因数，快速幂

AcWing 197. 阶乘分解

1.题意：

2.题解：

之前做过，
正常会有以下几种想法：
1，计算n！，再筛，—显然不行n!太大了，而且没有mod
2，将n！遍历，得出每个数的质因子，将对应的质因子的指数++，枚举每个数的所有质因子大概是根号n，所以总时间复杂度是n根号n，1e9显然也不行
so 换种思路，遍历所有质数，将n！中所有质数的倍数个数相加
比如p1,首先是n中p1的倍数，然后是p1²的倍数。。。相加得到n！中p1的指数

3.ac代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;
bool st[N];
int primes[N];
int tot;
void get_primes(int n){
    for(int i=2;i<=n;i++){
        if(!st[i]) primes[tot++]=i;
        for(int j=0;primes[j]<=n/i;j++){
            st[primes[j]*i]=1;
            if(i%primes[j]==0) break;
        }
    }
}
int main(){
    int n;
    cin>>n;
    get_primes(N);
    for(int i=0;i<tot;i++){
        int t=n;
        int s=0;
        while(t){
            s+=t/primes[i];
            t/=primes[i];
        }
        if(s!=0)
        cout<<primes[i]<<" "<<s<<endl;
    }
    
}

AcWing 1289. 序列的第k个数

1.题意：

2.题解：

快速幂如果参数和返回值都是int，中间计算值记得带longlong
给出前三项一定可以判断是等差还是等比，等值数列都一样

3.ac代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;
const int mod= 200907 ;
int qpow(int a,int b,int mod){
   int res=1;
   while(b){
       if(b&1) res=(ll)res*a%mod;
       b>>=1;
       a=(ll)a*a%mod;
   }
   return res%mod;
}
int main(){
    int t,k,a,b,c;
    cin>>t;
    while(t--){
        cin>>a>>b>>c>>k;
        if(a+c==2*b){
            
            cout<<((ll)a+(ll)(c-b)*(k-1)%mod)%mod<<endl;
        }else{
            
            cout<<(ll)a*qpow(b/a,k-1,mod)%mod<<endl;
        }
    }
}

AcWing 1290. 越狱

1.题意：

2.题解：

记得在每一次可能超的计算都余除mod，因为余除mod后就不一定谁大谁小，所以最后结果加上mod再取余，使得最终结果是正数

3.ac代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;
const int mod=100003  ;
ll n,m;
ll qpow(ll a,ll b){
    ll res=1;
    while(b){
        if(b&1) res=res*a%mod;
        b>>=1;
        a=a*a%mod;
    }
    return res%mod;
}
int main(){
    cin>>m>>n;
    cout<<(qpow(m,n)-m*qpow(m-1,n-1)%mod+mod)%mod<<endl;
}