首页 > 代码库 > POJ 3670 && POJ 3671 (dp)

POJ 3670 && POJ 3671 (dp)


最长不下降子序列的应用嘛。两题都是一样的。

POJ 3670:求给定序列按递增或递减排列时,所需改变的最小的数字的数目。

POJ 3671:求给定序列按递增排列时,所需改变的最小的数字的数目。

 

思路就是求最长不下降子序列,然后剩下的就是需要改变的字母。

最长不下降子序列:(我之前有写过,不懂请戳)http://blog.csdn.net/darwin_/article/details/38360997


POJ 3670:

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 30000+10
using namespace std;

int f[N],a[N];
int bin(int x,int r)   ////二分求比x大的第一个数字的位置
{
    int l=1,mid=(l+r)>>1,pos;
    while(l<=r)
    {
        mid=(l+r)>>1;
        if(a[mid]>x)
        {
            r=mid-1;
            pos=mid;
        }
        else
            l=mid+1;
    }
    return pos;
}
int main()
{
    int n;
    scanf("%d",&n);
    for(int i=0;i<n;i++)
        scanf("%d",&f[i]);
    int ca=0,cb=0;
    memset(a,0,sizeof(a));
    for(int i=0;i<n;i++)
    {
        if(i==0) a[++ca]=f[i];
        else if(f[i]>=a[ca]) a[++ca]=f[i];
        else
        {
            int pos=bin(f[i],ca);
            a[pos]=f[i];
        }
    }
    int ans=n-ca;
    reverse(f,f+n);  //翻转一次,就是求原序列的最长不上升子序列
    memset(a,0,sizeof(a));
    for(int i=0;i<n;i++)
    {
        if(i==0) a[++cb]=f[i];
        else if(f[i]>=a[cb]) a[++cb]=f[i];
        else
        {
            int pos=bin(f[i],cb);
            a[pos]=f[i];
        }
    }
    ans=min(ans,n-cb);  //求两者的最小值,即是答案。
    printf("%d\n",ans);
    return 0;
}

POJ 3671:

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 30000+10
using namespace std;

int f[N];
int bin(int x,int r)   //二分求比x大的第一个数字的位置
{
    int l=1,mid=(l+r)>>1,pos;
    while(l<=r)
    {
        mid=(l+r)>>1;
        if(f[mid]>x)
        {
            r=mid-1;
            pos=mid;
        }
        else
            l=mid+1;
    }
    return pos;
}
int main()
{
    int n;
    int cnt=0;
    scanf("%d",&n);
    for(int i=0;i<n;i++)
    {
        int t;
        scanf("%d",&t);
        if(i==0) f[++cnt]=t;
        else if(t>=f[cnt]) f[++cnt]=t;
        else
        {
            int pos=bin(t,cnt);
            f[pos]=t;
        }
    }
    printf("%d\n",n-(cnt));
    return 0;
}