线性DP

LCIS

单调优化思想

    a[0]=b[0]=-9999;
    FOR(i,1,n)
    {
        int val=0;
        if(b[0]<a[i]) val=f[i-1][0];
        FOR(j,1,n) 
        {
            if(a[i]==b[j]) f[i][j]=val+1;
            else f[i][j]=f[i-1][j];
            if(b[j]<a[i]) val=max(val,f[i-1][j]);
        }
    }
    FOR(i,1,n) ans=max(f[n][i],ans);
    printf("%d",ans);

CH5102 Mobile Service http://contest-hunter.org:83/contest/0x50%E3%80%8C%E5%8A%A8%E6%80%81%E8%A7%84%E5%88%92%E3%80%8D%E4%BE%8B%E9%A2%98/5102%20Mobile%20Service

三人的相对顺序不重要，考虑只记录两个不确定的位置

inline int rd()
{
    rg int x(0),w(1);
    rg char c(gc);
    while(c<'0'||c>'9') {if(c=='-')w=-1;c=gc;}
    while(c>='0'&&c<='9') x=x*10+c-48,c=gc;
    return x*w;
}


void scan()
{
    l=rd();n=rd();
    FOR(i,1,l) FOR(j,1,l) c[i][j]=rd();
    FOR(i,1,n) q[i]=rd();
    memset(f,8,sizeof(f));
}


void dp()
{
    q[0]=3;f[0][1][2]=0;
    FOR(i,0,n-1) FOR(j,1,l) FOR(k,1,l) if(j!=k&&f[i][j][k]<=100000000)
    {
        if(j!=k   &&q[i+1]!=j   &&q[i+1]!=k) gmin(f[i+1][j][k],f[i][j][k]+c[q[i]][q[i+1]]);
        if(q[i]!=k&&q[i+1]!=q[i]&&q[i+1]!=k) gmin(f[i+1][q[i]][k],f[i][j][k]+c[j][q[i+1]]);
        if(q[i]!=j&&q[i+1]!=q[i]&&q[i+1]!=j)gmin(f[i+1][j][q[i]],f[i][j][k]+c[k][q[i+1]]);
    }
}

void print()
{
    FOR(i,1,l) FOR(j,1,l) ans=min(ans,f[n][i][j]);
    printf("%d",ans);
}

P1006 传纸条https://www.luogu.org/problemnew/show/P1006

两条路同时走，将能用1、2、3维确定的第4维压掉

FOR(i,1,m) FOR(j,1,n) FOR(k,max(i+j-n,1),min(i+j-1,m))
    {
        if(i==1&&j==1) f[i][j][k]=0;
        else if(i==k  && !(i==m&&j==n)) f[i][j][k]=-9999999;
        else f[i][j][k]=a[i][j]+a[k][i+j-k]+
        max(f[i-1][j][k-1],max(f[i][j-1][k-1],max(f[i-1][j][k],f[i][j-1][k])));
    }
    cout<<f[m][n][m];

CH5104 I-country http://contest-hunter.org:83/contest/0x50%E3%80%8C%E5%8A%A8%E6%80%81%E8%A7%84%E5%88%92%E3%80%8D%E4%BE%8B%E9%A2%98/5104%20I-country

凸连痛快问题，利用两个维度分别记录左右边界的单调性

const int N=16,K=230;
struct hhh{int ii,jj,ll,rr,xx,yy;}pre[N][K][N][N][2][2];
int f[N][K][N][N][2][2],a[N][N];
int n,m,k;

void scan()
{
    n=rd();m=rd();k=rd(); 
    FOR(i,1,n) FOR(j,1,m) a[i][j]=rd(); 
    //memset(f,-8,sizeof(f));
    //FOR(i,1,n) f[i][0][0][0][1][0];
    //FOR(i,1,n){FOR(j,1,m)printf("%d",a[i][j]);printf("
");}
}

