bzoj2547: [Ctsc2002]玩具兵

划了一天水，其实我还是有点愧疚的。

传送门

其实是水题，然而我真是太蠢了。。。

首先不考虑天兵，其他兵要到一个点去一定是通过它-另一种兵-它……这样多次交换的，并且交换对象是无所谓的，和它换的兵最终会是在原位置。

而且骑兵和步兵的数量相等，就不存在一个兵找不到人跟它换的情况，那么一个确定的方案下换的次数最多的兵换的次数就是答案。

先spfa找出每个兵到每个位置的最小换的次数，然后二分答案，把每个兵向它可以到的地方连边，跑最大匹配，考虑到还有天兵，而天兵是可以一步登天的超神存在，所以如果最大匹配数+现在check的ans>=总兵数即OK。

因为一直在划水，断断续续地写这个代码 ，又困得要命，仿佛活在梦中，结果竟然1A了，神奇。

//Achen
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<queue>
#include<cmath>
const int N=;
typedef long long LL;
using namespace std;
int n,m,k,t,R[N],h[N][N],dis[N][N],dp[N][N][],vis[N][N][];
int dx[]={,,,-},dy[]={,-,,};

template<typename T>void read(T &x)  {
    char ch=getchar(); x=; T f=;
    while(ch!='-'&&(ch<''||ch>'')) ch=getchar();
    if(ch=='-') f=-,ch=getchar();
    for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;
}

struct node    {
    int x,y,o;
    node(){}
    node(int x,int y,int o):x(x),y(y),o(o){}
}b[N],p[N];

queue<node>que;
void spfa(int sx,int sy,int k) {
    memset(dp,,sizeof(dp));
    if(k) dp[sx][sy][]=,que.push(node(sx,sy,));
    else dp[sx][sy][]=,que.push(node(sx,sy,));
    while(!que.empty()) {
        node now=que.front();
        que.pop();
        int x=now.x,y=now.y,o=now.o;
        vis[x][y][o]=;
        for(int i=;i<;i++) {
            int xx=now.x+dx[i],yy=now.y+dy[i],oo,t;
            if(xx>=&&xx<=n&&yy>=&&yy<=m) {
                if((o==&&h[xx][yy]>=h[x][y])||(o==&&h[xx][yy]<=h[x][y])) t=dp[x][y][o],oo=o;
                else t=dp[x][y][o]+,oo=o^;
                if(dp[xx][yy][oo]>t) {
                    dp[xx][yy][oo]=t;
                    if(!vis[xx][yy][oo]) {
                        que.push(node(xx,yy,oo));
                        vis[xx][yy][oo]=;
                    }
                }
            }
        }
    }
}

int ok[N][N],vs[N],pr[N],sum[N];
int find(int x) {
    for(int i=;i<=*k;i++) if(ok[x][i]&&!vs[i]) {
        vs[i]=;
        if(!pr[i]||find(pr[i])) {
            pr[i]=x;
            return ;
        }
    }
    return ;
}

int ck(int ans) {
     int res=;
     for(int i=;i<=*k;i++)
         for(int j=;j<=t;j++) {
            if(dis[i][j]<=ans)
                for(int kk=;kk<=R[j];kk++)
                    ok[i][sum[j-]+kk]=;
            else for(int kk=;kk<=R[j];kk++)
                ok[i][sum[j-]+kk]=;
        }
    memset(pr,,sizeof(pr));
    for(int i=;i<=*k;i++) {
        memset(vs,,sizeof(vs));
        if(find(i)) res++;
    }
    if(res+ans>=*k) return ;
    else return ;
}

int main() {
    read(n); read(m); read(k); read(t);
    for(int i=;i<=*k+;i++) {
        read(b[i].x);
        read(b[i].y);
    }
    for(int i=;i<=t;i++) {
        read(p[i].x);
        read(p[i].y);
        read(R[i]); sum[i]=sum[i-]+R[i];
    }
    for(int i=;i<=n;i++) for(int j=;j<=m;j++) read(h[i][j]);
    memset(dis,,sizeof(dis));
    for(int i=;i<=*k;i++) {
        spfa(b[i].x,b[i].y,(i<=k));
        for(int j=;j<=t;j++)
            dis[i][j]=min(dp[p[j].x][p[j].y][],dp[p[j].x][p[j].y][]);
    }
    if(ck()) {
        puts(""); return ;
    }
    int l=,r=1e9,ans=r;
    while(l<=r) {
        int mid=((l+r)>>);
        if(ck(mid)) ans=mid,r=mid-;
        else l=mid+;
    }
    printf("%d\n",ans);
    return ;
}
/*
4 6 2 5
1 1   1 5   4 1   4 5    3 3
1 2 1   2 6 1   3 2 1   3 6 1   4 3 1
3 2 6 1 3 5
2 1 7 4 4 6
2 3 1 4 3 4
4 3 4 3 2 3
*/