题目链接:点击这里
很早就开的一道题目,然而一直被卡爆了。
知道这样的点对有 nlogn对,但是没有转换成二维偏序问题。
而是用主席树保存了所有的历史状态,空间复杂度nlognlogn直接爆掉。
转换成二维偏序问题就很好做了。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<set>
#include<deque>
#include<map>
#include<vector>
#include<cmath>
#define ll long long
#define llu unsigned ll
using namespace std;
const double eps = 1e-8;
const ll lnf=0x3f3f3f3f3f3f3f3f;
const int inf = 0x3f3f3f3f;
const int maxn=100010;
const int mod=1e9+7;
struct node
{
int l,r,id;
node(){}
node(int a,int b,int c)
{
l=a,r=b,id=c;
}
}q[maxn*20];
int a[maxn],pos[maxn],ans[maxn];
int sum[maxn];
void add(int x,int val)
{
for(;x<maxn;x+=(x&-x))
sum[x]+=val;
}
int ask(int x)
{
int ans=0;
for(;x;x-=(x&-x))
ans+=sum[x];
return ans;
}
bool cmp(const node &a,const node &b)
{
if(a.l!=b.l) return a.l>b.l;
else return a.id<b.id;
}
int main(void)
{
int n,m;
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
scanf("%d",&a[i]),pos[a[i]]=i;
int cnt=0;
int l,r;
for(int i=1;i<=n;i++)
{
for(int j=2;j<=n/a[i];j++)
{
l=i;
r=pos[j*a[i]];
if(l>r) swap(l,r);
q[++cnt]=node(l,r,0);
}
}
for(int i=1;i<=m;i++)
{
cnt++;
scanf("%d%d",&q[cnt].l,&q[cnt].r);
q[cnt].id=i;
}
sort(q+1,q+cnt+1,cmp);
for(int i=1;i<=cnt;i++)
{
if(!q[i].id) add(q[i].r,1);
else ans[q[i].id]=ask(q[i].r);
}
for(int i=1;i<=m;i++)
printf("%d\n",ans[i]);
return 0;
}