JZOJ1241. Number

题目

Description

有N(2<=N<=15)个数A1,A2,….,An-1,An,如果在这N个数中，有且仅有一个数能整除m，那么整数m就是一个幸运数，你的任务就是在给定A1,A2,….,An-1,An的情况下，求出第k小的幸运数。

Input

第一行为一整数数N,K(2<=N<=15,1<=K<=2^31-1)，意义如上述。
接下来一行有N个整数，A1,A2,….,An-1,An,这N个整数均不超过2^31-1。

Output

输出一行，仅包含一个整数ans，表示第K小的幸运数。答案保证不超过10^15。

Sample Input

输入1：

2 4

2 3

输入2：

2 100

125 32767

Sample Output

输出1：

8

输出2：

12500

Hint

对于50%的数据，N<=5,ANS<=100000
对于80%的数据，N<=10,ANS<=10^15
对于100%的数据，N<=15,ANS<=10^15

分析

观察数据范围：
对于50%的数据是很简单的，只需要暴力。

而要通过后面的数据就有点小困难。

题解

很容易想到最后的答案肯定是A[]中间的某个数乘上一个整数。

看到K怎么大枚举每一个数很显然是不现实的。

我们考虑一下，可不可以把问题转换一下，
看一下可不可以变成一个判断问题。

• 看到K这么大，想到用二分。

现在问题就变成了求在区间[1..mid]有多少个符合条件的数。

• 有关整数倍数的问题很容易想到容斥原理。

code

#include <cstdio>#include <algorithm>#include <cstring>#include <string.h>#include <cmath>#include <stdlib.h>#include <math.h>#define ll long long using namespace std;ll a[20],l,r,ans,k,mid;int n,m;ll gcd(ll x,ll y){    if(x%y==0)return y;else return gcd(y,x%y);}void bfs(int x,ll sum,int deep){    if(sum>mid)return;    if(deep%2)ans+=floor(mid/sum)*deep;else ans-=floor(mid/sum)*deep;    for(int i=x+1;i<=n;i++)        bfs(i,sum/gcd(sum,a[i])*a[i],deep+1);}int main(){    scanf("%d%lld",&n,&k);    for(int i=1;i<=n;i++)        scanf("%lld",&a[i]);    l=1;    r=1000000000000000;    while(l<r)    {        mid=(l+r)/2;        ans=0;        for(int i=1;i<=n;i++)            bfs(i,a[i],1);        if(ans<k)l=mid+1;else r=mid;    }    printf("%lld\n",l);}