题目链接:
[codeforces 354C](C. Vasya and Beautiful Arrays)
题目大意:
给出一个数组,每个元素最多减少k,问能够调整出的整个序列的最大的gcd是多少。
题目分析:
- 容易想到我们可能得到的答案的上限就是
maxi=1na[i].
- 因为通过减k,一定能够调整为k+1的倍数,所以下限是k+1
- 如果枚举答案,最坏依旧要枚举n次,所以我们如何能够快速判断当前答案是否正确呢?
- 我们可以知道每个数一定在 [i⋅gcd,i⋅gcd+k](1≤i≥maxngcd) ,然后利用前缀和可以判断每个区间数的个数,那么我们可以枚举maxn/gcd次,判断这些区间的数的总数是不是数组中数的总数,如果是,那么就是当前这组case的解。
- 如果min <=k+1,那么一定可以所有数都可以通过不超过k的调整到是min的倍数,所以min就是答案。
- 总的复杂度是
O(n⋅log2n)
(通过调和级数计算出)
AC代码:
#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#define M 1000007
using namespace std;
int sum[M],a,n,k;
int main ( )
{
while ( ~scanf ("%d%d" , &n , &k ))
{
memset ( sum , 0 , sizeof ( sum ) );
int maxn = -1 ,minn = M;
for ( int i = 0 ; i < n ; i++ )
{
scanf ("%d" , &a );
sum[a]++;
maxn = max ( maxn , a );
minn = min ( minn , a );
}
for ( int i = 1 ; i <= maxn ; i++ )
sum[i]+= sum[i-1];
int i;
if ( minn <= k+1 )
{
printf ( "%d\n" , minn );
continue;
}
for ( i = maxn ; i >= k+1 ; i-- )
{
int m = 0;
for ( int j = 1 ; j <= maxn/i ; j++ )
{
int l = j*i;
int r = min(j*i+k,maxn);
m += sum[r] - sum[l-1];
}
if ( m == n ) break;
}
printf ( "%d\n" , i );
}
}