void work()
{
    FOR(i,1,n) FOR(j,1,k) FOR(l,1,m) for(rg int r=l;r<=m&&r-l<j;++r)
    {
        int t=0;FOR(ii,l,r) t+=a[i][ii];
        int cz=r-l+1;
        //only widen 
        if(j==cz) 
            f[i][j][l][r][1][0]=t+f[i-1][0][0][0][1][0];
        else 
        {
            //only widen
            FOR(p,l,r) FOR(q,p,r)
                if(f[i][j][l][r][1][0]<f[i-1][j-cz][p][q][1][0])
                    f[i][j][l][r][1][0]=f[i-1][j-cz][p][q][1][0],
                    pre[i][j][l][r][1][0]={i-1,j-cz,p,q,1,0};
            //both right
            FOR(p,1,l) FOR(q,l,r) FOR(x,0,1)
                if(f[i][j][l][r][0][0]<f[i-1][j-cz][p][q][x][0])
                    f[i][j][l][r][0][0]=f[i-1][j-cz][p][q][x][0],
                    pre[i][j][l][r][0][0]={i-1,j-cz,p,q,x,0};
            //both left
            FOR(p,l,r) FOR(q,r,m) FOR(x,0,1)
                if(f[i][j][l][r][1][1]<f[i-1][j-cz][p][q][1][x])
                    f[i][j][l][r][1][1]=f[i-1][j-cz][p][q][1][x],
                    pre[i][j][l][r][1][1]={i-1,j-cz,p,q,1,x};
            //both unwiden
            FOR(p,1,l) FOR(q,r,m) FOR(x,0,1) FOR(y,0,1)
                if(f[i][j][l][r][0][1]<f[i-1][j-cz][p][q][x][y])
                    f[i][j][l][r][0][1]=f[i-1][j-cz][p][q][x][y],
                    pre[i][j][l][r][0][1]={i-1,j-cz,p,q,x,y};
            FOR(x,0,1) FOR(y,0,1) f[i][j][l][r][x][y]+=t;
        }
        //printf("%d %d %d %d %d %d %d %d %d
",i,j,l,r,t,
        //f[i][j][l][r][1][0],f[i][j][l][r][0][0],f[i][j][l][r][1][1],f[i][j][l][r][0][1]);
    }
}

