Codeforces 671C Ultimate Weirdness of an Array 线段树 (看题解)
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 线段树 (看题解)相关推荐
- 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 ...
- Codeforces 1114F Please, another Queries on Array? 线段树
Please, another Queries on Array? 利用欧拉函数的计算方法, 用线段树搞一搞就好啦. #include<bits/stdc++.h> #define LL ...
- Codeforces 213E Two Permutations 线段树 (看题解)
Two Permutations 关键是没想到按大小顺序把第二个排列一个一个加入线段树, 然后线段树维护整体的hash值, 得到的hs值减去一个sub 之后与, 第一个排列的hash值比较. #inc ...
- Codeforces 311D Interval Cubing 数学 + 线段树 (看题解)
Interval Cubing 这种数学题谁顶得住啊. 因为 (3 ^ 48) % (mod - 1)为 1 , 所以48个一个循环节, 用线段树直接维护. #include<bits/stdc ...
- CodeForces - 1401 F Reverse and Swap(线段树, 区间翻转, 区间交换,清晰易懂)
CodeForces - 1401 F Reverse and Swap(线段树, 区间翻转, 区间交换) 首先一共有四个操作,第一个和第四个都是线段树的基本操作,直接用线段树实现. 第 ...
- 洛谷 P3373 【模板】线段树 2 题解
洛谷 P3373 [模板]线段树 2 题解 题面 题目链接:[戳这里](https://www.luogu.org/problemnew/show/P3373) 题目描述 输入输出格式 输入输出样例 ...
- CodeForces - 1539F Strange Array(线段树区间合并)
题目链接:点击查看 题目大意:给出一个长度为 nnn 的序列,规定位置 iii 的贡献是:设 x=a[i]x=a[i]x=a[i],选择一个包含 iii 的区间 [l,r][l,r][l,r],将其中 ...
- CF671C. Ultimate Weirdness of an Array
n<=200000个<=200000的数问所有的f(i,j)的和,表示去掉区间i到j后的剩余的数字中任选两个数的最大gcd. 数论日常不会.. 先试着计算一个数组:Hi表示f(l,r)&l ...
- Educational Codeforces Round 37-F.SUM and REPLACE (线段树,线性筛,收敛函数)
F. SUM and REPLACE time limit per test2 seconds memory limit per test256 megabytes inputstandard inp ...
最新文章
- 兼容ie8_兼容IE8的一些笔记
- 华为研发出了哪些芯片?
- 嵌入式中常见的几个死循环:
- matlab手写遗传算法解决一元函数最值问题(实例)
- 配置ODBC数据源——找不到SA账户的解决
- 云炬随笔20161224
- html选中radio隐藏div,radio类型的input标签选中后隐藏其他元素
- Transform数据权限浅析2之利用Java完成权限设置
- linux 重新分区 调整,重新调整磁盘分区linux / vm的大小
- rip v1不支持vlsm的原因
- Tree-AC训练实录
- android开发app初始化,Android 的 Application 初始化
- Perl获取前后任意月份月末
- 初识数据分析利器SPSS
- 基于PLC的智能化配电箱系统
- 总在说思科华为认证 可你真的清楚它们的区别吗?
- linux qt程序崩溃_Qt程式异常崩溃处理技巧(Win)
- 电脑开机是哪个键 电脑常用快捷键盘点
- 快速云:云计算供应商在合同谈判时可能拒绝的三个事项以及要求
- [教你做小游戏] 《五子棋》怎么存棋局信息?
热门文章
- 数学建模python教材推荐_数模竞赛专攻python应该准备什么?
- LeetCode 231. 2的幂
- 四步创建TCP客户端
- 云服务器可以文件服务器,云服务器 可以上传文件
- php中mysql_fetch_row_php中的mysql_fetch_row,mysql_fetch_array,mysql_fetch_object
- 26条C++的经典语录,哪几句戳中你的心!
- 图像锐化处理算法matlab,图像锐化matlab算法
- oc引导win方法_[OC更新]机械革命X1/X6TIS标压测试版更新
- 如何用耳机翻页_游戏耳机的经典之作—罗技(G)Astro A40体验
- java演出厅选票_高仿猫眼电影选座(选票)模块-b