常见部分数学问题

最大公约数与最小公倍数

欧几里得算法

int gcd(int a, int b){
    if(b == 0) return a;
    else return gcd(b, a % b);
}

更简洁的写法

int gcd(int a, int b){
    return !b ? a : gcd(b, a % b);
}

由最大公约数可以很快得到最小公倍数即ab/d，而ab在实际计算时有可能溢出，因此更恰当的写法是a/d*b。

素数问题

• 素数的判断

若在2~n-1中存在n的约数，不妨设为k，则必存在另一个约数n/k，且一定满足一个小于sqrt(n)、另一个大于sqrt(n)。所以只要判断n能否被2~sqrt(n)中的一个整除，即可判定n是否为素数，算法复杂度为O(sqrt(n))。

代码如下：

bool isPrime(int n){
    if(n <= 1)return false;
    int sqr = (int)sqrt(1.0 * n);
    for(int i = 2;i <= sqr;i++){
        if(n % i == 0)return false;
    }
    return true;
}

如果n没有接近int型变量范围上界，则可将i <= sqr写成 i*i <= n，当然还可以将i定义为long long型。

• 素数筛法

Eratosthenes筛法，复杂度为O(nloglogn)，算法从小到大枚举所有数，对每一个素数，筛去它的所有倍数。

const int main = 101;
int prime[maxn],pNum = 0;
bool p[maxn] = {0};
void Find_Prime(){
    for(int i = 2;i < maxn;i++){
        if(p[i] == false){
            prime[pNum++] = i;//把素数i存到prime数组中
            for(int j = i + i;j < maxn;j += i){
                //筛去所有i的倍数
                p[j] = true;
            }
        }
    }
}

质因子分解

对一个正整数n来说，如果它存在[2,n]范围内的质因子，要么这些质因子全部小于等于sqrt(n)，要么只存在一个大于sqrt(n)的质因子。

存放质因子及其个数的结构体：

struct factor{
    int x,cnt;
}fac[10];

考虑到质因子乘到29就已经超过了int范围，因此fac数组开10就可以。

if(n % prime[i] == 0){
    fac[num].x = prime[i];
    fac[num].cnt = 0;
    while(n % prime[i] == 0){
        fac[num].cnt++;
        n /= prime[i];
    }
    num++;
}

如果上述步骤结束后n仍然大于1，说明n有且仅有一个大于sqrt(n)的质因子：

if(n != 1){
    fac[num].x = n;
    fac[num++].cnt == 1;
}

例 PAT A1059

#include <cstdio>
#include <cstring>
#include <cstring>
#include <algorithm>
#include <cstdlib>
#include <math.h>
using namespace std;

const int maxn = 100010;

struct factor{
    int x,cnt;
}fac[10];

bool is_prime(int n){
    if (n == 1)return false;
    int sqr = (int)sqrt(1.0 * n);
    for (int i = 2; i <= sqr; i++) {
        if(n % i == 0)return false;
    }
    return true;
}

int prime[maxn],pNum = 0;

void find_prime(){
    for(int i = 1;i < maxn;i++){
        if(is_prime(i) == true){
            prime[pNum++] = i;
        }
    }
}

int main(){
    int input,num = 0;
    find_prime();
    scanf("%d",&input);
    if(input == 1){
        printf("1=1");
        return 0;
    }
    printf("%d=",input);
    for(int i = 0;prime[i] <= sqrt(input * 1.0) && i < pNum;i++){
        if(input % prime[i] == 0){
            fac[num].x = prime[i];
            fac[num].cnt = 0;
            while (input % prime[i] == 0){
                fac[num].cnt++;
                input /= prime[i];
            }
            num++;
        }
        if(input == 1)break;
    }
    if(input != 1){
        fac[num].x = input;
        fac[num++].cnt = 1;
    }
    for(int i = 0;i < num;i++){
        printf("%d",fac[i].x);
        if(fac[i].cnt > 1){
            printf("^%d",fac[i].cnt);
        }
        if(i+1 < num){
            printf("*");
        }
    }
    return 0;
}