CF301D.Yaroslav and Divisors(扫描线+树状数组)

简介: CF301D.Yaroslav and Divisors(扫描线+树状数组)

Yaroslav and Divisors

原题链接

思路:

首先,朴素的暴力做法很好想,对于每次询问的区间都O ( n 2 )的去遍历,计算数量;考虑怎么优化。

由于询问的点对是由两个数组成的,所以可以先离线,对右端点排序,进行降维操作,然后可以枚举右端点,计算前面的贡献。

每次枚举到某个右端点x时,先将x前面的贡献加上,再对每个右端点为x的询问进行查询。

但是,对于每一个x,如果每次都暴力找前面和他成因子的数,时间复杂度又会退化成n 2,可以先预处理出每个数的因子位置,每次修改的时候直接修改这些位置。

要用到的是单点修改+区间查询,树状数组即可维护。

具体做法是:

1.对于给定的序列的每个数都统计比这个数的下标小的因子的位置;

2.将询问离线下来,对于每个右端点,存储询问的左端点跟询问下标;存储下标是因为这里进行了离线操作。

3.从小到大枚举右端点,对于每一个右端点,先将下标比他小的因子个数加上,即将每个因子的位置都+ 1,表示该右端点与该因子对答案的贡献是1;再询问[ l , r ]的答案,即q u e r y ( r ) − q u e r y ( l − 1 );

代码:

// Problem: D. Yaroslav and Divisors
// Contest: Codeforces - Codeforces Round #182 (Div. 1)
// URL: http://codeforces.com/problemset/problem/301/D
// Memory Limit: 256 MB
// Time Limit: 2000 ms
// Author:Cutele
// 
// Powered by CP Editor (https://cpeditor.org)
#pragma GCC optimize(2)
#pragma GCC optimize(3,"Ofast","inline")
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll, ll>PLL;
typedef pair<int, int>PII;
typedef pair<double, double>PDD;
#define I_int ll
inline ll read(){ll x = 0, f = 1;char ch = getchar();while(ch < '0' || ch > '9'){if(ch == '-')f = -1;ch = getchar();}while(ch >= '0' && ch <= '9'){x = x * 10 + ch - '0';ch = getchar();}return x * f;}
inline void write(ll x){if (x < 0) x = ~x + 1, putchar('-');if (x > 9) write(x / 10);putchar(x % 10 + '0');}
#define read read()
#define closeSync ios::sync_with_stdio(0);cin.tie(0);cout.tie(0)
#define multiCase int T;cin>>T;for(int t=1;t<=T;t++)
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define repp(i,a,b) for(int i=(a);i<(b);i++)
#define per(i,a,b) for(int i=(a);i>=(b);i--)
#define perr(i,a,b) for(int i=(a);i>(b);i--)
ll ksm(ll a, ll b,ll mod){ll res = 1;while(b){if(b&1)res=res*a%mod;a=a*a%mod;b>>=1;}return res;}
const int maxn=2e5+7;
int a[maxn],id[maxn],n,m,ans[maxn],tr[maxn];
vector<int>g[maxn];
vector<PII>q[maxn];
int lowbit(int x){
  return x&-x;
}
void update(int pos,int val){
  while(pos<=n){
    tr[pos]+=val;
    pos+=lowbit(pos);
  }
}
int query(int pos){
  int res=0;
  while(pos){
    res+=tr[pos];
    pos-=lowbit(pos);
  }
  return res;
}
int main(){
  n=read,m=read;
  rep(i,1,n) a[i]=read,id[a[i]]=i;
  for(int i=1;i<=n;i++){
    for(int j=i;j<=n;j+=i){
      int u=id[i],v=id[j];
      if(u<v) swap(u,v);
      g[u].push_back(v);
    }
  }
  rep(i,1,m){
    int l=read,r=read;
    q[r].push_back({l,i});
  }
  rep(i,1,n){
    for(int j=0;j<g[i].size();j++){
      int t=g[i][j];
      update(t,1);
    }
    for(int j=0;j<q[i].size();j++){
      PII t=q[i][j];
      int l=t.first,now=t.second;
      ans[now]=query(i)-query(l-1);
    }
  }
  rep(i,1,m) printf("%d\n",ans[i]);
  return 0;
}
目录
相关文章
|
7月前
|
人工智能 移动开发 算法
【动态规划】【C++算法】LeetCoce996正方形数组的数目
【动态规划】【C++算法】LeetCoce996正方形数组的数目
|
算法
【算法专题突破】双指针 - 有效三角形的个数(5)
【算法专题突破】双指针 - 有效三角形的个数(5)
35 0
|
7月前
|
算法 C++
(C++)有效三角形的个数--双指针法
(C++)有效三角形的个数--双指针法
70 0
|
算法 C++
剑指offer(C++)-JZ29:顺时针打印矩阵(算法-模拟)
剑指offer(C++)-JZ29:顺时针打印矩阵(算法-模拟)
|
算法 C++
剑指offer(C++)-JZ40:最小的K个数(算法-排序)
剑指offer(C++)-JZ40:最小的K个数(算法-排序)
【LeetCode】BM1 反转链表、NC21 链表内指定区间反转
BM1 反转链表 描述: 给定一个单链表的头结点pHead(该头节点是有值的,比如在下图,它的val是1),长度为n,反转该链表后,返回新链表的表头。
54 0
|
Go vr&ar
CF中的线段树题目
CF中的线段树题目
84 0
剑指offer_数组---顺时针打印矩阵
剑指offer_数组---顺时针打印矩阵
53 0
|
算法 程序员
【牛客算法BM2】 链表内指定区间反转
你好,欢迎来到我的博客!作为一名程序员,我经常刷LeetCode题目来提升自己的编程能力。在我的博客里,我会分享一些我自己做过的题目和解题思路,希望能够帮助到大家。今天,我想和大家分享一道挑战性较高的题目,让我们一起来挑战一下吧!作者也是在学习的过程中刷到有意思的题目就会选择与大家分享,并且提供较优解,关于力扣 的文章全部放在博客,如果大家喜欢记得支持作者。🤓
AC牛客 BM3链表中的节点每k个一组翻转
AC牛客 BM3链表中的节点每k个一组翻转
86 0
AC牛客 BM3链表中的节点每k个一组翻转