codeforces 354C C. Vasya and Beautiful Arrays( 数论+枚举)

题目链接：

[codeforces 354C](C. Vasya and Beautiful Arrays)

---

题目大意：

给出一个数组，每个元素最多减少k，问能够调整出的整个序列的最大的gcd是多少。

---

题目分析：

• 容易想到我们可能得到的答案的上限就是 $$\max_{i=1}^n a_{[i]}$$ .
• 因为通过减k,一定能够调整为k+1的倍数，所以下限是k+1
• 如果枚举答案，最坏依旧要枚举n次，所以我们如何能够快速判断当前答案是否正确呢？
• 我们可以知道每个数一定在$[i\cdot gcd, i \cdot gcd +k]( 1 \leq i \geq \frac{maxn}{gcd} )$,然后利用前缀和可以判断每个区间数的个数，那么我们可以枚举maxn/gcd次，判断这些区间的数的总数是不是数组中数的总数，如果是，那么就是当前这组case的解。
• 如果min <=k+1,那么一定可以所有数都可以通过不超过k的调整到是min的倍数，所以min就是答案。
• 总的复杂度是$\mathcal{O}(n \cdot log_2n)$(通过调和级数计算出)
  

  ---

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 );
    }
}