SSL 1082 2317 拦截导弹#匈牙利算法#

题目

某国为了防御敌国的导弹袭击，发展出一种导弹拦截系统，虽然它的第一发炮弹能够到达任意的高度，但是以后每一发炮弹都不能高于前一发的高度。计算这套系统最多能拦截多少导弹，如果要拦截所有导弹最少要配备多少套这种导弹拦截系统。

---

分析

对于第一个问题，用动态规划
状态转移方程：
$f[i]=max(f[i],f[j]+1)(i>j且a[i]>=a[j])$
对于第二个问题，虽然也可以用动态规划，但是因为数据太水，贪心也能过，为了介绍最大匹配，所以这样贴代码

---

代码

#include <cstdio>
#include <cstring>
using namespace std;
struct node{int y,next;}e[10001]; int s,ans,cover[501],a[501],n=1,m,ls[501],link[501];
int f[501],max;
void add(int i,int j){e[++m].y=j; e[m].next=ls[i]; ls[i]=m;}
bool find(int x){
	int t=ls[x];
	while (t){
		if (!cover[e[t].y]){
			cover[e[t].y]=1;
			int q=link[e[t].y];
			link[e[t].y]=x;
			if (!q||find(q)) return 1;
			link[e[t].y]=q;
		}
		t=e[t].next;
	}
	return 0;
}
int main(){
	while (scanf("%d",&a[n])==1) n++; n--; 
	f[n]=1;
	for (int i=n-1;i>=1;i--){//动态规划
	s=1;
	for (int j=i+1;j<=n;j++)
	if (a[i]>a[j]&&f[j]+1>s)
	s=f[j]+1,add(i,j);
	else if (a[i]>=a[j]) add(i,j);//构图
	f[i]=s;
	if (max<f[i]) max=f[i];
	}
	printf("%d\n",max);
	for (int i=1;i<=n;i++)
	ans+=find(i),memset(cover,0,sizeof(cover));
	printf("%d",n-ans);
    return 0;
}