题意定义f(l,r)为去掉[l,r]部分后剩下的数任意两个数的最大公约数的最大值

现在求f(l,r)的和

由于每个数ai最大只有200000,因此我们穷举因子x,记录以其为因子的a[i]的i值并按i升序。

下面我们从大到小穷举每个因子x,我们依次计算以f(l,r)=x的区间数,

有了上述的维护信息,我们很容易三种情况的区间满足f(l,r)=x(我就不具体细说了)

令pre[i]表示当前以第i个位置结尾的区间满足任意j(1<=j<=pre[i]),f(j,i)<x

显然f(l,r)=x的区间数=当前还没计算过的区间数-∑pre[i]

这就要维护∑pre[i],每次修改是将一段区间pre[i]>x的数赋值为x

显然pre[i]是单调不上升,因此可以用线段树维护之

 #include<bits/stdc++.h>

 using namespace std;
typedef long long ll;
int laz[*],mi[*],mx[*];
ll tr[*];
vector<int> g[];
int n,m;
void build(int i,int l,int r)
{
laz[i]=-;
if (l==r) mx[i]=mi[i]=tr[i]=l;
else {
int m=(l+r)>>;
build(i*,l,m);
build(i*+,m+,r);
tr[i]=tr[i*]+tr[i*+];
mx[i]=max(mx[i*],mx[i*+]);
mi[i]=min(mi[i*],mi[i*+]);
}
} void work(int i,int l,int r,int x,int y,int z)
{
// if (i==1) cout <<z<<endl;
if (x>y) return;
if (mx[i]<=z) return;
if (x<=l&&y>=r&&mi[i]>z)
{
mi[i]=mx[i]=laz[i]=z;
tr[i]=(ll)(r-l+)*z;
}
else {
int m=(l+r)>>;
if (laz[i]>-)
{
laz[i*]=laz[i*+]=laz[i];
mi[i*]=mx[i*]=laz[i];
mi[i*+]=mx[i*+]=laz[i];
tr[i*]=(ll)laz[i]*(m-l+);
tr[i*+]=(ll)laz[i]*(r-m);
laz[i]=-;
}
if (x<=m) work(i*,l,m,x,y,z);
if (y>m) work(i*+,m+,r,x,y,z);
tr[i]=tr[i*]+tr[i*+];
mx[i]=max(mx[i*],mx[i*+]);
mi[i]=min(mi[i*],mi[i*+]);
}
} int main()
{
scanf("%d",&n);
for (int i=; i<=n; i++)
{
int x;scanf("%d",&x);
for (int j=; j*j<=x; j++)
if (x%j==)
{
g[x/j].push_back(i);
if (j*j!=x) g[j].push_back(i);
}
m=max(m,x);
}
build(,,n);
ll pre=(ll)n*(n+)/;
ll ans=;
for (int i=m; i; i--)
{
if (g[i].size()<) continue;
int j=g[i].size()-;
work(,,n,g[i][]+,g[i][j]-,g[i][]);
work(,,n,,g[i][j-]-,);
work(,,n,g[i][]+,n,g[i][]);
ans+=(pre-tr[])*(ll)i;
// cout <<tr[1]<<endl;
pre=tr[];
}
printf("%lld\n",ans);
}

