POJ 3179 离散化+二维前缀和+枚举(二分？)

离散化和前缀和以前做过，但是不熟，所以借鉴的lyd的代码(不过好像他也没用二分查找,虽然书上这么写的）不过代码中有一些剪枝和为下一步预处理的的操作可能优化了时间，反正62ms过了。。。

附上代码:

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<set>
#include<vector>
#define INF 1<<30
using namespace std;
vector<int> X,Y;
const int maxn=510;
//typedef point pair<int,int>;
int sum[maxn][maxn];
vector<pair<int,int> >p; 
int main(){
    int C,n,x,y;
    int a,b;
    scanf("%d%d",&C,&n);
    for(int i=1;i<=n;i++){
        scanf("%d%d",&a,&b);
        p.push_back(make_pair(a,b));
        X.push_back(a);
        Y.push_back(b); 
    }
    //排序 
    sort(p.begin(),p.end());
    sort(X.begin(),X.end());
    X.erase(unique(X.begin(),X.end()),X.end());//去重 
    sort(Y.begin(),Y.end()); 
    Y.erase(unique(X.begin(),X.end()),X.end());
    x=X.size(),y=Y.size();
    X.push_back(INF),Y.push_back(INF);
    //求二维前缀和 
    //sum[i][j]是点X[i-1][j-1]的前缀和 
    int pos=0;
    for(int i=1;i<=x;i++){
        for(int j=1;j<=y;j++){
            sum[i][j]=sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1];
            while(pos<n&&p[pos].first==X[i-1]&&p[pos].second==Y[j-1]){//判断该点是否有草 
                sum[i][j]++; 
                pos++;
            }
        }
    }
    //查找 
    int ans=INF;
    for(int i=0;i<x;i++){//枚举起点x坐标 
        int i1=i+1;//终点x坐标 
        int s=0;//边长 
        int j=0;//起点y坐标 
        int j1=1;//终点y坐标 
        int val=sum[i1][j1]-sum[i][j1]-sum[i1][j]+sum[i][j];
        while(1){
            while(val<C&&(i1<x||j1<y)){
                s=min(X[i1]-X[i],Y[j1]-Y[j]);
                while(X[i1]-X[i]<=s)i1++;//因为sum[i][j]是坐标(x-1,y-1)的前缀和 
                while(Y[j1]-Y[j]<=s)j1++;
                val=sum[i1][j1]-sum[i][j1]-sum[i1][j]+sum[i][j];
            }
            if(val<C)break;//剪枝 
            ans=min(ans,s+1); 
            j++;//y坐标+1 
            if(j==y)break;//y坐标枚举到头了,退出循环 
            if(j1<=j){ 
                j1=j+1;
                s=0;
            }
            else{
                s-=Y[j]-Y[j-1];//边长也相应的减少 
            }
            while(X[i1-1]-X[i]>s) i1--;//为下一轮计算做准备
            val=sum[i1][j1]-sum[i][j1]-sum[i1][j]+sum[i][j];
        }
    }
    printf("%d\n",ans);
}