noip模拟赛 柜（暴力）

分析：暴力的方法是非常显然的，从起点走一次，从终点走一次，路径相交的点即为所求，但是这样存图都很难存下，而且如果数据极端可能要求R*C次，时间空间都受不了.如果不需要记录整张图，并且一次能移动很多步就好了,标程用了树状数组+set,我不是很懂，只能打一个暴力了.

暴力：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;
const int maxn = ;
const int dx[] = { , , , - };
const int dy[] = { , -, ,  };
int n, m, A, B, ansx, ansy, cnt;
bool vis[maxn][maxn]; //512M的内存这数组应该是极限了.
int a[maxn][maxn];

struct node
{
    int x, y;
}e[maxn];

bool cmp(node a, node b)
{
    if (a.x == b.x)
        return a.y < b.y;
    return a.x < b.x;
}

void solve1()
{
    int x = , y = , dir = ;
    while ()
    {
        if (a[x][y] == )
            dir =  - dir;
        if (a[x][y] == )
            dir = (dir + ) % ;
        vis[x][y] = ;
        x = x + dx[dir], y = y + dy[dir];
        if (x > n || y > m || x <  || y < )
        {
            ansx = x;
            ansy = y;
            return;
        }
    }
}

void solve2()
{
    int x = n, y = m, dir = ;
    while ()
    {
        if (a[x][y] == )
            dir =  - dir;
        if (a[x][y] == )
            dir = (dir + ) % ;
        if (vis[x][y] && !a[x][y])
        {
            e[++cnt].x = x;
            e[cnt].y = y;
            vis[x][y] = ;
        }
        x += dx[dir], y += dy[dir];
        if (x > n || y > m || x <  || y < )
            return;
    }
}

int main()
{
    while (scanf("%d", &n) != EOF)
    {
        scanf("%d%d%d", &m, &A, &B);
        cnt = ;
        ansx = ;
        ansy = ;
        //memset(vis, 0, sizeof(vis));
        for (int i = ; i <= n; i++)
            for (int j = ; j <= m; j++)
                vis[i][j] = a[i][j] = ;
        //memset(a, 0, sizeof(a));
        for (int i = ; i <= A; i++)
        {
            int x, y;
            scanf("%d%d", &x, &y);
            a[x][y] = ;
        }
        for (int i = ; i <= B; i++)
        {
            int x, y;
            scanf("%d%d", &x, &y);
            a[x][y] = ;
        }
        solve1();
        if (ansx == n && ansy == m + )
        {
            printf("0\n");
            continue;
        }
        solve2();
        if (cnt == )
        {
            printf("-1\n");
            continue;
        }
        sort(e + , e +  + cnt, cmp);
        printf("%d %d %d\n", cnt, e[].x, e[].y);
    }

    return ;
}

std:

#include<cstdio>
#include<set>
#include<vector>
#include<algorithm>
using namespace std;
#define rep(i,n) for (int i=0;i<n;++i)
#define pb push_back
#define mk make_pair
#define X first
#define Y second
#define tree int t,int l,int r
#define left t*2,l,mid
#define right t*2+1,mid+1,r
#define M int mid=l+r>>1
const int N=;
typedef pair<int,int> pr;
typedef vector<pair<int,pr> > seq;
set<pr> a[N],b[N]; seq f1,g1,f2,g2; int Case,n,m,R,C,x,y,ll,rr,c[N]; long long ans;
int get(int x){int res=; for (;x;x-=x&-x) res+=c[x]; return res;}
void add(int x,int v){for (;x<=C;x+=x&-x) c[x]+=v;}
void ins(int side){scanf("%d%d",&x,&y),a[x].insert(mk(y,side)),b[y].insert(mk(x,side));}
bool track(int x,int y,int d,seq &f,seq &g)
{
    f.clear(),g.clear(); set<pr> :: iterator it;
    for (;;){
        if (d&){
            if (d==){
                it=b[y].upper_bound(mk(x,)); f.pb(mk(x+,mk(y,)));
                if (it==b[y].end()) return f.pb(mk(R+,mk(y,-))),;
                f.pb(mk(it->X,mk(y,-))),x=it->X,d=it->Y?:;
            }else{
                it=b[y].lower_bound(mk(x,)); f.pb(mk(x,mk(y,-)));
                if (it==b[y].begin()) return f.pb(mk(,mk(y,))),; --it;
                f.pb(mk(it->X+,mk(y,))),x=it->X,d=it->Y?:;
            }
        }else{
            if (d==){
                it=a[x].upper_bound(mk(y,));
                if (it==a[x].end()) return g.pb(mk(x,mk(y+,C))),x==R;
                g.pb(mk(x,mk(y+,it->X-))),y=it->X,d=it->Y?:;
            }else{
                it=a[x].lower_bound(mk(y,));
                if (it==a[x].begin()) return g.pb(mk(x,mk(,y-))),; --it;
                g.pb(mk(x,mk(it->X+,y-))),y=it->X,d=it->Y?:;
            }
        }
    }
}
void work(seq &f,seq &g)
{
    sort(f.begin(),f.end()),sort(g.begin(),g.end());
    int m=f.size(),n=g.size(),j=;
    rep(i,n){
        while (j<m && f[j].X<=g[i].X) add(f[j].Y.X,f[j].Y.Y),++j;
        ll=g[i].Y.X,rr=g[i].Y.Y; int res=get(rr)-get(ll-); ans+=res;
        if (g[i].X<x && res){
            x=g[i].X,y=ll;
            for (int j=;j>=;--j)
                if (y+(<<j)<=rr && !(get(y-+(<<j))-get(y-))) y+=<<j;
            }
        }
    while (j<m) add(f[j].Y.X,f[j].Y.Y),++j;
}
int main()
{
    freopen("safe.in", "r", stdin);
    freopen("safe.out", "w", stdout);

    while (scanf("%d%d%d%d",&R,&C,&n,&m)!=EOF){
        rep(i,R+) a[i].clear(); rep(j,C+) b[j].clear();
        rep(i,n+m) ins(i<n);
        if (track(,,,f1,g1)){puts(""); continue;} track(R,C+,,f2,g2);
        ans=,x=R+,work(f1,g2),work(f2,g1);
        if (ans) printf("%lld %d %d\n",ans,x,y); else puts("-1");
    }
    return ;
}