Codeforces Round #757 (Div. 2) D1 + D2. Divan and Kostomuksha

D1

给定一个数组a，对其重新排序，求

的最大值。
思路（看了题解来的：如果某个gcd值为x，为了让和最大，应该让x的倍数都在x前面。
因此，可以枚举因子进行dp。枚举每个数的 gcd值，前面一个gcd值是它的倍数。
代码：

#include<bits/stdc++.h>
using namespace std;
#define int long long
#define mem(a, b) memset(a, b, sizeof a)

const int N = 5e6;
int n, cnt[N + 5], f[N + 5];

int read(){
    int x = 0, f = 1;
    char c = 0;
    while(!isdigit(c)){if(c == '-') f = -1; c = getchar(); }
    while(isdigit(c)){x = x * 10 + c - '0'; c = getchar(); }
    return x * f;
}
signed main(){
    n = read();
    for(int i = 1; i <= n; i ++){
        int u;
        cnt[u = read()] ++;
    }
    for(int i = 1; i <= N; i ++)
        for(int j = i * 2; j <= N; j += i)
            cnt[i] += cnt[j];
    f[1] = n;
    for(int i = 1; i <= N; i ++)
        for(int j = i * 2; j <= N; j += i)
            f[j] = max(f[j], f[i] + (j - i) * cnt[j]);
    cout << *max_element(f + 1, f + N + 1) << endl;
}

D2

和 D1 相同，但是范围到了2e7。

思路：在进行dp时，上一个倍数j，j / i如果不是素数，那么不是最优的，因为它们中间还可以塞进去其他i的倍数，使得和更大，所以在进行dp时，只需枚举i的素数倍，这样就能卡常卡过去了。

注意不要#define int long long，会变得不幸。

代码：

#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define mem(a, b) memset(a, b, sizeof a)

const int N = 2e7 + 10;
int n;
ll cnt[N], f[N];
bool vis[N];
int primes[N], k;

void get_primes(){
    for(int i = 2; i < N; i ++){
        if(!vis[i]) primes[++ k] = i;
        for(int j = 1; j <= k && i * primes[j] < N; j ++){
            vis[i * primes[j]] = 1;
            if(i % primes[j] == 0) break;
        }
    }
}
inline int read(){
    ll x = 0, f = 1;
    char c = 0;
    while(!isdigit(c)){if(c == '-') f = -1; c = getchar(); }
    while(isdigit(c)){x = x * 10 + c - '0'; c = getchar(); }
    return x * f;
}
inline void write(ll n){
    if(n < 0){ putchar('-'); n = -n; } 
    if(n > 9) write(n / 10);
    putchar(char(n % 10 + '0'));
}
signed main(){
    get_primes();
    n = read();
    for(int i = 1; i <= n; i ++){
        int u = read();
        for(int j = 1; j * j <= u; j ++)
            if(u % j == 0){
                cnt[j] ++;
                if(j != u / j) cnt[u / j] ++;
            }
    }
    for(int i = N - 5; i; i --){
        f[i] = 1LL * i * cnt[i];
        for(int j = 1; j <= k && i * primes[j] < N; j ++)
            f[i] = max(f[i], f[i * primes[j]] + 1LL * (cnt[i] - cnt[i * primes[j]]) * i);
    }
    ll ans = 0;
    for(int i = 1; i < N; i ++) if(cnt[i] == n) ans = max(ans, f[i]);
    write(ans);
}

做了无用的快读快输， 其实去掉也可（