ARC136E
首先显然节点 1 必选,下面就忽略节点 1 了。
之后先考虑可达性,定义 f(n) 为 n 的最小质因数。
那么 x 能到 y 当且仅当满足以下四种条件之一:
- x≡0(mod2) 且 y≡0(mod2)
- x≡1(mod2) 且 y≡0(mod2),并且 x+f(x)≤y
- x≡0(mod2) 且 y≡1(mod2),并且 x≤y−f(y)
- x≡0(mod1) 且 y≡0(mod1),并且 x+f(x)≤y−f(y)
第一种情况显然。
对于节点 x 它能走到的最小可达点是 x+f(x),如果 x 是奇数,那么 x+f(x) 是偶数,所以只要 x+f(x)≤y 就能走到,这就是第二种情况。
同理 y−f(y) 是能走到 y 的最大可达点,对应第三种情况。
最后一种情况的充分性显然,下面证明必要性:
假设 y−f(y)<x+f(x),但是 x 依旧能走到 y。那么由上面的分析得 x 要是能走到 y 肯定是直达的。因此设 x=ac,y=bc,a,b,c 都是奇数,且 a<b,c>1。
那么 c 显然是不小于 f(x) 和 f(y) 的,那么有 x+f(x)≤x+c,y−c≤y−f(y),因为 a<b 且 a,b 均是奇数,所以 b−a≥2,所以 x+c≤y−c,所以 x+f(x)≤y−f(y),与假设矛盾。证毕。
由上面的分析可以想到构造区间
- x≡0(mod2):l(x)=r(x)=x
- x≡1(mod2):l(x)=x−f(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;
}
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列01:轻松3步本地部署deepseek,普通电脑可用
· 按钮权限的设计及实现
· 25岁的心里话