The Preliminary Contest for ICPC Asia Xuzhou 2019:I. query(二维偏序):

题目链接:点击这里

很早就开的一道题目,然而一直被卡爆了。
知道这样的点对有 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;
}

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值