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);
}
做了无用的快读快输, 其实去掉也可(