动态规划：线形 树形 矩阵

1.动态规划的基本类型

1.1线形：（分枝或分层清晰）

a.一推一型。比如等差数列，a(n)=a(n-1)+d,d为公差；等比数列，q(n)=q(n-1)*p,p为公比；累加和，s(n)=s(n-1)+n;累乘积，f(n)=f(n-1)*n;

b.多推一型。比如兔子数列f(n)=f(n-1)+f(n-2);最长不降子序列；

1.2 树形：（分枝和分层都清晰）

a.父推子（一推多）

b.子推父（多推一）

1.3 矩阵：（分层清晰）

a.一推多（我为人人）:数塔问题;

b.多推一（人人为我）：数塔问题，背包问题，最长公共子序列；

2.编程实例

2.1 等差数列

#include <bits/stdc++.h>
using namespace std;
int n, a[101], d; 
int main()
{
	cin>>n>>d;
	a[1]=100;
	for(int i=2;i<=n;i++)
	{
		a[i]=a[i-1]+d;
	}
	cout<<a[n];
	return 0;
}

2.2 等比数列

#include <bits/stdc++.h>
using namespace std;
int n, q[101], p; 
int main()
{
	cin>>n>>p;
	q[1]=1;
	for(int i=2;i<=n;i++)
	{
		q[i]=q[i-1]*p;
	}
	cout<<q[n];
	return 0;
}

2.3 累加和

#include <bits/stdc++.h>
using namespace std;
int n, s[101]; 
int main()
{
	cin>>n;
	s[1]=1;
	for(int i=2;i<=n;i++)
	{
		s[i]=s[i-1]+i;
	}
	cout<<s[n];
	return 0;
}

2.4 累乘积

#include <bits/stdc++.h>
using namespace std;
int n, f[101]; 
int main()
{
	cin>>n;
	f[1]=1;
	for(int i=2;i<=n;i++)
	{
		f[i]=f[i-1]*i;
	}
	cout<<f[n];
	return 0;
}

2.5 兔子数列

#include <bits/stdc++.h>
using namespace std;
int n, f[101]; 
int main()
{
	cin>>n;
	f[1]=1;
	f[2]=1;
	for(int i=3;i<=n;i++)
	{
		f[i]=f[i-1]+f[i-2];
	}
	cout<<f[n];
	return 0;
}

2.6 数塔问题

 此问题的解决方案有四种：逆推两种（我为人人，人人为我），顺推两种（我为人人，人人为我）

#include <bits/stdc++.h>
using namespace std;
int n;
int a[1001][1001],f[1001][1001];
int main()
{
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=i;j++)
		{
			cin>>a[i][j];
		}
	}
	//自下向上，人人为我
	for(int i=1;i<=n;i++)
	{
		f[n][i]=a[n][i];
	}	 
	for(int i=n-1;i>=1;i--)
	{
		for(int j=1;j<=i;j++)
		{
			f[i][j]=max(f[i+1][j],f[i+1][j+1])+a[i][j];
		}
	}
	cout<<f[1][1];
	return 0;
}

2.7 最长不降子序列

信息学奥赛一本通（C++版）在线评测系统 (ssoier.cn)

#include <bits/stdc++.h>
using namespace std;
#define N 1005
int n;
int a[N],f[N],p[N];
int maxL,maxi;
int main()
{
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		cin>>a[i];
		f[i]=1;
		p[i]=i;
	}
	for(int i=1;i<=n;i++)
	{
		//第一步选最大 
		int max=0;
		int k=0;
		for(int j=1;j<i;j++)
		{
			if(a[j]<=a[i] && f[j]>max)
			{
				max=f[j];
				k=j;
			}
		}
		//第二步更新 
			f[i]=max+1;
			p[i]=k;//k点更新了i点的f值 
	}
	//在f数组中选出最大值就是最终答案
	for(int i=1;i<=n;i++)
	{
		if(f[i]>maxL)
		{
			maxL=f[i];
			maxi=i;
		}
	} 
	cout<< "max="<<maxL<<endl;
	stack<int> s;
	int i=maxi;
	do
	{
		s.push(a[i]);//顺藤摸瓜 
		i=p[i];
	}
	while(p[i]!=i);
	while(!s.empty())
	{
		cout<<s.top()<<" ";
		s.pop();
	}
	return 0;
}

2.8 最长公共子序列

信息学奥赛一本通（C++版）在线评测系统

#include <bits/stdc++.h>
using namespace std;
string s,t;
int f[1010][1010] ;
void sg(string s,string t)
{
	memset(f,0,sizeof(f));
	for(int i=1;i<=s.size();i++)
	{
		for(int j=1;j<=t.size();j++) 
		{
			f[i][j]=max(f[i-1][j],f[i][j-1]);
			if(s[i-1]==t[j-1])
			{
				f[i][j]=max(f[i][j],f[i-1][j-1]+1);
			}
		}
	}
	cout<<f[s.size()][t.size()]<<endl;
}
int main()
{
	while(cin>>s>>t)
	{
		sg(s,t);
	}
	return 0;
}