DP复习1

 

反正是复习，就直接写例题了。

1.Mr. Youngs Picture Permutations

要求每组数据输出一个可行的方案数。

样例数据

input

1
30
5
1 1 1 1 1
3
3 2 1
4
5 3 3 1
5
6 5 4 3 2
2
15 15
0

output

1
1
16
4158
141892608
9694845

 

可以用一个数组（a1....an）来表示每一行已经安排的学生人数，然后就只要考虑新加入一个学生所要满足的安排条件：1.ai<ni；

2.i=1或ai-1>ai./

因此我们可以用f[n][n][n][n][n]来表示从左端开始数分别占了a1,a2,a3,a4,a5人时，合影的方案数。

边界f[0][0][0][0][0]=1; 其余都为0.     目标f[n1][n2][n3][n4][n5];

DP方程：

1:若ai<ni则f[a1+1][a2][a3][a4][a5]+=f[a1][a2][a3][a4][a5];

2:若a2<n2并且a1>a2,则f[a1][a2+1][a3][a4][a5]+=f[a1][a2][a3][a4][a5];第3~5排同理；

#include<bits/stdc++.h>
using namespace std;
int k,n[6],a[6];
int main()
{
	for(;;)
	{
		scanf("%d",&k);
		if(k==0) break;
		memset(n,0,sizeof(n)); 
		for(int i=1;i<=k;i++) scanf("%d",&n[i]);
		long long f[n[1]+1][n[2]+1][n[3]+1][n[4]+1][n[5]+1];
		memset(f,0,sizeof(f));
		f[0][0][0][0][0]=1;
		for(a[1]=0;a[1]<=n[1];a[1]++)
		{
			for(a[2]=0;a[2]<=n[2];a[2]++)
			{
				for(a[3]=0;a[3]<=n[3];a[3]++)
				{
					for(a[4]=0;a[4]<=n[4];a[4]++)
					{
						for(a[5]=0;a[5]<=n[5];a[5]++)
						{
							if(a[1]<n[1]) f[a[1]+1][a[2]][a[3]][a[4]][a[5]]+=f[a[1]][a[2]][a[3]][a[4]][a[5]];
							if(a[2]<n[2]&&a[1]>a[2]) f[a[1]][a[2]+1][a[3]][a[4]][a[5]]+=f[a[1]][a[2]][a[3]][a[4]][a[5]];
							if(a[3]<n[3]&&a[2]>a[3]) f[a[1]][a[2]][a[3]+1][a[4]][a[5]]+=f[a[1]][a[2]][a[3]][a[4]][a[5]];
							if(a[4]<n[4]&&a[3]>a[4]) f[a[1]][a[2]][a[3]][a[4]+1][a[5]]+=f[a[1]][a[2]][a[3]][a[4]][a[5]];
							if(a[5]<n[5]&&a[4]>a[5]) f[a[1]][a[2]][a[3]][a[4]][a[5]+1]+=f[a[1]][a[2]][a[3]][a[4]][a[5]];
						}
					}
				}
			}
		}
		cout<<f[n[1]][n[2]][n[3]][n[4]][n[5]]<<endl;
	}
	return 0;
} 

 

 

 

2.LCIS

LCIS，就是LIS和LCS的结合。

可以用f[i][j]表示a1~ai和b1~bj可以构成以bj为结尾的的LCIS的长度。

当ai不等于bj时，f[i][j]=f[i-1][j];

当ai=bj时，f[i][j]=max (0<=k<j,bk<bj) {f[i-1][k]}+1=max (0<=k<j,bk<ai) {f[i-1][k]}+1;

就可以用三重循环来求解:

for(int i=1;i<=n;i++)
{
	for(int j=1;j<=m;j++)
	{
		if(a[i]==b[j])
		{
			for(int k=0;k<j;k++)
			{
				f[i][j]=max(f[i][j],f[i-1][k]+1);
			}
		}
		else f[i][j]=f[i-1][j];
	}
}

但是三重循环任然不行。。。。

