Codeforces 685D

1.枚举x的取值范围[l,l+k-1],把所有位于这个范围的点存起来;
2.枚举所有的y取值范围[l,l+k-1],然后计数更新答案
3.复杂度仔细想想,还真是对的

#include<cstdio>
#include<iostream>
#include<vector>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1e5+100;
typedef __int64 LL;
#define pb push_back

int sta[N*2],top;
vector<int> vv[N*2];
int yy[N*2];
int tt[N*2],cnt;

LL ans[N];
void fun(int co,int k3,int n){
    int i=0,j=0,l=-2e9,r=l+k3-1;
    int sum=0;
    while(1){
        if(j>=n||sta[tt[i]]+1-l<sta[tt[j]]-r){
            sum-=yy[tt[i]];
            l=sta[tt[i]]+1;
            r=l+k3-1;
            i++;
        }
        else if(sta[tt[i]]+1-l>sta[tt[j]]-r){
            sum+=yy[tt[j]];
            r=sta[tt[j]];
            l=r-k3+1;
            j++;
        }
        else {
            sum-=yy[tt[i]];
            sum+=yy[tt[j]];
            r=sta[tt[j]];
            l=r-k3+1;
            i++;j++;
        }
        //printf("i%d j%d l%d r%d sum%d\n",i,j,l,r,sum);
        if(i>=n&&j>=n)break;
        int coo;
        if(j>=n||sta[tt[i]]+1-l<sta[tt[j]]-r){
            coo=sta[tt[i]]+1-l;
        }
        else {
            coo=sta[tt[j]]-r;
        }
        //printf("%dcoo\n",coo);
        ans[sum]+=(LL)co*coo;
    }
}

int aa[N],bb[N];
int main(){
    #ifdef DouBi
    freopen("in.cpp","r",stdin);
    #endif // DouBi
    int n,k3;
    while(scanf("%d%d",&n,&k3)!=EOF){
        for(int i=0;i<n;i++){
            scanf("%d%d",&aa[i],&bb[i]);
            sta[top++]=aa[i];sta[top++]=bb[i];
        }
        sort(sta,sta+top);top=unique(sta,sta+top)-sta;
        for(int i=0;i<top;i++){
            vv[i].clear();
        }
        for(int i=0;i<n;i++){
            aa[i]=lower_bound(sta,sta+top,aa[i])-sta;
            bb[i]=lower_bound(sta,sta+top,bb[i])-sta;
            vv[aa[i]].pb(bb[i]);
        }

        memset(yy,0,sizeof(yy));
        memset(ans,0,sizeof(ans));

        int i=0,j=0,l=-2e9,r=l+k3-1;
        //for(int j=0;j<top;j++)printf("%d ",sta[j]);printf("\n");
        while(1){
            if(j>=top||sta[i]+1-l<sta[j]-r){
                l=sta[i]+1;
                r=l+k3-1;
                i++;
            }
            else if(sta[i]+1-l>sta[j]-r){
                r=sta[j];
                l=r-k3+1;
                j++;
            }
            else {
                r=sta[j];
                l=r-k3+1;
                i++;j++;
            }
            //printf("i%d j%d l%d r%d\n",i,j,l,r);
            cnt=0;
            if(i>=top&&j>=top)break;
            for(int k=i;k<j;k++){
                for(int k1=0;k1<vv[k].size();k1++){
                    tt[cnt++]=vv[k][k1];
                    yy[vv[k][k1]]=0;
                }
            }
            //for(int k=0;k<cnt;k++)printf("%d ",tt[k]);printf("\n");
            for(int k=0;k<cnt;k++){
                yy[tt[k]]++;
            }
            sort(tt,tt+cnt);cnt=unique(tt,tt+cnt)-tt;
            //for(int k=0;k<cnt;k++)printf("%d ",sta[tt[k]]);printf("\n");
            int co;
            if(j>=top||sta[i]+1-l<sta[j]-r){
                co=sta[i]+1-l;
            }
            else {
                co=sta[j]-r;
            }
            //printf("co%d\n",co);
            fun(co,k3,cnt);
            //printf("..........\n");
        }
        for(int i=1;i<=n;i++){
            printf("%I64d",ans[i]);
            if(i<n)printf(" ");
            else printf("\n");
        }
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值