BZOJ[2044]三维导弹拦截 匈牙利算法

传送门ber~

对于问题一，对x排序，$n^2$暴力就可以
对于问题二，将一个点对所有能接着走的点($x,y,z$均比自己大)连边，再跑最小路径覆盖就

最小路径覆盖可以用匈牙利算法

代码如下：

#include<algorithm>
#include<ctype.h>
#include<cstring>
#include<cstdio>
#define N 1050
using namespace std;
inline int read(){
    int x=0,f=1;char c;
    do c=getchar(),f=c=='-'?-1:f; while(!isdigit(c));
    do x=(x<<3)+(x<<1)+c-'0',c=getchar(); while(isdigit(c));
    return x*f;
}
struct Data{
    int x,y,z;
}a[N];
int n,ans,Timer,top;
int match[N],tim[N],fir[N],f[N];
struct Edge{
    int to,nex;
    Edge(int _=0,int __=0):to(_),nex(__){}
}nex[N*N];
inline void add(int x,int y){
    nex[++top]=Edge(y,fir[x]);
    fir[x]=top;
}
inline bool cmp(Data a,Data b){
    return a.x<b.x;
}
bool Hungary(int x){
    for(int i=fir[x];i;i=nex[i].nex)
        if(tim[nex[i].to]!=Timer){
            tim[nex[i].to]=Timer;
            if(!match[nex[i].to] || Hungary(match[nex[i].to])){
                match[nex[i].to]=x;
                return true;
            }
        }
    return false;
}
int main(){
    n=read();
    for(int i=1;i<=n;i++){
        a[i].x=read();a[i].y=read();a[i].z=read();
    }
    sort(a+1,a+n+1,cmp);
    for(int i=1;i<=n;i++){
        f[i]=1;
        for(int j=1;j<i;j++)
            if(a[j].x<a[i].x && a[j].y<a[i].y && a[j].z<a[i].z)
                f[i]=max(f[i],f[j]+1),add(j,i);
        ans=max(ans,f[i]);
    }
    printf("%d\n",ans);
    ans=0;
    for(int i=1;i<=n;i++){
        ++Timer;
        if(Hungary(i)) ans++;
    }
    printf("%d",n-ans);
return 0;
}