算法实践学习笔记3（dp）

学习day3

基础代码

矩阵连乘

问题描述

代码实现

#include <iostream>
#include <cstring>
const int Max = 110;

int matrix[30][3] = { 0 };//输入矩阵，每组数据只会用一次
int x[Max][Max] = { 0 };//m[][],多次使用，记录最佳乘法数目
char Matrix_list[Max] = {};//矩阵表达式
void input(int n)//矩阵个数为n,输入矩阵
{
	char list;
	for (int i = 1; i <= n; i++)
	{
		std::cin >> list;
		matrix[i][0] = (int)list;
		std::cin >> matrix[i][1];
		std::cin >> matrix[i][2];
	}
}
int input_list()//输入一个序列,返回序列长度,从a[0]开始
{
	memset(Matrix_list, 0, sizeof(Matrix_list));//清零
	std::cin >> Matrix_list;
	int len = strlen(Matrix_list);
	return len;

}
int find(char b,int n)//n为a中存的矩阵个数
{
	for (int i = 1; i <= n; i++)
	{
		if (matrix[i][0] == (int)b)
			return i;
	}

}
bool check(int m,int n)//判断是否可以连乘,m为表达式长度,n为存的矩阵个数
{
	for (int i = 0; i <m-1; i++)
	{
		if (matrix[find(Matrix_list[i],n)][2] != matrix[find(Matrix_list[i + 1],n)][1])
		{
			return false;
		}
	}
	return true;
}
int MatrixChain(int n,int y)//n为表达式长度,y为存的矩阵个数
{
	
	int p[Max] = {};
	p[0] = matrix[find(Matrix_list[0],y)][1];
	for (int i = 1; i <= n; i++)
	{
		p[i] = matrix[find(Matrix_list[i-1],y)][2];
	}
	for (int i = 1; i <=n; i++)
	{
		x[i][i] = 0;
	}
	for (int r = 2; r <= n; r++)
	{
		for (int i = 1; i <=n-r+1; i++)
		{
			int j = i + r - 1;
			x[i][j] = x[i][i]+x[i + 1][j] +p[i-1]*p[i]*p[j];
			for (int k = i + 1; k < j; k++)
			{
				int t=x[i][k] +x[k+1][j]+ p[i - 1] * p[k] * p[j];
				if (t < x[i][j])
				{
					x[i][j] = t;
				}
			}
		}
	}
	return x[1][n];
}
void run(int length,int n)//length为表达式长度,n为存的矩阵个数
{
	if (check(length,n) == false)
	{
		std::cout << "MengMengDa" << std::endl;
	}
	else
	{
		std::cout << MatrixChain(length,n) << std::endl;
	}
}
void clear()
{
	memset(matrix, 0, sizeof(matrix));
	memset(x, 0, sizeof(x));//清零
}
int main()
{
	int n, m;
	while (std::cin >> n >> m)
	{
		input(n);
		for (int i = 0; i < m; i++)//输入m个矩阵序列
		{
			int len = input_list();
			run(len,n);
		}
		clear();
	}
	
	return 0;
}

总结

• 一个关键点在于如何找到表达式中矩阵的相关信息在那个位置（find函数）
  思路：匹配ASCII码值

• 这题要求输入多组数据，记得清空全局变量（clear函数）
  这里用到了memset(void *a, int c,size_t n)函数，它的作用：将a 的前 n 个字节的值设为值 c（！！该函数赋值是以字节为单位的），可以用于任何维度数组的清零，需要加上#include <cstring>

• 输入未知长度的字符串，并且返回该字符串的长度
  char a[10]就可以用cin>>a;来输入任意不超过数组设定长度的字符串，而长度可以用strlen(a)（头文件：string.h或cstring）

合并沙堆问题

问题描述

代码实现

#include <iostream>
const int Max = 310;
int a[Max] = {};//存沙子重量
int m[Max][Max] = { 0 };//记录最小代价
int p[Max][Max] = { 0 };//p[i][j]为a[i]到a[j]的和

int Minimum_cost(int n)
{
 for (int i = 1; i <= n; i++)
 {
  m[i][i] = 0;
 }
 for (int i = 1; i <=n; i++)
 {
  int x = 0;
  for (int j = i ; j <= n; j++)
  {
   p[i][j] = x + a[j];
   x = p[i][j];
   //std::cout << "!!"<< x << " ";(debug)
  }
  //std::cout << std::endl;(debug)
 }
 for (int r = 2; r <= n; r++)
 {
  for (int j = 1; j <= n - r + 1; j++)//若r=2,则依次算m12,m23,m34...
  {
   int k = j + r - 1;
   m[j][k] = m[j][j]+p[j][j] + m[j + 1][k]+p[j+1][k];
   for (int x = j + 1; x <k; x++)
   {
    //std::cout << "@"<<m[j][x] << " " << m[x + 1][k] << " " << p[j][k] << " " << p[x + 1][k] << std::endl;(debug)
    int t = m[j][x] + p[j][x] + m[x + 1][k] + p[x + 1][k];
    if (m[j][k] > t)
    {
     m[j][k] = t;
    }
   }
   //std::cout << "!" << m[j][k] << std::endl;(debug)
  }
 }
 return m[1][n];
}

void input(int n)
{
 for (int i = 1; i <= n; i++)
 {
  std::cin >> a[i];
 }
}

int main()
{
 int n;
 std::cin >> n;
 input(n);
 std::cout<<Minimum_cost(n);
 return 0;
}

总结

• 将问题的解S分成两个子问题S1和S2，S=S1+S2+合并S1和S2的代价
  合并S1和S2的代价即S1的沙堆和S2的沙堆的总质量