1.题目描述
给定a[n]数组,重新排列a[n],使得最大,求最大值。
2.状态表示:
2.1目标:
若最终答案对应的数组a'[n],
gcd(a'[1]) = a'[1],记作i。
则最终答案 = dp[i] + c[i] * i;(c[i]为i的倍数的出现的次数)
2.2状态表示
参考2.1;
3.状态转移:
3.1转移形式:一个已知状态应该更新出哪些后续阶段的状态。
3.2具体解释:
已知状态:dp[i]
遍历后续阶段的状态:i 的倍数(除了i自己),记作j;
更新:dp[j] = max(dp[j], dp[i] + (c[i] - c[j]) * i);
4.代码
#include <bits/stdc++.h>
using namespace std;
#define endl "\n"
#define x first
#define y second
typedef pair<int,int> PII;
typedef long long LL;
typedef unsigned long long ULL;
const int mod = 1e9 + 7;
const int N = 5e6 + 10, M = 2 * N, P = 131;
int n;
int a[N];
LL dp[N];
int c[N];
void solve()
{
cin >> n;
int maxv = 0;
for(int i = 1; i <= n; i ++)
{
cin >> a[i];
maxv = max(a[i], maxv);
c[a[i]] ++;
}
for(int i = 1; i <= maxv; i ++)
for(int j = 2 * i; j <= maxv; j += i)
c[i] += c[j];
LL ans = 0;
for(int i = 1; i <= maxv; i ++)
for(int j = i; j <= maxv; j += i)
{
dp[j] = max(dp[j], dp[i] + (LL)(c[i] - c[j]) * (LL)i);
ans = max(ans, dp[j] + (LL)c[j] * j);
}
cout << ans << endl;
}
int main(){
std::ios::sync_with_stdio(false);std::cin.tie(0);
solve();}