ARC136E

首先显然节点 1 必选,下面就忽略节点 1 了。

之后先考虑可达性,定义 f(n)n 的最小质因数。

那么 x 能到 y 当且仅当满足以下四种条件之一:

  • x0(mod2)y0(mod2)
  • x1(mod2)y0(mod2),并且 x+f(x)y
  • x0(mod2)y1(mod2),并且 xyf(y)
  • x0(mod1)y0(mod1),并且 x+f(x)yf(y)

第一种情况显然。

对于节点 x 它能走到的最小可达点是 x+f(x),如果 x 是奇数,那么 x+f(x) 是偶数,所以只要 x+f(x)y 就能走到,这就是第二种情况。

同理 yf(y) 是能走到 y 的最大可达点,对应第三种情况。

最后一种情况的充分性显然,下面证明必要性:

假设 yf(y)<x+f(x),但是 x 依旧能走到 y。那么由上面的分析得 x 要是能走到 y 肯定是直达的。因此设 x=ac,y=bca,b,c 都是奇数,且 a<b,c>1

那么 c 显然是不小于 f(x)f(y) 的,那么有 x+f(x)x+cycyf(y),因为 a<ba,b 均是奇数,所以 ba2,所以 x+cyc,所以 x+f(x)yf(y),与假设矛盾。证毕。

由上面的分析可以想到构造区间

  • x0(mod2)l(x)=r(x)=x
  • x1(mod2)l(x)=xf(x)+1,r(x)=x+f(x)1

那么 x 走不到 y 就说明 [l(x),r(x)][l(y),r(y)],那么差分一下,枚举交集是哪个点即可。

Code:

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 1000000;
int n;
int a[N+5]; ll sum[N+5];
bool vis[N+5];
int pri[N], tot;
int mn[N+5];

void init() {
	for (int i = 2; i <= N; ++i) {
		if (!vis[i]) pri[++tot] = i, mn[i] = i;
		for (int j = 1; j <= tot && i * pri[j] <= N; ++j) {
			vis[i * pri[j]] = 1;
			mn[i * pri[j]] = pri[j];
			if (i % pri[j] == 0) break;
		}
	}
}

int main() {
	init();
	scanf("%d", &n);
	for (int i = 1; i <= n; ++i) scanf("%d", &a[i]);
	for (int i = 1; i <= n; ++i) {
		int l, r;
		if (i == 1) l = 1, r = n;
		else if (i % 2 == 0) l = r = i;
		else l = i - mn[i] + 1, r = i + mn[i] - 1;
		sum[max(1, l)] += a[i], sum[min(n + 1, r + 1)] -= a[i];
	}
	for (int i = 1; i <= n; ++i) sum[i] += sum[i - 1];
	ll ans = 0;
	for (int i = 1; i <= n; ++i) ans = max(ans, sum[i]);
	printf("%lld", ans);
	return 0;
}
posted @   Kobe303  阅读(28)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列01:轻松3步本地部署deepseek,普通电脑可用
· 按钮权限的设计及实现
· 25岁的心里话
点击右上角即可分享
微信分享提示