首页 > 代码库 > 【BZOJ 1592】[Usaco2008 Feb]Making the Grade 路面修整 dp优化之转移变状态
【BZOJ 1592】[Usaco2008 Feb]Making the Grade 路面修整 dp优化之转移变状态
我们感性可证离散(不离散没法做),于是我们就有了状态转移的思路(我们只考虑单不减另一个同理),f[i][j]到了第i块高度为j的最小话费,于是我们就可以发现f[i][j]=Min(f[i-1][k])+|a[i]-j|(k<=j),于是我们的思路就去了各种数据结构…….然后我们发现对于这些转移就是在记录小于等于,那么我们直接带状态里体现这一点就可以了,而不是在转移的时候,我们f[i][j]表示到了第i个点小于等于j的高度的最小花费,这样我们就n^2了。
#include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #define N 2010 using namespace std; inline int read() { register int sum=0;register char ch=getchar(); while(ch<‘0‘||ch>‘9‘)ch=getchar(); while(ch>=‘0‘&&ch<=‘9‘)sum=(sum<<1)+(sum<<3)+ch-‘0‘,ch=getchar(); return sum; } int a[N],f[N][N],ans,n,pos[N],len,Hash[N]; int comp(const int x,const int y){ return a[x]<a[y]; } inline int Min(int x,int y){ return x<y?x:y; } inline int Abs(int x){ return x<0?-x:x; } inline int get_Min(){ memset(f,0x7f,sizeof(f));register int ans=0x7fffffff; for(register int i=1;i<=len;i++)f[1][i]=Min(f[1][i-1],Abs(Hash[i]-Hash[a[1]])); for(register int i=2;i<=n;i++) for(register int j=1;j<=len;j++) f[i][j]=Min(f[i][j-1],f[i-1][j]+Abs(Hash[j]-Hash[a[i]])); for(register int i=1;i<=len;i++)ans=Min(ans,f[n][i]); return ans; } inline int get_Max(){ memset(f,0x7f,sizeof(f));register int ans=0x7fffffff; for(register int i=1;i<=len;i++)f[1][i]=Min(f[1][i+1],Abs(Hash[i]-Hash[a[1]])); for(register int i=2;i<=n;i++) for(register int j=len;j>0;j--) f[i][j]=Min(f[i][j+1],f[i-1][j]+Abs(Hash[j]-Hash[a[i]])); for(register int i=1;i<=len;i++)ans=Min(ans,f[n][i]); return ans; } int main(){ n=read();for(register int i=1;i<=n;i++)a[i]=read(),pos[i]=i; sort(pos+1,pos+n+1,comp); for(register int i=1;i<=n;i++) if(i==1||a[pos[i]]!=a[pos[i-1]])Hash[++len]=a[pos[i]],a[pos[i]]=len; else a[pos[i]]=len; printf("%d",Min(get_Min(),get_Max())); }
【BZOJ 1592】[Usaco2008 Feb]Making the Grade 路面修整 dp优化之转移变状态
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。