void dg(hhh ans)
{
    if(!ans.jj) return;
    dg(pre[ans.ii][ans.jj][ans.ll][ans.rr][ans.xx][ans.yy]);
    FOR(i,ans.ll,ans.rr) printf("%d %d
",ans.ii,i);
}

void print()
{
    int ans=-99999999;hhh as;
    FOR(i,1,n) FOR(l,1,m) FOR(r,l,m) FOR(x,0,1) FOR(y,0,1) 
        if(ans<f[i][k][l][r][x][y])
            ans=f[i][k][l][r][x][y],
            as={i,k,l,r,x,y};
    printf("Oil : %d
",ans);
    dg(as);
}

CH5105 Cookies http://contest-hunter.org:83/contest/0x50%E3%80%8C%E5%8A%A8%E6%80%81%E8%A7%84%E5%88%92%E3%80%8D%E4%BE%8B%E9%A2%98/5105%20Cookies

代价倒序排序，由于最后每人饼干数单调递减，每次转移

a.给每人发一个

b.给i前面的k个人分别发一个

const int N=32,M=5002;
int g[N],f[N][M],a[N][M],b[N][M],c[N],ans[N],ran[N];
int n,m;

bool cmp(int x,int y){return g[x]>g[y];}

void scan()
{
    n=rd();m=rd();
    FOR(i,1,n) g[i]=rd(),ran[i]=i;
    sort(ran+1,ran+n+1,cmp);
    //FOR(i,1,n)printf("%d ",g[i]);
    FOR(i,1,n) c[i]=c[i-1]+g[ran[i]];
    memset(f,8,sizeof(f));
    f[0][0]=0;
}

void work()
{
    FOR(i,1,n) FOR(j,i,m)
    {
        f[i][j]=f[i][j-i];
        a[i][j]=i,b[i][j]=j-i;
        FOR(k,0,i-1) 
            if(f[i][j]>f[k][j-(i-k)]+k*(c[i]-c[k]))
                f[i][j]=f[k][j-(i-k)]+k*(c[i]-c[k]),
                a[i][j]=k,b[i][j]=j-(i-k);
    }
}

void dg(int n,int m)
{
    if(n==0) return;
    dg(a[n][m],b[n][m]);
    if(a[n][m]==n)
        FOR(i,1,n) ans[ran[i]]++;
    else
        FOR(i,a[n][m]+1,n) ans[ran[i]]=1;
}

void print()
{
    printf("%d
",f[n][m]);
    dg(n,m);
    FOR(i,1,n) printf("%d ",ans[i]);
}

P1541 乌龟棋 https://www.luogu.org/problemnew/show/P1541

　　FOR(i,0,c[1]) FOR(j,0,c[2]) FOR(k,0,c[3]) FOR(l,0,c[4])
    {
        if(j) f[j][k][l]=max(f[j][k][l],f[j-1][k][l]);
        if(k) f[j][k][l]=max(f[j][k][l],f[j][k-1][l]);
        if(l) f[j][k][l]=max(f[j][k][l],f[j][k][l-1]);
        f[j][k][l]+=a[i+(j<<1)+k+(k<<1)+(l<<2)];
    }
    printf("%d",f[c[2]][c[3]][c[4]]);

P1854 花店橱窗布置 https://www.luogu.org/problemnew/show/P1854

void scan()
{
    scanf("%d%d",&n,&m);
    FOR(i,1,n) FOR(j,1,m) scanf("%d",&a[i][j]);
    memset(f,-8,sizeof(f));
    FOR(i,0,n+1)f[i][0]=0;
}

void dp()
{
    FOR(i,1,m) FOR(j,1,n)
    {
        //choose
        f[i][j]=f[i-1][j],d1[i][j]=i-1,d2[i][j]=j;
        //do not choose
        if(f[i-1][j-1]+a[j][i]>f[i][j])
            f[i][j]=f[i-1][j-1]+a[j][i],
            d1[i][j]=i-1,d2[i][j]=j-1,flag[i][j]=1;
    }
}

void findpath(int i,int j)
{
    if(i==0||j==0) return;
    if(flag[i][j])FOR(ii,d1[i][j]+1,i) fl[ii]=j;
    findpath(d1[i][j],d2[i][j]);
}

void print()
{
    printf("%d
",f[m][n]);
    findpath(m,n);
    FOR(i,1,m) vase[fl[i]]=i;
    FOR(i,1,n) printf("%d ",vase[i]);
}

poj1952 BUY LOW, BUY LOWER http://poj.org/problem?id=1952

void scan()
{
    n=rd();a[0]=2e9;d[0]=1;
    FOR(i,1,n) a[i]=rd();
}

void dp()
{
    FOR(i,1,n) 
    {
        ll ans=-1;
        For(j,i-1,0) if(a[j]>a[i]) gmax(ans,f[j]);
        f[i]=ans+1;
        For(j,i-1,0) 
        {
            if(a[j]>a[i]&&f[j]==ans) d[i]+=d[j];
            else if(a[i]==a[j]&&f[i]==f[j]) break;
        }
        
    }
}

void print()
{
    ll ans=0,num;FOR(i,1,n) 
    if(ans<f[i]){ans=f[i];num=d[i];}
    else if(ans==f[i]) num+=d[i];
    printf("%lld %lld",ans,num);
}

luogu SP33 TRIP - Trip https://www.luogu.org/problemnew/show/SP33

const int N=82;
char s1[N],s2[N];
int f[N][N],d1[N][26],d2[N][26];
vector < string > s;
int n,m;

void create(int a1,int a2,int l,string ss)
{
    if(l==0) {s.push_back(ss);return;}
    FOR(i,0,25) if(f[d1[a1][i]][d2[a2][i]]==l)
        create(d1[a1][i]-1,d2[a2][i]-1,l-1,(char)('a'+i)+ss);
}

void scan()
{
    memset(s1,'',sizeof(s1));
    memset(s2,'',sizeof(s2));
    memset(f,0,sizeof(f));
    memset(d1,0,sizeof(d1));
    memset(d2,0,sizeof(d2));
    s.clear();
    scanf("%s",s1+1);scanf("%s",s2+1);
    n=strlen(s1+1),m=strlen(s2+1);
}

void dp()
{
    FOR(i,1,n) FOR(j,1,m)
    {
        if(s1[i]==s2[j]) gmax(f[i][j],f[i-1][j-1]+1);
        gmax3(f[i][j],f[i-1][j],f[i][j-1]);
    }
    FOR(i,1,n) FOR(j,0,25)
        if(s1[i]-'a'==j) d1[i][j]=i;
        else d1[i][j]=d1[i-1][j];
    FOR(i,1,m) FOR(j,0,25)
        if(s2[i]-'a'==j) d2[i][j]=i;
        else d2[i][j]=d2[i-1][j];
    create(n,m,f[n][m],"");
}

void print()
{
    sort(s.begin() , s.end());
    FOR(i,0,s.size()-1) cout<<s[i]<<endl;
    printf("
");
}