1、矩阵连乘,tyvj 1198(最优矩阵连乘),关键是写出动态转移方程。用DP[i][j]表示矩阵Ai乘到Aj的最优解,P[]用来存储矩阵的行和列,M[i-1]表示矩阵i的行,M[i]表示矩阵i的列。
当i==j时,DP[i][i]=0;
当i<j时,可以假设他从k位置隔开。比方说从A3..A8。那么A1A2(A3A4A5A6A7A8)A9。我们知道若A矩阵为r1*c的矩阵,A2为c*r2的矩阵,得到的将是一个r1*r2的矩阵,且连乘次数为r1*c*r2。假设我从5位置隔开,就可以分为A1A2((A3A4A5)(A6A7A8))A9。那么这个时候的次数为:DP[3][5]+DP[6][8]+M[2]*M[5]*M[8](A3的行*A5的列*A8的列)。那么我只需要模举这其中的组合然后取最小值即可,这个时候的动态转移方程为:DP[i][j]=min(DP[i][k]+DP[k+1][j]+M[i-1]*M[k]*M[j])(i<=k<j)。
问题是:你要知道DP[i][k]和DP[k+1][j]的值,你才可以和DP[i][j]做比较,才能求出DP[i][j]的值。
思想是先求出相邻组的连乘次数,如一开始算相邻两个相乘,如图:
那么每两个相乘的结果我们可以得到。同理继续求三个相邻的结果。
计算的话比方说要求A1A2A3,可以化为(A1)(A2A3)=DP[1][1]+DP[2][3]+M[0]*M[1]*M[3],依次类推结果为:DP[i][j]=DP[i][i]+DP[i+1][j]+M[i-1]*M[i]*M[j]。
只有这样做后我才可以求任意组合的次数,求出三个一组的结果后,我就可以求任意的四个矩阵连乘的最小值。如:A4A5A6A7,随便怎么组合,通过上面是计算都可以求出相应的结果。
现在讨论下范围:用r保存要求的每每相邻的组数(2<=r<=n)(从2开始是这里前提条件是i<j),i的范围为:(1<=i<=n-r+1,稍作推理就知道,当r=3时,即三个三个一组时i<=9-3+1=7),j的值固定j=r+i-1,具体见代码:
- #include<iostream>
- #include<cstring>
- #include<cstdio>
- using namespace std;
- const int MAX=110;
- #define min(a,b) (a)<(b)?(a):(b)
- #define CLR(arr,val) memset(arr,val,sizeof(arr))
- int n,M[MAX],DP[MAX][MAX];
- void Matrix()
- { for(int r=2;r<=n;r++)
- for(int i=1;i<=n-r+1;i++)
- { int j=r+i-1;
- DP[i][j]=DP[i+1][j]+M[i-1]*M[i]*M[j];
- for(int k=i;k<j;k++)
- DP[i][j]=min(DP[i][k]+DP[k+1][j]+M[i-1]*M[k]*M[j],DP[i][j]);
- }
- }
- int main()
- { scanf("%d",&n);
- for(int i=0;i<=n;i++)
- scanf("%d",&M[i]);
- CLR(DP,0);
- Matrix();
- cout<<DP[1][n]<<endl;
- return 0;
- }
路径输出:
- #include<iostream>
- #include<cstring>
- #include<cstdio>
- using namespace std;
- const int MAX=110;
- #define CLR(arr,val) memset(arr,val,sizeof(arr))
- int n,M[MAX],DP[MAX][MAX],s[MAX][MAX];
- void Matrix()
- { for(int r=2;r<=n;r++)
- for(int i=1;i<=n-r+1;i++)
- { int j=r+i-1;
- DP[i][j]=DP[i+1][j]+M[i-1]*M[i]*M[j];
- s[i][j]=i; //i位置断开
- for(int k=i;k<j;k++)
- { int temp=DP[i][k]+DP[k+1][j]+M[i-1]*M[k]*M[j];
- if(temp<DP[i][j])
- { s[i][j]=k;//k位置断开
- DP[i][j]=temp;
- }
- }
- }
- }
- void Find(int i,int j)
- { if(i==j)
- { cout<<'A'<<i;
- return ;
- }
- if(i<s[i][j]) cout<<"(";
- Find(i,s[i][j]);
- if(i<s[i][j]) cout<<")";
- if(s[i][j]+1<j) cout<<"(";
- Find(s[i][j]+1,j);
- if(s[i][j]+1<j) cout<<")";
- }
- int main()
- { scanf("%d",&n);
- for(int i=0;i<=n;i++)
- scanf("%d",&M[i]);
- CLR(DP,0);
- Matrix();
- Find(1,n);
- cout<<endl<<DP[1][n]<<endl;
- return 0;
- }
2、NYOJ 15(括号匹配),令DP[i][j]表示字符串s的第i个字符到第j个字符需要添加的最少括号数。
①、当i==j时,只有一个符号,需再添加一个字符,所以DP[i][i]=1;
②、当i<j时,若是s[i]与S[j]配对,那么DP[i][j]=DP[i+1][j-1];
若是不配对,那么从中任选一个字符作为分界点,i<=k<j,就有:DP[i][j]=min(DP[i][j],DP[i][k]+DP[k][j])。具体实现和上例一样:
- #include<iostream>
- #include<string>
- #include<cstring>
- using namespace std;
- const int MAX=1010;
- #define min(a,b) (a)<(b)?(a):(b)
- #define CLR(arr,val) memset(arr,val,sizeof(arr))
- string s;
- int n,DP[MAX][MAX];
- bool pipei(char c1,char c2)
- { if(c1=='('&&c2==')'||c1=='['&&c2==']') return true;
- return false;
- }
- int main()
- { cin>>n;
- cin.get();
- while(n--)
- { cin>>s;
- CLR(DP,0);
- int len=s.length();
- for(int i=1;i<=len;i++) DP[i][i]=1;
- for(int r=2;r<=len;r++)
- for(int i=1;i<=len-r+1;i++)
- { int j=r+i-1;
- DP[i][j]=MAX; //一定要赋值
- if(pipei(s[i-1],s[j-1])) DP[i][j]=DP[i+1][j-1];
- for(int k=i;k<j;k++)
- DP[i][j]=min(DP[i][k]+DP[k+1][j],DP[i][j]);
- }
- cout<<DP[1][len]<<endl;
- }
- return 0;
- }
3、石子合并:Tyvj 1055(石子合并),和这个题目类似的是Tyvj 1062(合并傻子)。
- #include<iostream>
- #include<cstring>
- using namespace std;
- const int MAX=1010;
- #define CLR(arr,val) memset(arr,val,sizeof(arr))
- int n,a[MAX],sum[MAX][MAX],DP[MAX][MAX];
- int Sand()
- { for(int r=2;r<=n;r++)
- for(int i=1;i<=n-r+1;i++)
- { int j=i+r-1;
- DP[i][j]=DP[i+1][j]+sum[i][j];
- for(int k=i;k<j;k++)
- DP[i][j]=min(DP[i][j],DP[i][k]+DP[k+1][j]+sum[i][j]);
- }
- return DP[1][n];
- }
- int main()
- { cin>>n;
- for(int i=1;i<=n;i++)
- cin>>a[i];
- CLR(DP,0);
- CLR(sum,0);
- for(int i=1;i<=n;i++)
- for(int j=i;j<=n;j++)
- sum[i][j]=sum[i][j-1]+a[j];
- cout<<Sand()<<endl;
- return 0;
- }
4、加分二叉树。Tyvj 1073(加分二叉树)。
- #include<iostream>
- #include<cstring>
- using namespace std;
- const int MAX=50;
- #define CLR(arr,val) memset(arr,val,sizeof(arr))
- long n,a[MAX],DP[MAX][MAX],s[MAX][MAX];
- int flag=1;
- long Max()
- { for(int r=2;r<=n;r++)
- for(int i=1;i<=n-r+1;i++)
- { int j=i+r-1;
- DP[i][j]=DP[i+1][j]+a[i];
- s[i][j]=i;
- for(int k=i+1;k<j;k++)
- { int temp=DP[i][k-1]*DP[k+1][j]+a[k];
- if(temp>DP[i][j])
- { DP[i][j]=temp;
- s[i][j]=k;
- }
- }
- }
- return DP[1][n];
- }
- void Search(int L,int R)
- { if(L>R) return ;
- if(flag) flag=0;
- else cout<<" ";
- cout<<s[L][R];
- if(L==R) return ;
- Search(L,s[L][R]-1);
- Search(s[L][R]+1,R);
- }
- int main()
- { cin>>n;
- CLR(DP,0);
- for(int i=1;i<=n;i++)
- { cin>>a[i];
- DP[i][i]=a[i];
- }
- cout<<Max()<<endl;
- for(int i=1;i<=n;i++)
- s[i][i]=i; //子树为空,本身就为树的根节点
- Search(1,n);
- cout<<endl;
- return 0;
- }
题5:牢房:
描述
SB王国中有一个奇怪的监狱,这个监狱里一共有P间牢房,这些牢房一字排开,从左到右按1到P进行编号,第i间后面紧挨着第(i+1)间(最后一间除外)。现在有P名罪犯被关押在这P间牢房里。某一天,上级下发了一个释放名单,要求每天释放名单上的一个人。这可把监狱中的看守们吓得不轻,因为看守们知道,现在牢房中的P个人,可以相互之间传话。如果某个人离开了,那么原来和这个人能够传上话的人都会很气愤,导致他们那天会一直大吼大叫,搞得看守们很是头疼,但是如果给这些要发火的人吃上肉,他们就会安静。现在看守们想知道,如何安排释放的顺序,才能使得他们花费的肉钱最少。1<=p<=1000,1<=Q<=100
输入
第一行两个数P和Q,Q表示要释放名单上的人数;
第二行Q个数,表示释放哪些人
输出
仅一行,表示最少给多少人次送肉.
样例输入
20 3
3 6 14
样例输出
35
提示:(样例说明)
先释放14号监狱中的罪犯,要给1到13号监狱和15到20号监狱中的19人送肉吃;再释放6号监狱中的罪犯,要给1到5号监狱和7到13号监狱中的12人送肉吃;最后释放3号监狱中的罪犯,要给1到2号监狱和4到5号监狱中的4人送肉吃。
分析:用DP[i][j]表示在[i,j]区间所有该被释放的人都被释放的最小支付代价。
若是释放k,那么代价为DP[i][k-1]+DP[k+1][j]+(j-i),所以动态转移方程为:DP[i][j]=min(DP[i][k-1]+DP[k+1][j]+(j-i)).
- #include<iostream>
- #include<cstring>
- #include<cstdio>
- #include<algorithm>
- using namespace std;
- const int MAX=1010;
- const int Inf=100000000;
- int n,m,a[MAX],sum[MAX],DP[MAX][MAX];
- int main()
- { for(int i=0;i<MAX;i++)
- { fill(DP[i],DP[i]+MAX,Inf);
- DP[i][i]=0;
- }
- scanf("%d%d",&n,&m);
- for(int i=1;i<=m;i++)
- scanf("%d",&a[i]);
- sort(a+1,a+m+1);
- for(int i=1;i<=m;i++)
- sum[i]=a[i]-a[i-1]-1;
- sum[m+1]=n-a[m];
- for(int i=1;i<=m+1;i++)
- sum[i]=sum[i]+sum[i-1];
- for(int i=m+1;i>=1;i--)
- for(int j=i+1;j<=m+1;j++)
- for(int k=i;k<j;k++)
- DP[i][j]=min(DP[i][j],DP[i][k]+DP[k+1][j]+sum[j]-sum[i-1]+j-i-1);
- printf("%d\n",DP[1][m+1]);
- return 0;
- }