CodeForces 73 E.Morrowindows（数论）

Description

给出$n$个不同的正整数$a_i$，可以用这些数去测试一个数，即对于一个未知数$y$，可以得知$\lfloor\frac{y}{a_i}\rfloor$，问从这$n$个数中是否可以选出若干数字使得他们可以区分出$2$~$x$中所有的数，如果可以则要求选取最少数量的数，否则输出$-1$

Input

第一行两个整数$n,x$，之后输入$n$个整数$a_i(0\le n\le 10^5,2\le x\le 10^9,1\le a_i\le 10^9)$

Output

如果从这$n$个数中可以选出若干数字使得他们可以区分出$2$~$x$中所有的数，则输出所选数的最少数量，否则输出$-1$

Sample Input

1 4
2

Sample Output

-1

Solution

$x=2$时只有一个数字不用区分，答案是$0$，如果存在某个$a_i=1$则每个数都知道了也可以区分，答案是$1$

对于$2\le y<x$，要区分$y$和$y+1$，那么必然存在一个$a_i$使得$\frac{y}{a_i}<\frac{y+1}{a_i}$，由于$a_i>1$，故$a_i|y$，即对于任意$2$~$x-1$中任一数$y$，需存在一个$a_i$使得$a_i|y$，所以首先$2$~$x-1$中所有素数需要在$a_1,...,a_n$中， 且注意到只要有这些素数，必然可以区分出所有数字，缺一个素数$p$则区分不出$p$和$p+1$

由于$n\le 10^5$，而当$x\ge 1300000$时$2$~$x-1$里的素数个数已经超过$10^5$，故当$x>130000$时无解，否则直接一个个判断素数是否存在即可

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f,maxn=1300005;
int p[maxn],res;
bool vis[maxn];
void get_prime(int n=1300000)
{
    for(int i=2;i<=n;i++)
    {
        if(!vis[i])p[res++]=i;
        for(int j=0;j<res&&i*p[j]<=n;j++)
        {
            vis[i*p[j]]=1;
            if(i%p[j]==0)break;
        }
    }
}
int n,x,a[maxn],mark[maxn];
int main()
{
    get_prime();
    while(~scanf("%d%d",&n,&x))
    {
        memset(mark,0,sizeof(mark));
        int flag=0;
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
            if(a[i]==1)flag=1;
        }
        if(x==2)printf("0\n");
        else if(flag)printf("1\n");
        else if(x>=1300000)printf("-1\n");
        else
        {
            for(int i=1;i<=n;i++)
                if(a[i]<x)mark[a[i]]=1;
            int ans=0;
            for(int i=0;i<res&&p[i]<x;i++)
                if(!mark[p[i]])
                {
                    ans=-1;
                    break;
                }
                else ans++;
            printf("%d\n",ans);
        }
    }
    return 0;
}