Ultimate Weirdness of an Array

写不出来, 日常好菜啊。。

考虑枚举GCD, 算出一共有多少个对 f(l, r) <= GCD, 我们用fuc[ i ] 表示的是在 l = i 这个位置开始, 最小的合法的 R,

可以发现这个函数随 i 单调不下降, 枚举GCD 的时候, 找到GCD 的倍数的位置, 用线段树更新最大值。

#include<bits/stdc++.h>
#define LL long long
#define LD long double
#define ull unsigned long long
#define fi first
#define se second
#define mk make_pair
#define PLL pair<LL, LL>
#define PLI pair<LL, int>
#define PII pair<int, int>
#define SZ(x) ((int)x.size())
#define ALL(x) (x).begin(), (x).end()
#define fio ios::sync_with_stdio(false); cin.tie(0);using namespace std;const int N = 2e5 + 7;
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int mod = 998244353;
const double eps = 1e-8;
const double PI = acos(-1);template<class T, class S> inline void add(T& a, S b) {a += b; if(a >= mod) a -= mod;}
template<class T, class S> inline void sub(T& a, S b) {a -= b; if(a < 0) a += mod;}
template<class T, class S> inline bool chkmax(T& a, S b) {return a < b ? a = b, true : false;}
template<class T, class S> inline bool chkmin(T& a, S b) {return a > b ? a = b, true : false;}int n, a[N];
int maxPos[N][2];
int minPos[N][2];
int cnt[N];
LL H[N];
int mx[2], mn[2], tot;#define lson l, mid, rt << 1
#define rson mid + 1, r, rt << 1 | 1
struct segmentTree {LL sum[N << 2]; int lazy[N << 2], mn[N << 2], mx[N << 2];inline void pull(int rt) {sum[rt] = sum[rt << 1] + sum[rt << 1 | 1];mx[rt] = max(mx[rt << 1], mx[rt << 1 | 1]);mn[rt] = min(mn[rt << 1], mn[rt << 1 | 1]);}inline void push(int rt, int l, int r) {if(~lazy[rt]) {int mid = l + r >> 1;sum[rt << 1] = 1LL * (mid - l + 1) * lazy[rt];sum[rt << 1 | 1] = 1LL * (r - mid) * lazy[rt];lazy[rt << 1] = lazy[rt << 1 | 1] = lazy[rt];mx[rt << 1] = mn[rt << 1] = lazy[rt];mx[rt << 1 | 1] = mn[rt << 1 | 1] = lazy[rt];lazy[rt] = -1;}}void build(int l, int r, int rt) {lazy[rt] = -1;if(l == r) {sum[rt] = l;mx[rt] = l;mn[rt] = l;return;}int mid = l + r >> 1;build(lson); build(rson);pull(rt);}void update(int L, int R, int val, int l, int r, int rt) {if(R < l || r < L || R < L) return;if(mn[rt] >= val) return;if(L <= l && r <= R && mx[rt] <= val) {sum[rt] = 1LL * (r - l + 1) * val;mx[rt] = val;mn[rt] = val;lazy[rt] = val;return;}if(l == r) return;push(rt, l, r);int mid = l + r >> 1;update(L, R, val, lson);update(L, R, val, rson);pull(rt);}
} Tree;int main() {memset(maxPos, 0xc0, sizeof(maxPos));memset(minPos, 0x3f, sizeof(minPos));scanf("%d", &n);for(int i = 1; i <= n; i++) {scanf("%d", &a[i]);cnt[a[i]]++;if(i >= maxPos[a[i]][0]) {maxPos[a[i]][1] =  maxPos[a[i]][0];maxPos[a[i]][0] = i;} else if(i > maxPos[a[i]][1]) maxPos[a[i]][1] = i;if(i <= minPos[a[i]][0]) {minPos[a[i]][1] =  minPos[a[i]][0];minPos[a[i]][0] = i;} else if(i < minPos[a[i]][1]) minPos[a[i]][1] = i;}Tree.build(1, n, 1);for(int v = 200000; v >= 0; v--) {H[v] = 1LL * n * n - Tree.sum[1] + n;if(!v) break;mx[0] = mx[1] = -inf - 1;mn[0] = mn[1] = inf;tot = 0;for(int w = v; w <= 200000; w += v) {if(maxPos[w][0] == -inf - 1) continue;if(maxPos[w][0] >= mx[0]) mx[1] = mx[0], mx[0] = maxPos[w][0];else if(maxPos[w][0] > mx[1]) mx[1] = maxPos[w][0];if(maxPos[w][1] >= mx[0]) mx[1] = mx[0], mx[0] = maxPos[w][1];else if(maxPos[w][1] > mx[1]) mx[1] = maxPos[w][1];if(minPos[w][0] <= mn[0]) mn[1] = mn[0], mn[0] = minPos[w][0];else if(minPos[w][0] < mn[1]) mn[1] = minPos[w][0];if(minPos[w][1] <= mn[0]) mn[1] = mn[0], mn[0] = minPos[w][1];else if(minPos[w][1] < mn[1]) mn[1] = minPos[w][1];tot += cnt[w];}if(tot == 2) {int p1 = mn[0], p2 = mn[1];Tree.update(1, p1, p1, 1, n, 1);Tree.update(p1 + 1, p2, p2, 1, n, 1);Tree.update(p2 + 1, n, n + 1, 1, n, 1);} else if(tot == 3) {int p1 = mn[0], p2 = mn[1], p3 = mx[0];Tree.update(1, p1, p2, 1, n, 1);Tree.update(p1 + 1, p2, p3, 1, n, 1);Tree.update(p2 + 1, n, n + 1, 1, n, 1);} else if(tot > 3){int p1 = mn[0], p2 = mn[1], p3 = mx[1], p4 = mx[0];Tree.update(p2 + 1, n, n + 1, 1, n, 1);Tree.update(p1 + 1, p2, p4, 1, n, 1);Tree.update(1, p1, p3, 1, n, 1);}}LL ans = 0;for(int i = 1; i <= 200000; i++)ans += (H[i] - H[i - 1]) * i;printf("%lld\n", ans);return 0;
}/*
*/

