描述
佳佳最近又迷上了某款类似于虚拟人生的网络游戏。在游戏中，佳佳是某旅行团的团长，他需要安排客户住进旅馆。旅馆给了佳佳的旅行团一个房间数的限制。每一个房间有不同的容纳人数和价钱（这个价格是房间的总价格，不是每个人付的）。佳佳决定找到最小的花费，安排参加旅行的人住在这里。但是他遇到了这么一个问题：两个不同性别的人不能住在同一个房间里，除非他们是夫妻；一对夫妻如果在一起住，那么别的人就不能再住进去。你不必让所有的夫妻都单独住在一起。也就是说：
1.给你一些房间，告诉你这些房间的容纳人数和价格
2.安排一定数量的人住到旅馆里，满足：
a.不同性别的人如果不是夫妻那么不能住一起。
b.夫妻如果住在一起，那么房间不能安排其他的人进去。
你来写一个程序帮助佳佳找到安排这些来参加旅行的人住进旅馆所需要的最小花费。
格式
输入格式
第一行有4个用空格隔开的整数m,f,r,c，分别表示参加旅行的男性人数、参加旅行的女性人数、旅馆的房间数、这些男女中有多少对夫妻。注意每一个人不是单身就是和他/她唯一的妻子/丈夫一起参加旅行。
接下来有r行，每行描述了一个房间。每行有两个整数Bi，Pi，它们分别表示每一个房子的容纳人数和价格（无论住多少人，房间的价格不变）。
对于30%的数据，0<=m,f,r<=50；
对于100%的数据，0<=m,f,r<=300，0<=c<=Min(m,f)，0<=Bi,Pi<=10。
输出格式
输出为旅行的人订购房间所需要的最小花费。如果没有这样的安排，请输出“Impossible”代替。
限制
各测试点5秒
-----------------------------------------
总存在一个最优解使得只有一对couple在一起或者没有，因为假如有两队，可以男女交换
f[i][j][k]和g[i][j][k]分别表示没有和有一对，转移时除了f向f，g向g外，再考虑g[i][j][k]=min(g[i][j][k],f[i-1][j-1][k-1]+p[i])
【WARN】：1.注意这个转移g[i][j][k]=min(g[i][j][k],f[i-1][j-1][k-1]+p[i]) j和k>=1
2.房间不一定注满人，如果再加一层枚举房间住几个人会很慢，可以把j和k<0的情况从不合法改成从0转移，想当于考虑了房间不满，最多男女各有一个房间不满
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=,INF=1e9;
int m,fm,r,c,b[N],p[N];
int f[N][N][N],g[N][N][N]; void dp(){
for(int j=;j<=m;j++) for(int k=;k<=fm;k++) f[][j][k]=g[][j][k]=INF;
f[][][]=;
for(int i=;i<=r;i++)
for(int j=;j<=m;j++)
for(int k=;k<=fm;k++){
f[i][j][k]=f[i-][j][k]; g[i][j][k]=g[i-][j][k];
if(b[i]>=&&j>=&&k>=) g[i][j][k]=min(g[i][j][k],f[i-][j-][k-]+p[i]);
int sj=max(,j-b[i]),sk=max(,k-b[i]);
f[i][j][k]=min(f[i][j][k],f[i-][sj][k]+p[i]),
g[i][j][k]=min(g[i][j][k],g[i-][sj][k]+p[i]);
f[i][j][k]=min(f[i][j][k],f[i-][j][sk]+p[i]),
g[i][j][k]=min(g[i][j][k],g[i-][j][sk]+p[i]);
}
}
int main(){
scanf("%d%d%d%d",&m,&fm,&r,&c);
for(int i=;i<=r;i++) scanf("%d%d",&b[i],&p[i]);
dp(); int ans;
if(c==) ans=f[r][m][fm];
else ans=min(f[r][m][fm],g[r][m][fm]);
if(ans!=INF) cout<<ans;
else cout<<"Impossible";
}
//题解中神犇的滚动数组做法 边读边做
//快了近三倍
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
typedef long long LL;
const int N = + ;
int m,f,r,c;
int dp[N][N][];
int main()
{
scanf("%d%d%d%d",&m,&f,&r,&c); memset(dp,0x3f3f3f3f,sizeof(dp));
dp[m][f][] = dp[m][f][] = ;
for(int i = ; i < r; i++)
{
int a,b;
scanf("%d%d",&a,&b);
for(int i = ; i <= m; i++)
{
for(int j = ; j <= f; j++)
{
int x1 = max(,i-a);
int x2 = max(,j-a);
dp[x1][j][] = min(dp[x1][j][],dp[i][j][] + b);
dp[i][x2][] = min(dp[i][x2][],dp[i][j][] + b);
dp[x1][j][] = min(dp[x1][j][],dp[i][j][] + b);
dp[i][x2][] = min(dp[i][x2][],dp[i][j][] + b);
if(i > && j > && a >= )
dp[i-][j-][] = min(dp[i-][j-][],dp[i][j][] + b);
}
}
}
int d = min(dp[][][],dp[][][]);
if(d == 0x3f3f3f3f)
printf("Impossible");
else
printf("%d",d); }