poj 3141 Distant Galaxy

这道题其实处理起复杂度有点高，我们采用的是枚举所有的矩形，判断落在边上的点有多少个，这道题我们看到坐标会很大，所以离散化是必不可少的，离散化之后，地图就会被划为100*100大小的，就可以方便操作了。首先要进行一次初始化，co[k][j][i][0]中存的是第K行的第j列到第i列中的点的个数，同理，co[k][j][i][1]中存的是第列中第j行到第i行的线段上点的个数。

之后我采用的是四层循环枚举，复杂度远远高于大神们的三次方，在此继续膜拜大神。

/*#########################################################################
# File Name: d.cpp
# Author: CaoLei
# Created Time: 2015/7/15 13:52:11
#########################################################################*/

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <set>
#include <queue>
#include <map>
using namespace std;
#define MAX(x,y) (((x)>(y))?(x):(y))
#define MIN(x,y) (((x)<(y))?(x):(y))
#define N 500010
#define pi acos(-1.0)
#define inf 100000000
typedef long long ll;
typedef unsigned long long ull;
struct node {
    int x1,y1;
    int x,y;
}a[105];

int cmp( const void *a ,const void *b){
return (*(node *)a).x1 > (*(node *)b).x1 ? 1 : -1;
}

int cmp1( const void *a ,const void *b){
return (*(node *)a).y1 > (*(node *)b).y1 ? 1 : -1;
}
int m[105][105];
int s[105];
int co[105][105][105][2];
int main(){
    int ca=0;
    int n;
    //freopen("in.txt","r",stdin);
    while(~scanf("%d",&n)){
        if(n==0) break;
        printf("Case %d: ",++ca);
        if(n==1){ 
            int aa,aaa;
            scanf("%d%d",&aa,&aaa);
            printf("1\n");
            continue;
        }
        memset(a,0,sizeof(a));
        memset(co,0,sizeof(co));
        memset(m,0,sizeof(m));
        for(int i=1;i<=n;i++){
            scanf("%d%d",&a[i].x1,&a[i].y1);
        }
        qsort(a+1,n,sizeof(a[0]),cmp);
        int ex=1;
        a[1].x=1;
        for(int i=2;i<=n;i++){
            if(a[i].x1==a[i-1].x1) {
                a[i].x=ex;
            }
            else{
                a[i].x=++ex;
            }
        }
        qsort(a+1,n,sizeof(a[0]),cmp1);
        a[1].y=1;
        int ey=1;
        for(int i=2;i<=n;i++){
            if(a[i].y1==a[i-1].y1) {
                a[i].y=ey;
            }
            else{
                a[i].y=++ey;
            }
            m[a[i].x][a[i].y]=1;
        }
        m[a[1].x][a[1].y]=1;
        for(int k=1;k<=ex;k++){
            memset(s,0,sizeof(s));
            for(int i=1;i<=ey;i++)
              if(m[k][i]) s[i]=s[i-1]+1;
              else s[i]=s[i-1];
            for(int i=1;i<=ey;i++){
                for(int j=1;j<i;j++){
                    co[k][j][i][0]=s[i]-s[j];
                    if(m[k][j]) co[k][j][i][0]++;
                }
            }
        }
        for(int k=1;k<=ey;k++){
            memset(s,0,sizeof(s));
            for(int i=1;i<=ex;i++)
              if(m[i][k]) s[i]=s[i-1]+1;
              else s[i]=s[i-1];
            for(int i=1;i<=ex;i++){
                for(int j=1;j<i;j++){
                    co[k][j][i][1]=s[i]-s[j];
                    if(m[j][k]) co[k][j][i][1]++;
                }
            }
        }
        if(ex==1){
            printf("%d\n",co[1][1][ey][0]);
            continue;
        }
        if(ey==1){
            printf("%d\n",co[1][1][ex][1]);
            continue;
        }
        int max=0;
        for(int x1=1;x1<=ex;x1++){
            for(int y1=1;y1<=ey;y1++){
                for(int x2=x1+1;x2<=ex;x2++){
                    for(int y2=y1+1;y2<=ey;y2++){
                        int ans=co[x1][y1][y2][0]+co[x2][y1][y2][0]+co[y1][x1][x2][1]+co[y2][x1][x2][1];
                        if(m[x1][y1]) ans--;
                        if(m[x1][y2]) ans--;
                        if(m[x2][y1]) ans--;
                        if(m[x2][y2]) ans--;
                        if(ans>max){
                            max=ans;
                        }
                    }
                }
            }
        }
        printf("%d\n",max);
    }
    return 0;
}