动态规划：线行，树形，矩阵

1.动态规划的基本类型

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

a.一推一型：比如等比等差数列，q(n)=p*q(n-1),p为公比；a(n)=a(n-1)+d,d为公差；累加和：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[100],d;
int main()
{
	cin>>n>>d;
	a[1]=1;
	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[100],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[100];
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;
long long f[100];
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;
long long f[100];
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 数塔问题

此题有四种解法，顺推两种：人人为我，我为人人。逆推两种：人人为我，我为人人

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

#include<bits/stdc++.h>
int n,a[100][100],f[100][100]; 
using namespace std;
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;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;
}
/*
5
13
11 8
12 7 26
6 14 15 8
12 7 13 24 11
*/

2.7 最长不降子序列

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

#include<bits/stdc++.h>
using namespace std;
//人人为我动规 
#define N 2001 
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 maxf=0,maxfi=0;
		for(int j=1;j<i;j++)
		{
			if(f[j]>maxf && a[j]<=a[i])
			{
				maxf=f[j];
				maxfi=j;
			}		
		}
		//第二步更新 
		f[i]=maxf+1;
		p[i]=maxfi;
		//更新最长不降子序列的长度 
		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;
#define N 1005
string x,y;
int f[N][N];
void getst(string s,string t);
int main()
{
	while(cin>>x>>y)
	{
		getst(x,y);
	}
	return 0;
}
void getst(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;
}