首页 > 代码库 > Mark一下， dp状态转移方程写对，但是写代码都错，poj 1651 poj 1179
Mark一下， dp状态转移方程写对，但是写代码都错，poj 1651 poj 1179
dp题：
1、写状态转移方程;
2、考虑初始化边界，有意义的赋定值，还没计算的赋边界值；
3、怎么写代码自底向上计算最优值
今天做了几个基础dp，全部是dp方程写对但是初始化以及计算写错
先是poj 1651 其实就是个赤裸裸的矩阵连乘，dp方程很容易写出
dp[i][j]=min(dp[i][k]+dp[k+1][j]+r[i]*c[k]*c[j],dp[i][j]);
先贴两个个二逼的代码，mark下自己多么的二逼：
二逼一：在计算的时候使用了还没有算出来的值，模拟下就知道第一重循环里算dp[0][2]=dp[0][0]+dp[1][2];
!!!!!dp[1][2]这时候根本没有计算出来还，Tmd我就用了这个值！！！！！！！！！！！！！！！！！！！！！！！二逼啊！！！！！！！！！
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int SIZE = 100+10;
const int INF = 10000000;
int dp[SIZE][SIZE];
int r[SIZE],c[SIZE];
int main()
{
int n;
while(~scanf("%d",&n))
{
for(int i=0;i<n-1;i++)
{
scanf("%d",&r[i]);
if(i>0)c[i-1]=r[i];
}
scanf("%d",&c[n-2]);
n--;
int ans=0;
for(int i=0;i<n;i++)dp[i][i]=0;
for(int i=0;i<=n;i++)
for(int j=i;j<n;j++)
{
for(int k=i;k<j;k++)
dp[i][j]=min(dp[i][k]+dp[k+1][j]+r[i]*c[k]*c[j],dp[i][j]);
}
for(int i=0;i<n;i++)
{
for(int j=0;j<n;j++)
printf("%d ",dp[i][j]);
putchar('\n');
}
printf("%d\n",dp[0][n-1]);
}
return 0;
}
二逼二：来回换思路的时候各种思路弄混，看注释
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int SIZE = 100+10;
const int INF = 10000000;
int dp[SIZE][SIZE];
int r[SIZE],c[SIZE];
int main()
{
int n;
while(~scanf("%d",&n))
{
for(int i=0;i<n-1;i++)
{
scanf("%d",&r[i]);
if(i>0)c[i-1]=r[i];
}
scanf("%d",&c[n-2]);
n--;
int ans=0;
for(int i=0;i<n;i++)dp[i][i]=0;
int j;
for(int l=2;l<=n;l++)
for(int i=0;i<l;i++)//此处i<l及其傻逼，请勿模仿
{
j=i+l-1;
dp[i][j]=INF;
for(int k=i;k<j;k++)
dp[i][j]=min(dp[i][k]+dp[k+1][j]+r[i]*c[k]*c[j],dp[i][j]);
}
/*for(int i=0;i<n;i++)
{
for(int j=0;j<n;j++)
printf("%d ",dp[i][j]);
putchar('\n');
}*/
printf("%d\n",dp[0][n-1]);
}
return 0;
}
正确代码：
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define ll long long
const int SIZE = 100+10;
const int INF = 200000000;
ll dp[SIZE][SIZE];
int r[SIZE],c[SIZE];
int main()
{
int n;
while(~scanf("%d",&n))
{
for(int i=0;i<n-1;i++)
{
scanf("%d",&r[i]);
if(i>0)c[i-1]=r[i];
}
scanf("%d",&c[n-2]);
n--;
int ans=0;
for(int i=0;i<=n;i++)dp[i][i]=0;
int j;
for(int l=2;l<=n;l++)
for(int i=0;i<n-l+1;i++)
{
j=i+l-1;
dp[i][j]=INF;
for(int k=i;k<j;k++)
dp[i][j]=min(dp[i][k]+dp[k+1][j]+r[i]*c[k]*c[j],dp[i][j]);
}
printf("%lld\n",dp[0][n-1]);
}
return 0;
}
在做这个poj 1651之前 先做的poj 1179，然后过了poj 1651之后改了poj1179原来的循环次序，其实就是类比这个poj1651,第一重循环弄成长度即可
注意 负数乘以负数是正数 考虑最大最小值得时候一定注意
不过重新检查下发现代码在初始化的时候冗余，精简了下。
贴代码：
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define ll long long
const int SIZE =55;
const int INF = (-2147483647 ) ;
ll dp[SIZE][SIZE],dpmin[SIZE][SIZE];
int pos[SIZE];
bool op[SIZE];
int main()
{
//freopen("poj1179.txt","r",stdin);
int n;
ll tmp1,tmp2;
while(scanf("%d",&n)!=EOF)
{
for(int i=0;i<n;i++)dp[i][0]=dpmin[i][0]=0;
for(int i=0;i<n;i++)
{
scanf(" %c %lld",&op[i],&dp[i][1]);
dpmin[i][1]=dp[i][1];
}
for(int j=1;j<=n;j++)
for(int i=0;i<n;i++)
{
if(j>1)dp[i][j]=INF,dpmin[i][j]=-INF;
for(int k=1;k<j;k++)
{
if(op[(i+k)%n ] == 't')
{
dp[i%n][j]=max(dp[i%n][k]+dp[(i+k)%n ][j-k],dp[i%n][j]);
dpmin[i%n][j]=min(dpmin[i%n][k]+dpmin[(i+k)%n ][j-k], dpmin[i%n][j]);
}
else
{
tmp1=max(dp[i%n][k]*dp[(i+k)%n ][j-k],dpmin[i%n][k]*dpmin[(i+k)%n ][j-k]);
tmp2=max(dp[i%n][k]*dpmin[(i+k)%n ][j-k], dpmin[i%n][k]*dp[(i+k)%n ][j-k]);
tmp1=max(tmp1,tmp2);
dp[i%n][j]=max(tmp1,dp[i%n][j]);
tmp2=min(dp[i%n][k]*dp[(i+k)%n ][j-k],dpmin[i%n][k]*dpmin[(i+k)%n ][j-k]);
tmp1=min(dp[i%n][k]*dpmin[(i+k)%n ][j-k], dpmin[i%n][k]*dp[(i+k)%n ][j-k]);
tmp2=min(tmp1,tmp2);
dpmin[i%n][j]=min(dpmin[i%n][j],tmp2);
}
}
}
int cnt=0;
ll ans=INF;
for(int i=0;i<n;i++)//开始的时候，这里i<=n，于是初始化的时候必须优dp[n][...]=INF,才能AC，奇怪的是计算过程中第一个下标我都%n了，最后才看到这里的问题
{
ans=max(ans,dp[i][n]);
}
for(int i=0;i<n;i++)
if(ans == dp[i][n])
pos[cnt++]=i+1;
printf("%lld\n",ans);
for(int i=0;i<cnt;i++)
if(i == cnt-1)printf("%d\n",pos[i]);
else printf("%d ",pos[i]);
}
return 0;
}
声明：以上内容来自用户投稿及互联网公开渠道收集整理发布，本网站不拥有所有权，未作人工编辑处理，也不承担相关法律责任，若内容有误或涉及侵权可进行投诉： 投诉/举报 工作人员会在5个工作日内联系你，一经查实，本站将立刻删除涉嫌侵权内容。