poj2429（miller_robin算法和pollard分解质因数）

/*
translation:
	给出两个数的gcd,lcm的值，求两个数分别是多少？有多解时输出两数和最小的一组。
solution:
	miller-robin和pollard素因子分解算法
	这道题可以参考poj1811（模板），根据规律a*b == gcd(a,b)*lcm(a，b)很容易想出分解质因数
	然后凑出a,b来暴力比较求出答案。但是由于数据范围太大。所以常规的素数分解算法不行。套上pollard
	算法的模板即可。然后dfs求出最靠近sqrt(y)但是不超过y的数字，这个数字就是所求的第一个数。
	注意一开始x,y相乘可能溢出，根据x,y的特性，可以用分解y/x的质因数来代替。因为y*x和y/x只有次数不同。
note:
	# 注意原来pollard模板不能对1分解质因数，会RE。在原来的基础上在一开始加上了对n==1的情况的处理即可
date:
	2016.10.26
*/
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <vector>
#include <ctime>
#include <utility>

using namespace std;
const int max_prime = 1e6;
const int s = 20;
typedef long long ll;

vector<ll> factor;
vector<ll> prime_factor;
ll x, y;	//x:gcd, y:lcm

ll multi_mod(ll a, ll b, ll n)
{
    a %= n;
    b %= n;
    ll res = 0;
    while(b)
    {