[CF744C]Hongcow Buys a Deck of Cards

题目

传送门 to luogu

思路

看到 $n\le 16$ 应该没有别的想法。但问题是，$r_i,b_i$ 很大，如果用背包的思路，容量开不下……

所以什么东西的数据范围很小呢？是 隐藏条件 $\mu =1$ 。这个 $\mu$ 是 最大静摩擦系数 一张卡将会带来的花费减免。也就是说，不管 $r_i,b_i$ 多大，我的实际花费跟 $\sum r$ 和 $\sum b$ 不会偏差太多。

所以用 $f(S,r^{\prime })=b^{\prime }$ 表示 红宝石实际花费比 $\sum r$ 小 $r^{\prime }$ 的情况下，蓝宝石节省（当然，作为取值本身的定义就无所谓了，你也可以定义成最少蓝宝石花费）$b^{\prime }$ 个。

状态数 $\mathcal{O}(2^n{n}^2)$ ，转移是 $\mathcal{O}(n)$ 的。还是挺慢的……但是时限 $2\mathrm{s}$ 帮你略去了卡常环节。

如果你的定义是蓝宝石节省数量，你可以开$\sout{\mathrm{s}\mathrm{h}\mathrm{o}\mathrm{r}\mathrm{t}}$节省空间！

代码

#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;
inline int readint(){
	int a = 0; char c = getchar(), f = 1;
	for(; c<'0'||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}

const int MaxN = 16;
int dp[1<<MaxN][MaxN*MaxN];
int R[1<<MaxN], B[1<<MaxN]; // cnt of card
int r[MaxN], b[MaxN]; // cost

char inBuf[MaxN];
int main(){
	int n = readint();
	int allR = 0, allB = 0;
	for(int i=0; i<n; ++i){
		scanf("%s",inBuf);
		if(*inBuf == 'R')
			R[1<<i] = 1;
		else B[1<<i] = 1;
		allR += (r[i] = readint());
		allB += (b[i] = readint());
	}
	memset(dp,-1,(1<<MaxN)*MaxN*MaxN<<2);
	dp[0][0] = 0; // nothing here
	for(int S=1; S<(1<<n); ++S){
		R[S] = R[S-(S&-S)]+R[S&-S];
		B[S] = B[S-(S&-S)]+B[S&-S];
		for(int i=0; i<n; ++i){ // last bought
			if(!(S>>i&1)) continue; // invalid
			int w = R[S-(1<<i)]; // red saving
			if(w > r[i]) w = r[i];
			int v = B[S-(1<<i)]; // value
			if(v > b[i]) v = b[i];
			int xez = R[S]*(R[S]-1)/2 // max w
				+B[S]*R[S]; // 0+1+...+R-1 + R*B 
			for(int j=xez; j>=w; --j)
				if(dp[S-(1<<i)][j-w] != -1)
					dp[S][j] = max(dp[S][j],
						dp[S-(1<<i)][j-w]+v);
		}
	}
	int ans = max(allB,allR);
	for(int i=0; i<MaxN*MaxN; ++i)
		if(dp[(1<<n)-1][i] != -1)
//			printf("dp[%d] = %d\n",i,dp[(1<<n)-1][i]),
			ans = min(ans,max(allR-i,
				allB-dp[(1<<n)-1][i]));
	printf("%d\n",ans+n);
	return 0;
}