codeforces 671C的更多相关文章

  1. Codeforces 671C - Ultimate Weirdness of an Array(线段树维护+找性质)

    Codeforces 题目传送门 & 洛谷题目传送门 *2800 的 DS,不过还是被我自己想出来了 u1s1 这个 D1C 比某些 D1D 不知道难到什么地方去了 首先碰到这类问题我们肯定考 ...

  2. codeforces 671C Ultimate Weirdness of an Array 线段树+构造

    题解上说的很清楚了,我照着写的,表示膜拜题解 然后时间复杂度我觉得应该是O(nlogn),虽然常数略大,预处理和倒着扫,都是O(nlogn) #include <stdio.h> #inc ...

  3. 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 ...

  4. Ultimate Weirdness of an Array CodeForces - 671C (gcd,线段树)

    大意: 定义一个数列的特征值为两个数gcd的最大值, $f(l,r)$表示数列删除区间$[l,r]$的元素后剩余元素的特征值, 求$\sum_{i=1}^n\sum_{j=i}^n{f(i,j)}$ ...

  5. Codeforces 671C. Ultimate Weirdness of an Array(数论+线段树)

    看见$a_i\leq 200000$和gcd,就大概知道是要枚举gcd也就是答案了... 因为答案是max,可以发现我们很容易算出<=i的答案,但是很难求出单个i的答案,所以我们可以运用差分的思 ...

  6. python爬虫学习(5) —— 扒一下codeforces题面

    上一次我们拿学校的URP做了个小小的demo.... 其实我们还可以把每个学生的证件照爬下来做成一个证件照校花校草评比 另外也可以写一个物理实验自动选课... 但是出于多种原因,,还是绕开这些敏感话题 ...

  7. 【Codeforces 738D】Sea Battle(贪心)

    http://codeforces.com/contest/738/problem/D Galya is playing one-dimensional Sea Battle on a 1 × n g ...

  8. 【Codeforces 738C】Road to Cinema

    http://codeforces.com/contest/738/problem/C Vasya is currently at a car rental service, and he wants ...

  9. 【Codeforces 738A】Interview with Oleg

    http://codeforces.com/contest/738/problem/A Polycarp has interviewed Oleg and has written the interv ...

随机推荐

  1. window下对samba的清理操作

    windows清除访问samba局域网密码缓存 1.在dos窗口中输入control userpasswords2或者control keymgr.dll,然后[高级]/[密码管理],删掉保存的该机器 ...

  2. QT启动一个工程

    功能描述: 模拟如下页面. 当输入一个字符串时打开对应的应用程序. 实现方法: 1. 建立工程 2. 界面编辑: 3. 在test1.h中添加slot声明 4. test1.cpp中添加slot定义 ...

  3. ps学习笔记(二)

    1)选择所有图层: Ctrl+Alt+A2)查找层:ctrl+alt+shift+f,需要在层面板输入查找层名,可自动查找层:3)隔离层:可将选择图层,更改为隔离,只对选择的层编辑:注:图层面板中有一 ...

  4. Hibernate映射关系之_多对多

    多对多关系由于效率的原因,实际中会拆成相互的一对多的关系,不再累述

  5. 使用Kibana

    Kibana基本使用 https://www.elastic.co/guide/en/kibana/6.x/tutorial-load-dataset.html https://www.elastic ...

  6. Citrix Netscaler负载均衡算法

    Citrix Netscaler负载均衡算法 http://blog.51cto.com/caojin/1926308 众所周知,作为新一代应用交付产品的Citrix Netscaler具有业内领先的 ...

  7. BZOJ1103 [POI2007]大都市meg 【树剖】

    1103: [POI2007]大都市meg Time Limit: 10 Sec  Memory Limit: 162 MB Submit: 3038  Solved: 1593 [Submit][S ...

  8. POJ3450 Corporate Identity 【后缀数组】

    Corporate Identity Time Limit: 3000MS   Memory Limit: 65536K Total Submissions: 7662   Accepted: 264 ...

  9. 微信小程序使用原生WebSokcet实现断线重连及数据拼接

    以前做小程序为了应急找了个插件去链接WebSokcet,文章传送门. 回过头在新项目中再次使用时出现了些许问题,不一一赘述.遂决定好好用一下原生的WebSokcet. 一.说明 1.小程序原生的Web ...

  10. mysql共享表空间和独立表空间

    innodb这种引擎,与MYISAM引擎的区别很大.特别是它的数据存储格式等. 对于innodb的数据结构,首先要解决两个概念性的问题: 共享表空间以及独占表空间. 什么是共享表空间和独占表空间 共享 ...