其实在j从1增加到m的时候，i不变，所以bk<ai一只满足，当j增加一时，k的范围也增大一，从0<=k<=j变成了0<=k<=j（这个j是加过一的），所以，只要判断一下bj<ai，就可以得出当bj>=ai时，s[i][j+1]=s[i][j];当bj<ai时，s[i][j+1]=s[i][j]U{j}.

所以，只要两重循环就可以完成这道题了。

#include<bits/stdc++.h>
using namespace std;
int n,a[3100],b[3100],f[3100][3100],maxx,ans=0; 
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	for(int i=1;i<=n;i++) scanf("%d",&b[i]);
	for(int i=1;i<=n;i++)
    {
        maxx=0;
        for(int j=1;j<=n;j++) 
		{
            if(a[i]!=b[j]) f[i][j]=f[i-1][j];
            else f[i][j]=maxx+1;
            if(a[i]>b[j]) maxx=max(maxx,f[i-1][j]);
        }
    }
	for(int i=1;i<=n;i++)
	{
		ans=max(ans,f[n][i]);
	}
	printf("%d",ans);
	return 0;
}

 

 

 

3.修路Making the Grade

Sample Input

7
1
3
2
4
5
3
9

Sample Output

3

这题其实不是严格的DP，我就是用贪心水过去的， 既然要求递增或者递减，那么把两总情况的成本算出来，看看那个小就输出那个就好了。

#include<bits/stdc++.h>
using namespace std;
priority_queue<int> a;
priority_queue<int> b;
int e[2100],n,maxx=-1,ans1=0,ans2=0;
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d",&e[i]);
	for(int i=1;i<=n;i++)
	{
		a.push(e[i]);
		maxx=a.top();
		if(maxx>e[i])
		{
			ans1+=maxx-e[i];
			a.pop();
			a.push(e[i]);
		}
	}
	for(int i=n;i>=1;i--)
	{
		b.push(e[i]);
		maxx=b.top();
		if(maxx>e[i])
		{
			ans2+=maxx-e[i];
			b.pop();
			b.push(e[i]);
		}
	}
	cout<<min(ans1,ans2);
	
	return 0;
} 

 

 

 

4.Mobile Service

这题肯定会先想到用f[i][x][y][z]表示第i个请求时，三个员工的位置分别为x，y，z。

但是这样就会爆空间。

题目有这样易得条件，就是一个员工必须到请求的位置才能完成请求，所以就可以把f压成三维。

f[i][x][y]表示第i个请求的时候，三个员工分别在x，y，和第i个请求所在的位置。

方程 请求坐标用p表示，c表示花费：

f[i+1][x][y]=min(f[i+1][x][y],f[i][x][y]+c(pi，pi+1));

f[i+1][pi][y]=min(f[i+1][pi][y],f[i][x][y]+c(x,pi+1));

f[i+1][x][pi]=min(f[i+1][x][pi],f[i][x][y]+c(y,pi+1));

初始值p0=3，f[0][1][2]=0,目标max(f[n][][]);

#include<bits/stdc++.h>
using namespace std;
int l,n,c[201][201],p[1100],f[1100][201][201];
int ans=10000000;
int main()
{
    memset(f,0x3f3f3f3f,sizeof(f));
	scanf("%d%d",&l,&n);
	f[0][1][2]=0; p[0]=3; 
	for(int i=1;i<=l;i++)
	{
		for(int j=1;j<=l;j++) scanf("%d",&c[i][j]);
	}
	for(int i=1;i<=n;i++) scanf("%d",&p[i]);
	
	for(int i=0;i<n;i++)
	{
		for(int x=1;x<=l;x++)
		{
			for(int y=1;y<=l;y++)
			{
				f[i+1][x][y]=min(f[i+1][x][y],f[i][x][y]+c[p[i]][p[i+1]]);
				f[i+1][p[i]][y]=min(f[i+1][p[i]][y],f[i][x][y]+c[x][p[i+1]]);
				f[i+1][x][p[i]]=min(f[i+1][x][p[i]],f[i][x][y]+c[y][p[i+1]]);
			}
		}
	}
	for(int i=1;i<=l;i++)
	{
		for(int j=1;j<=l;j++)
		{
			ans=min(ans,f[n][i][j]);
		}
	}
	cout<<ans;
	return 0;
} 

 虽然还有几题（不，你不想写）；