洛谷2018-7月月赛

U29473

【题目链接】

　　　　https://www.luogu.org/problemnew/show/U27076

【算法】

　　　　我的算法比较渣，最多n=1e5个数，于是离散化，记录每个数出现次数。

　　　　然后判断是否是质数，每个数范围为1~1e12于是只需判断是否能能被1~1e6范围内且小于它的1/2次方的质数整除即可，所以打了一个1~1e6的质数表。

　　　　看了标程自己是个zz，根号ai不超过1e6，没必要大表，不过快一点。。。。而且，排个序判断下就好，实乃水题，溜了，我太渣。。。。

【注意】

　　　　stl::unique() 函数如果需要计算重复数字的话尽量不要用。。。

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e6+10;
typedef long long ll;
vector<int> a;
bool is_prime[maxn];
int t,m,n,tot;
ll b[100010],c[100010];
int d[100010];
void prime()
{
    is_prime[2] = true;
    for(int i = 3; i < maxn; i += 2)
        is_prime[i] = true;
    for(int i = 3; i < sqrt(maxn+0.1); i++) {
        if(is_prime[i]) {
            for(int j = i*i; j < maxn; j += 2 * i)
                is_prime[j] = false;
        }
    }
    for(int i = 2; i < maxn; i++) if(is_prime[i]) a.push_back(i);
}
void discrete()
{
    sort(b+1,b+n+1);
    for(int i = 1; i <= n; i++) {
        if(b[i] != b[i-1]) c[++tot] = b[i];
        d[tot]++;
    }
}
int query(ll x)
{
    return lower_bound(c+1,c+tot+1,x) - c;
}
int main()
{
    prime();
    cin>>t;
    while(t--) {
        memset(d,0,sizeof(d));
        tot = 0;
        cin>>n>>m;
        for(int i = 1; i <= n; i++) {
            cin>>b[i];
        }
        discrete();
        for(int i = 1; i <= m; i++) {
            ll cur;
            cin>>cur;
            d[query(cur)]--;
        }

        ll times=0, rec=0;
        for(int i = 1; i <= tot; i++) {
            if(d[i] > 0 && c[i] != 1) {
                times += d[i];
                rec = c[i];
            }
        }
        if(times != 1)
            cout<<"NO"<<endl;
        else {
            bool flag = true;
            for(int i = 0; a[i] <= sqrt(rec+1.0) && i < a.size(); i++)
                if(rec % a[i] == 0) { flag = false; break; }
            if(flag) cout<<"YES"<<endl;
            else cout<<"NO"<<endl;
        }
    }
    return 0;
}