BZOJ.5311.贞鱼(DP 决策单调)

题目链接

很容易写出\(O(n^2k)\)的DP方程。然后显然决策点是单调的，于是维护决策点就可以了。。
这个过程看代码或者别的博客吧我不写了。。（其实是忘了）
这样复杂度\(O(nk\log n)\)。但是在BZOJ T了=-=。
\(k\)可以带权二分优化到\(O(n\log k\log n)\)就能过了吧。
不想改了。

我特么学的是假的单调。。
又是zz错误浪费半下午（╯‵□′）╯︵┴─┴
辣鸡题还卡时间
不过就不过吧mmp

Upd: Codeforces 321E.Ciel and Gondolas上过了。1028ms 130224KB。。

//130968kb  3000ms
#include <cstdio>
#include <cctype>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 5000000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
const int N=4003;

int n,K,A[N][N],sum[N][N],f[2][N],Now;
char IN[MAXIN],*SS=IN,*TT=IN;
struct Node{
    int l,r,pos;//pos是区间[l,r]的最优转移点 
}q[N];

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now;
}
inline int Cost(int i,int p){//在i之前，分割p,p+1处 
    return f[Now][p]+sum[p+1][i];
}
inline int Find(Node t,int x)
{
    int l=t.l, r=t.r, mid;
    while(l<=r)//l==r时应再Check一次？
        if(mid=l+r>>1, Cost(mid,x)<Cost(mid,t.pos)) r=mid-1;//!
        else l=mid+1;
    return l;
}

int main()
{
    n=read(),K=read();
    for(int i=1; i<=n; ++i)
        for(int j=1; j<=n; ++j) A[i][j]=A[i][j-1]+read();
    for(int i=1; i<=n; ++i)
        for(int j=i+1; j<=n; ++j)
            sum[i][j]=sum[i][j-1]+A[j][j]-A[j][i-1];
    for(int i=1; i<=n; ++i) f[1][i]=sum[1][i];
    Now=1;
    for(int j=1; j<K; ++j, Now^=1)
    {
        int h=1,t=1; q[1]=(Node){1,n,1};
        for(int i=2; i<=n; ++i)
        {
            if(i>q[h].r) ++h;
            f[Now^1][i]=Cost(i,q[h].pos);
            if(Cost(n,i)<Cost(n,q[t].pos))//为什么要拿n比？不太明白。
            {
                while(h<=t && Cost(q[t].l,i)<Cost(q[t].l,q[t].pos)) --t;//队尾区间的l用i都比pos更优了，而决策点是单调的，所以[l,r]肯定都要不选pos而选i了 
                if(h>t) q[++t]=(Node){i,n,i};
                else
                {
                    int Pos=Find(q[t],i);
                    q[t].r=Pos-1, q[++t]=(Node){Pos,n,i};
                }
            }
        }
    }
    printf("%d",f[Now][n]);

    return 0;
}

转载于:https://www.cnblogs.com/SovietPower/p/9163540.html