CodeForces 799D :Field expansion DP

传送门

题目描述

小L有一块大小为h × w的农田,但是这个农田太小了。幸运的是有 n 个延伸道具。第 i 个道具可以让农田的长或宽乘上 ai。每个延伸道具最多只能用一次。

现在他想在农田里面放一个大小为 a × b 的东西(可以横放也可以竖放),问最少用几个延伸道具。

分析

需要分析一波,我们最多使用34个道具,因为2^17已经大于1e5了
然后我们肯定是尽量选大的,然后用f[i][j]表示在使用第i个道具的情况下一条边长为j的时候,另一条边的最大值是多少,然后进行状态转移即可

最后需要注意一下将两条边互换一下,再求一次

代码

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <queue>
#include <cstring>
#define debug(x) cout<<#x<<":"<<x<<endl;
#define _CRT_SECURE_NO_WARNINGS
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> PII;
const int INF = 0x3f3f3f3f;
const int N = 1e5 + 10;
ll a,b,x,y;
int n;
ll w[N];
ll f[40][N];
inline char nc() {
    
    
	static char buf[1000000], *p1 = buf, *p2 = buf;
	return p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 1000000, stdin), p1 == p2) ? EOF : *p1++;
}
template <typename _Tp> inline void read(_Tp&sum) {
    
    
	char ch = nc(); sum = 0;
	while (!(ch >= '0'&&ch <= '9')) ch = nc();
	while (ch >= '0'&&ch <= '9') sum = (sum << 3) + (sum << 1) + (ch - 48), ch = nc();
}

bool cmp(ll A,ll B){
    
    
	return A > B;
}


int main(){
    
    
	read(a),read(b),read(x),read(y);
    read(n);
    for(int i = 1;i <= n;i++) read(w[i]);
    if((a <= x && b <= y) ||(a <= y && b <= x)){
    
    
    	puts("0");
    	return 0;
    }
    sort(w + 1,w + 1 + n,cmp);
    memset(f,-1,sizeof f);
    f[0][min(a,x)] = min(b,y);
    for(int i = 1;i <= min(35,n);i++)
    	for(ll j = 1;j <= a;j++){
    
    
    		if(f[i - 1][j] == -1) continue;
    		ll p1 = min(a,j * w[i]);
    		ll p2 = min(b,f[i - 1][j] * w[i]);
    		f[i][j] = max(f[i][j],p2);
    		f[i][p1] = max(f[i][p1],f[i - 1][j]);
    	}
    int ans = INF;
    for(int i = 1;i <= 35;i++){
    
    
    	if(f[i][a] >= b){
    
    
    		ans = i;
    		break;
    	}
    }
    swap(a,b);
    memset(f,-1,sizeof f);
    f[0][min(a,x)] = min(b,y);
    for(int i = 1;i <= min(35,n);i++)
    	for(ll j = 1;j <= a;j++){
    
    
    		if(f[i - 1][j] == -1) continue;
    		ll p1 = min(a,j * w[i]);
    		ll p2 = min(b,f[i - 1][j] * w[i]);
    		f[i][j] = max(f[i][j],p2);
    		f[i][p1] = max(f[i][p1],f[i - 1][j]);
    	}
    for(int i = 1;i <= 35;i++){
    
    
    	if(f[i][a] >= b){
    
    
    		ans = min(ans,i);
    		break;
    	}
    }
    if(ans == INF) ans = -1;
    printf("%d\n",ans);
    return 0;
}


/**3 3 2 4 4
2 5 4 10
*/

/**
*  ┏┓   ┏┓+ +
* ┏┛┻━━━┛┻┓ + +
* ┃       ┃
* ┃   ━   ┃ ++ + + +
*  ████━████+
*  ◥██◤ ◥██◤ +
* ┃   ┻   ┃
* ┃       ┃ + +
* ┗━┓   ┏━┛
*   ┃   ┃ + + + +Code is far away from  
*   ┃   ┃ + bug with the animal protecting
*   ┃    ┗━━━┓ 神兽保佑,代码无bug 
*   ┃        ┣┓
*    ┃        ┏┛
*     ┗┓┓┏━┳┓┏┛ + + + +
*    ┃┫┫ ┃┫┫
*    ┗┻┛ ┗┻┛+ + + +
*/


猜你喜欢

转载自blog.csdn.net/tlyzxc/article/details/113066023
今日推荐