转载于:https://www.cnblogs.com/CJLHY/p/10978952.html

Codeforces 671C Ultimate Weirdness of an Array 线段树 (看题解)相关推荐

  1. CodeForces 671C - Ultimate Weirdness of an Array

    题意: 给以一个定义, F(l, r) 的值表示序列 A[1:n]的子序列 A[1....(l-1),(r+1)...n] 之中 任意两个数的最大公约数的最大值. 求 Sum=∑i=1N∑j=1N(F ...

  2. Codeforces 1114F Please, another Queries on Array? 线段树

    Please, another Queries on Array? 利用欧拉函数的计算方法, 用线段树搞一搞就好啦. #include<bits/stdc++.h> #define LL ...

  3. Codeforces 213E Two Permutations 线段树 (看题解)

    Two Permutations 关键是没想到按大小顺序把第二个排列一个一个加入线段树, 然后线段树维护整体的hash值, 得到的hs值减去一个sub 之后与, 第一个排列的hash值比较. #inc ...

  4. Codeforces 311D Interval Cubing 数学 + 线段树 (看题解)

    Interval Cubing 这种数学题谁顶得住啊. 因为 (3 ^ 48) % (mod - 1)为 1 , 所以48个一个循环节, 用线段树直接维护. #include<bits/stdc ...

  5. CodeForces - 1401 F Reverse and Swap(线段树, 区间翻转, 区间交换,清晰易懂)

    CodeForces - 1401 F Reverse and Swap(线段树, 区间翻转, 区间交换)   首先一共有四个操作,第一个和第四个都是线段树的基本操作,直接用线段树实现.      第 ...

  6. 洛谷 P3373 【模板】线段树 2 题解

    洛谷 P3373 [模板]线段树 2 题解 题面 题目链接:[戳这里](https://www.luogu.org/problemnew/show/P3373) 题目描述 输入输出格式 输入输出样例 ...

  7. CodeForces - 1539F Strange Array(线段树区间合并)

    题目链接:点击查看 题目大意:给出一个长度为 nnn 的序列,规定位置 iii 的贡献是:设 x=a[i]x=a[i]x=a[i],选择一个包含 iii 的区间 [l,r][l,r][l,r],将其中 ...

  8. CF671C. Ultimate Weirdness of an Array

    n<=200000个<=200000的数问所有的f(i,j)的和,表示去掉区间i到j后的剩余的数字中任选两个数的最大gcd. 数论日常不会.. 先试着计算一个数组:Hi表示f(l,r)&l ...

  9. Educational Codeforces Round 37-F.SUM and REPLACE (线段树,线性筛,收敛函数)

    F. SUM and REPLACE time limit per test2 seconds memory limit per test256 megabytes inputstandard inp ...

最新文章

  1. 兼容ie8_兼容IE8的一些笔记
  2. 华为研发出了哪些芯片?
  3. 嵌入式中常见的几个死循环:
  4. matlab手写遗传算法解决一元函数最值问题(实例)
  5. 配置ODBC数据源——找不到SA账户的解决
  6. 云炬随笔20161224
  7. html选中radio隐藏div,radio类型的input标签选中后隐藏其他元素
  8. Transform数据权限浅析2之利用Java完成权限设置
  9. linux 重新分区 调整,重新调整磁盘分区linux / vm的大小
  10. rip v1不支持vlsm的原因
  11. Tree-AC训练实录
  12. android开发app初始化,Android 的 Application 初始化
  13. Perl获取前后任意月份月末
  14. 初识数据分析利器SPSS
  15. 基于PLC的智能化配电箱系统
  16. 总在说思科华为认证 可你真的清楚它们的区别吗?
  17. linux qt程序崩溃_Qt程式异常崩溃处理技巧(Win)
  18. 电脑开机是哪个键 电脑常用快捷键盘点
  19. 快速云:云计算供应商在合同谈判时可能拒绝的三个事项以及要求
  20. [教你做小游戏] 《五子棋》怎么存棋局信息?

热门文章

  1. 数学建模python教材推荐_数模竞赛专攻python应该准备什么?
  2. LeetCode 231. 2的幂
  3. 四步创建TCP客户端
  4. 云服务器可以文件服务器,云服务器 可以上传文件
  5. php中mysql_fetch_row_php中的mysql_fetch_row,mysql_fetch_array,mysql_fetch_object
  6. 26条C++的经典语录,哪几句戳中你的心!
  7. 图像锐化处理算法matlab,图像锐化matlab算法
  8. oc引导win方法_[OC更新]机械革命X1/X6TIS标压测试版更新
  9. 如何用耳机翻页_游戏耳机的经典之作—罗技(G)Astro A40体验
  10. java演出厅选票_高仿猫眼电影选座(选票)模块-b