标签：cnt Kostomuksha int ll CF long Divan 复杂度 dp

• 题意：NKOJ
  CF
• 思路：首先发现贪心不了。因此dp。然后这题需要维护的就\(g_i\)和\(sum{g_i}\)
  状态：\(dp[i]\): 当前最后一个为\(g_i\)的最大值
  \(dp[i]= \max_{i|j}(dp[j]+(cnt[i]-cnt[j])*i)\)
  \(cnt[i]\): \(a[]\)中\(i\)的倍数的个数
  \(cnt[i]= \sum_{i|j}cnt[j]\)这个可以用狄利克雷后缀和
  然后复杂度：\(O(nlogn)\)其中\(n=10^7\),然后我们可以证明枚举倍数的时候乘上质数因此复杂度可以除以一个\(ln(n)\)时限是4s就可以过了.
  当然稳一点的话，就反过来枚举因数。这样可以\(O(nlogn)\)其中\(n=10^5\)
• code:

#include<bits/stdc++.h>
using namespace std;
const int N=2e7+5;
typedef long long ll;
bool is_p[N];
int a[N/100],cnt[N],up=2e7,p[N>>2],ptot;
ll dp[N];
int main() {
	int n;scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]),cnt[a[i]]++;
	for(int i=2;i<=up;i++) {
		if(!is_p[i]){p[++ptot]=i;}
		for(int j=1,x;j<=ptot&&i*p[j]<=up;j++) {
			is_p[i*p[j]]=1;
			if(i%p[j]==0)break;
		}
	}
	for(int i=1;i<=ptot;i++) {
		for(int j=up/p[i];j>=1;j--) cnt[j]+=cnt[j*p[i]];
	}
	ll ans=0;
	for(int i=up;i>=1;i--) {
		if(!cnt[i])continue;
		dp[i]=0;for(int j=1,x;(x=i*p[j])<=up;j++) if(cnt[x])dp[i]=max(dp[i],dp[x]-1ll*cnt[x]*i);
		dp[i]+=1ll*i*cnt[i];ans=max(ans,dp[i]);
	}
	printf("%lld\n",ans);
	return 0;
}

标签：cnt,Kostomuksha,int,ll,CF,long,Divan,复杂度,dp
来源： https://www.cnblogs.com/bestime/p/15846138.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。