HPU272 wzy的大冒险——出发前的清理（米勒罗宾板子）

                         272. wzy的大冒险——出发前的清理

单点时限: 1.0 sec

内存限制: 512 MB

由于上次学弟们没有ak，导致许多蚂蚁被留下了。wzy在出发冒险前请来了一只食蚁兽帮忙清理。
现在出现了一只食蚁兽。每个蚂蚁都有wzy给它的一个编号，食蚁兽要吃蚂蚁必须要确认蚂蚁的编号X 是否满足要求，如下：

1. X 的质因子的种类不超过13种
2. X 的食蚁数是个素数(食蚁数的定义见最下方提示)
3. X 的食蚁数是个回文数

如果都满足则输出 YES ,否则输出 NO

输入格式

第一行给定一个T ,(1≤T≤1e5)
接下来T 行，每行给出数字X ,(1≤X≤1e17)

输出格式

每行输出 YES 或者 NO

样例

Input

3
121666
111666
1312333

Output

NO
NO
YES

提示

补充：
1. 回文数至少为两位数，如 131 、 22
2. 如果X 的长度为 len1 ，令 len2=⌊len12⌋ ，则 X 的前 len2 位形成的数是它的食蚁数

 

以前写的米勒罗宾板子过不了，然后找了找其他人的博客过的：https://blog.csdn.net/lbperfect123/article/details/85057358

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define INF 0x3f3f3f3f
const int MAXN = 1e5 + 7;

ll prime[6] = {2, 3, 5, 233, 331};
ll qmul(ll x, ll y, ll mod) 
{ // 乘法防止溢出， 如果p * p不爆ll的话可以直接乘； O(1)乘法或者转化成二进制加法
    return (x * y - (long long)(x / (long double)mod * y + 1e-3) *mod + mod) % mod;
    /*
	ll ret = 0;
	while(y) {
		if(y & 1)
			ret = (ret + x) % mod;
		x = x * 2 % mod;
		y >>= 1;
	}
	return ret;
	*/
}
ll qpow(ll a, ll n, ll mod) 
{
    ll ret = 1;
    while(n) 
    {
        if(n & 1) ret = qmul(ret, a, mod);
        a = qmul(a, a, mod);
        n >>= 1;
    }
    return ret;
}
bool Miller_Rabin(ll p)
{
    if(p < 2) return 0;
    if(p != 2 && p % 2 == 0) return 0;
    ll s = p - 1;
    while(! (s & 1)) s >>= 1;
    for(int i = 0; i < 5; ++i) 
    {
        if(p == prime[i]) return 1;
        ll t = s, m = qpow(prime[i], s, p);
        while(t != p - 1 && m != 1 && m != p - 1)
        {
            m = qmul(m, m, p);
            t <<= 1;
        }
        if(m != p - 1 && !(t & 1)) return 0;
    }
    return 1;
}

int main()
{
    ios::sync_with_stdio(false);
    int t;
    cin>>t;
    string s1,s2;
    while(t--)
    {
        cin>>s1;
        s2 = s1 = s1.substr(0,s1.size()/2);
        reverse(s2.begin(), s2.end());
        if(s1 != s2 || s1.size() < 2)
        {
            cout<<"NO"<<endl;
            continue;
        }
        ll x = stoi(s2);//字符串转换为数字
        if(Miller_Rabin(x))
            cout<<"YES"<<endl;
        else
            cout<<"NO"<<endl;
    }

    return 0;
}

大佬说这道题用罗宾做就是暴力，正解应该是打表，这个等讲完课再补吧。。。