技巧：机试常用算法

一、素数判断：
1、直接判断 o(n^1/2)

bool isprime( int n )
{
    if( n < 2 )
        return false;
    for( int i = 2; i <= sqrt(n); ++i )
        if( n % i == 0 )
            return false;
    return true;
}

2、埃氏筛打标 0~n的素数o(nloglogn)

    vector<bool> isprime(n + 1, true);
    isprime[0] = isprime[0] = false;
    for( int i = 2; i <= n; ++i )
        if( isprime[i] )
            for(int j = 2; i * j <= n; ++j )
                isprime[i * j] = false;

二、最大公因数、最小公倍数
1、辗转相除法
注意：
（1）在分数化简等领域的应用中，我们总是希望求得最大公约数为一个正数，从而实现不改变数符的前提下的化简，因此，尽量对传入参数调用abs()函数
（2）应该无须要求a>b，具体暂时没空证明，但是应该是这样(摸鱼爽)

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

2、最小公倍数

int lcm(int a, int b)
{
	return a * b / gcd(a, b);
}

三、因式分解
1、质因数分解

struct node
{
	int f, e;
};
vector<node> factor;
r = n;
for( int i = 2; i <= sqrt(n) && r != 1; ++i )
	if( r % i == 0 )
	{
		node temp = {i, 0};
		while( r % i == 0 )
		{
			r /= i;
			++temp.e;
		}
		factor.push_back(temp);
	}
if( r != 1 )
	factor.push_back({r, 1});

2、求全部因数

vector<int> factor;
int i;
for( i = 1; i * i < n; ++i )
	if( n % i == 0 )
	{
		factor.push_back(i);
		factor.push_back(n / i);
	}
if( i * i == n )
	factor.push_back(i);

3、求n!的质因数

    vector<bool> isprime(n + 1, true);
    isprime[0] = isprime[0] = false;
    for( int i = 2; i <= n; ++i )
        if( isprime[i] )
            for(int j = 2; i * j <= n; ++j )
                isprime[i * j] = false;
	vector<int> cnt_f(n + 1, 0);
	for( int i = 2, temp; i <= n; ++i )
		if( isprime[i] )
		{
			temp = n;
			while(temp)
			{
				cnt_f[i] += temp / i;
				temp /= i;
			}
		}

三、组合数求解
1、分解+打表 C(k, n) = C(k, n - 1) + C(k - 1, n - 1)

vector<vector<int>> table(n + 1, vector<int>(n + 1, 0))
int C(int k, int n)
{
	if( !k || k == n )
		return 1;
	if( !table[k][n] )
		table[k][n] = C(k, n - 1) + C(k - 1, n - 1);
	return table[k][n];
}

2、定义变形

int C(int k, int n)
{
	int ans = 1;
	for( int i = 1; i <= k; ++i )
		ans = ans * (n - k + i) / i;
	return ans;
}

四、快速幂
1、递归 o(logk)

int bpow( int n, int k )
{
	if( !k )
		return 1;
	if( k % 2 )
		return bpow(n, k - 1) * n;
	int temp = bpow(n, k / 2);
	return temp * temp; 
}

2、迭代o(logk)

int bpow( int n, int k )
{
	int ans = 1;
	for( ; k; k /= 2, n *= n )
		if( k % 2 )
			ans *= n;
	return ans;
}

五、DFS
1、注意剪枝

六、Dijkstra
1、存在零权负权边时注意，距离更新时排除距离已确定的顶点
2、异值重复边取最小值

七、大整数
1.注意乘法溢出
2.注意除数为零