poj1548 Robots(最小路径覆盖/LIS)

这道题很显然是个最小路径覆盖问题。可以用二分图匹配来解决。
其实也可以用dp来解决。把所有点按第一关键字x,第二关键字y从小到大排序。则从1~n的点已经满足了第一维,从i->j(i< j)只需要满足p[i].y<=p[j].y即可。例如样例就是4 4 6 4 7.根据题意,这个偏序集的链是一个不降序列。我们现在就是要求这个偏序集的最小链划分数(也就是图中的最小链覆盖数),根据Dilworth定理,也就是求最长反链的长度。也就是求最长下降子数列的长度。可以dp解决。

最小路径覆盖

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <stack>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 600
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,h[N],num=1,ans=0,bf[N];
bool f[N];
struct edge{
    int to,next;
}data[200000];
struct point{
    int x,y;
}p[N];
inline bool cmp(point a,point b){return a.x==b.x?a.y<b.y:a.x<b.x;}
inline void add(int x,int y){
    data[++num].to=y;data[num].next=h[x];h[x]=num;
}
inline bool find(int x){
    for(int i=h[x];i;i=data[i].next){
        int y=data[i].to;if(f[y]) continue;f[y]=1;
        if(!bf[y]||find(bf[y])){bf[y]=x;return 1;}
    }return 0;
}
int main(){
//  freopen("a.in","r",stdin);
    while(1){
        memset(h,0,sizeof(h));num=1;n=0;
        memset(bf,0,sizeof(bf));ans=0;
        while(1){
            p[++n].x=read();p[n].y=read();
            if(p[n].x==-1) return 0;
            if(p[n].x==0) break;
        }--n;sort(p+1,p+n+1,cmp);
        for(int i=1;i<=n;++i)
            for(int j=i+1;j<=n;++j)
                if(p[j].y>=p[i].y) add(i,j);
        for(int i=1;i<=n;++i){
            memset(f,0,sizeof(f));if(find(i)) ans++;
        }printf("%d\n",n-ans);
    }
}

dp

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <stack>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 600
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,ans=0,dp[N];
struct point{
    int x,y;
}p[N];
inline bool cmp(point a,point b){return a.x==b.x?a.y<b.y:a.x<b.x;}
int main(){
//  freopen("a.in","r",stdin);
    while(1){
        n=0;memset(dp,0,sizeof(dp));ans=0;
        while(1){
            p[++n].x=read();p[n].y=read();
            if(p[n].x==-1) return 0;
            if(p[n].x==0) break;
        }--n;sort(p+1,p+n+1,cmp);
        for(int i=1;i<=n;++i)
            for(int j=1;j<i;++j)
                if(p[j].y>p[i].y) dp[i]=max(dp[i],dp[j]+1);
        for(int i=1;i<=n;++i) ans=max(ans,dp[i]);
        printf("%d\n",ans+1);
    }
}

手写二分…

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <stack>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 600
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,ans=0,d[N];
struct point{
    int x,y;
}p[N];
inline bool cmp(point a,point b){return a.x==b.x?a.y<b.y:a.x<b.x;}
inline int find(int l,int r,int x){//严格下降的数组,找第一个<=x的数 
    while(l<=r){
        int mid=l+r>>1;
        if(d[mid]<=x) r=mid-1;
        else l=mid+1;
    }return r+1;
}
int main(){
//  freopen("a.in","r",stdin);
    while(1){
        n=0;memset(d,0,sizeof(d));ans=0;
        while(1){
            p[++n].x=read();p[n].y=read();
            if(p[n].x==-1) return 0;
            if(p[n].x==0) break;
        }--n;sort(p+1,p+n+1,cmp);d[++ans]=p[1].y;
        for(int i=2;i<=n;++i){
            int x=p[i].y;
            if(x<d[ans]) d[++ans]=x;
            else{
                int id=find(1,ans,x);
                d[id]=x;
            }
        }printf("%d\n",ans);
    }
}
发布了1239 篇原创文章 · 获赞 36 · 访问量 28万+
展开阅读全文

没有更多推荐了,返回首页

©️2019 CSDN 皮肤主题: 编程工作室 设计师: CSDN官方博客

分享到微信朋友圈

×

扫一扫,手机浏览