博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
dp四边形优化
阅读量:5291 次
发布时间:2019-06-14

本文共 1519 字,大约阅读时间需要 5 分钟。

如果dp[i][j]=min(dp[i][k]+dp[k+1][j]+w[i][j]);且满足dp[a][c]+dp[b][d]<=dp[a][d]+dp[c][d](a<b<=c<d);

那么dp具有四边形不等式性质

另外如果可以证明w[i][j]满足单调性和四边形不等式性质,那么dp也具有四边形不等式性质

单调性:w[i][j]<=w[i][j+1]<=w[i+1][j+1]

 

http://www.51nod.com/Challenge/Problem.html#problemId=1022

 

#include
#define fi first#define se second#define INF 0x3f3f3f3f#define LNF 0x3f3f3f3f3f3f3f3f#define fio ios::sync_with_stdio(false);cin.tie(0);cout.tie(0)#define pqueue priority_queue#define NEW(a,b) memset(a,b,sizeof(a))const double pi=4.0*atan(1.0);const double e=exp(1.0);const int maxn=1e6+8;typedef long long LL;typedef unsigned long long ULL;//typedef pair
P;const LL mod=1e9+7;const ULL base=1e7+7;using namespace std;LL dp[2008][2008];LL s[2008][2008];LL a[2008];int main(){ int n; int x; scanf("%d",&n); for(int i=1;i<=n;i++){ scanf("%lld",&a[i]); a[i+n]=a[i]; } n*=2; for(int i=1;i<=n;i++){ a[i]+=a[i-1]; } memset(dp,LNF,sizeof(dp)); for(int i=1;i<=n;i++){ s[i][i]=i; dp[i][i]=0; } for(int i=n;i>=1;i--){ for(int j=i+1;j<=n;j++){ for(int k=s[i][j-1];k<=s[i+1][j];k++){ if(dp[i][j]>dp[i][k]+dp[k+1][j]+a[j]-a[i-1]){ dp[i][j]=dp[i][k]+dp[k+1][j]+a[j]-a[i-1]; s[i][j]=k; } } } } n/=2; LL ans=LNF; for(int i=1;i<=n;i++){ ans=min(ans,dp[i][i+n-1]); } printf("%lld\n",ans);}

 

转载于:https://www.cnblogs.com/Profish/p/11200467.html

你可能感兴趣的文章
程序员的“机械同感”
查看>>
在16aspx.com上下了一个简单商品房销售系统源码,怎么修改它的默认登录名和密码...
查看>>
c++回调函数
查看>>
linux下Rtree的安装
查看>>
【Java】 剑指offer(53-2) 0到n-1中缺失的数字
查看>>
Delphi中ListView类的用法
查看>>
bzoj3110: [Zjoi2013]K大数查询 【树套树,标记永久化】
查看>>
[原创]Java 的传值小例子
查看>>
【MySQL学习】安装和配置 服务无法启动 没有报告任何错误
查看>>
C# 修饰符
查看>>
JavaScript启示录
查看>>
我需要什么样的浏览器?
查看>>
取textaera里的值
查看>>
java设计模式1--工厂方法模式(Factory Method)
查看>>
博客第一弹—聊聊HTML的那些事
查看>>
上海2017QCon个人分享总结
查看>>
HIVE快速入门 分类: B4_HIVE 2015-...
查看>>
Mysql安装方法及安装问题解决
查看>>
Java动态代理的两种实现方式:
查看>>
PHP trait
查看>>