【状压DP】周长

题目描述

给出n个不同长度，宽度都为1的矩形，现在可以叠起来，然后要求最大的周长以及周长最长的方案数

样例输入

4
1 2 3 4

样例输出

20 8

思路

这道题DP还是很容易想到的吧
就是很麻烦的是怎么记录当前状态是什么矩形用过什么没用过
那么我们就用一个01串来处理
设$f_{S,i}$表示当前状态为S，当前矩形是i的最大周长（不记录上下面）
那么转移就可以从所有S中为0的矩形去转移

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define ll long long

using namespace std;

int n, m;
ll ans, ans_;
ll a[55], f[50005][55], sum[50005][55];

int main()
{
	scanf("%d", &n);
	for(int i = 1; i <= n; ++i) {
		scanf("%lld", &a[i]);
		ll j = (1 << (i - 1));
		f[j][i] = 2 * a[i];
		sum[j][i] = 1;//初始化
	}
	for(int S = 0; S < (1 << n); ++S)
	for(int i = 1; i <= n; ++i)
	if((S >> (i - 1)) & 1)
	{
		for(int j = 1; j <= n; ++j)//枚举下一个矩形
		{
			if(!(S >> (j - 1) & 1)) {//选取的下一个矩形一定是要没有用过的
				int t = S | (1 << (j - 1));
				ll val = f[S][i] + 2 * max((ll)0, a[j] - a[i]);
				if(f[t][j] < val)
				{
					f[t][j] = val;
					sum[t][j] = sum[S][i];
				}
				else if(f[t][j] == val)
					sum[t][j] += sum[S][i];//累计方案数
			}
		}
	}
	ll j = (1 << n) - 1;
	for(int i = 1; i <= n; ++i)
	{
		if(f[j][i] > ans)
		{
			ans = f[j][i];
			ans_ = sum[j][i];	
		}	
		else if(f[j][i] == ans)
			ans_ += sum[j][i];
	}
	printf("%lld %lld", ans + 2 * n, ans_);
	return 0;
}