2020 ICPC·小米邀请赛 决赛 J. Rikka with Book(状压dp)

2023-12-17 08:28

本文主要是介绍2020 ICPC·小米邀请赛 决赛 J. Rikka with Book(状压dp),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题目

登录—专业IT笔试面试备考平台_牛客网

n(n<=20)本书,放在桌子上,

第i本书的可以看成是li(li<=1e3)*1*1的物体,其中长为li,宽为1,高为1,

质量均匀分布,且为wi(wi<=1e3)

求n本书摞在一起,使得任意一本书都不掉下桌子时,书能伸出桌沿的长度的最大值是多少

思路来源

官方题解&申老师

题解

放的话肯定是从上往下放,这样已经放上去的可以看成是一个物体,

并且当b物品摞在a物品之上时,一定是把b物品的重心放到a物品的边沿上,

好比把a物品当成桌子,一定是放到桌沿上,

再将a和b看成同一个物品时,一定是放到下一个物品的边沿上,

一旦一个物体的质量和重心的位置确定了,这个物品的其他属性就无关紧要了,从而无后效性

所以状压每次往下垫的书是哪一本,确定放的顺序,关注的是伸出去的最大值

往下垫的书的重心位于l/2处,质量为a;

上面的书看成一体时,重心位于边沿,质量为b

那么新的重心,根据杠杆原理,位于距边沿a/(a+b)的位置,记add=a/(a+b)

记原来的伸出去的最大值为x,则新的最大值为x+add,

此外,可以旋转一下整个物体,使整个物体的重心仍落在边沿上不落下去,

但是伸出边沿的是往下垫的书的另外半边,也就是l-add这半边,二者取max即可

所以,如果最优解是第i本书伸的最远,最上面的书是1,最下面的书是n,

一定是对于j∈[1,i-1]来说,把[1,j]看成一体时,[1,j]的重心压在j+1的左边沿,

对于j∈[i+1,n]来说,将[1,j-1]看成一体时,[1,j-1]的重心压在j的右边沿

每次枚举的时候,旋转or不旋转二选一都试一下,显然可以覆盖这种情况

代码1

维护的是长度l、到左边沿的距离p、整体的质量w

// Problem: Rikka with Book
// Contest: NowCoder
// URL: https://ac.nowcoder.com/acm/contest/9328/J?&headNav=acm
// Memory Limit: 1048576 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(int i=(a);i<=(b);++i)
#define per(i,a,b) for(int i=(a);i>=(b);--i)
typedef long long ll;
typedef double db;
typedef pair<ll,int> P;
#define fi first
#define se second
#define pb push_back
#define dbg(x) cerr<<(#x)<<":"<<x<<" ";
#define dbg2(x) cerr<<(#x)<<":"<<x<<endl;
#define SZ(a) (int)(a.size())
#define sci(a) scanf("%d",&(a))
#define scll(a) scanf("%lld",&(a))
#define pt(a) printf("%d",a);
#define pte(a) printf("%d\n",a)
#define ptlle(a) printf("%lld\n",a)
#define debug(...) fprintf(stderr, __VA_ARGS__)
const int N=20,M=(1<<20)+5;
int n;
db dp[M];
int lb(int x){return x&(-x);
}
struct node{db l,p;//长度 离左侧距离int w;//质量node(){l=0;p=1e7;w=0;}db dis(){return l-p;}void show(int i=-1){printf("i:%d l:%lf p:%lf w:%d\n",i,l,p,w);}
}e[M];
//b放在a上
bool operator>(node a,node b){return a.dis()>b.dis();
}
node mer(node a,node b){db x=a.l-a.p;node c;db B=a.w*x/(a.w+b.w);if(b.p>a.l){//b更左c.l=b.l;c.p=b.p-B;//b.l-b.p+B}else{//a更左c.l=a.l+b.l-b.p;c.p=a.l-B;//b.l-b.p+B}c.w=a.w+b.w;if(c.p>c.l-c.p)c.p=c.l-c.p;//puts("");//a.show();b.show();c.show();//puts("");return c;
}
void sol(){sci(n);rep(i,0,n-1){int x=1<<i;scanf("%lf",&e[x].l);e[x].p=e[x].l/2.0;}rep(i,0,n-1){int x=1<<i;scanf("%d",&e[x].w);}int up=(1<<n)-1;rep(i,1,up){if(lb(i)==i)continue;//printf("i:%d\n",i);rep(j,0,n-1){if(!(i>>j&1))continue;int v=1<<j,oth=i^v;node w=mer(e[v],e[oth]);//只枚举最底下那个是什么if(w>e[i])e[i]=w;//e[oth].p=e[oth].l-e[oth].p;//w=mer(e[v],e[oth]);//if(w>e[i])e[i]=w;//w.show();}//if(e[i].p>e[i].l-e[i].p)e[i].p=e[i].l-e[i].p;//b[1].p=b[1].l-b[1].p;//e[i].show(i);}printf("%.10lf\n",e[up].dis());
}
int main(){sol();return 0;
}

代码2(三个顶俩代码)

发现无需维护长度l和距一端的位置p,只维护右半边伸出去的最大值即可

每次尝试一下翻或不翻

#include <bits/stdc++.h>using namespace std;
int n,l[30],w[30],x[1100000];
double f[1100000];int main()
{int i,j;scanf("%d",&n);for(i=0;i<n;i++)scanf("%d",&l[i]);for(i=0;i<n;i++)scanf("%d",&w[i]);for(i=1;i<(1<<n);i++){for(j=0;j<n;j++)if(i&(1<<j))break;x[i]=x[i-(1<<j)]+w[j];}for(i=1;i<(1<<n);i++){for(j=0;j<n;j++)if(i&(1<<j))f[i]=max(f[i],max(f[i^(1<<j)]+double(0.5*w[j]*l[j])/x[i],l[j]-double(0.5*w[j]*l[j])/x[i]));}printf("%.12lf\n",f[(1<<n)-1]);return 0;
}

这篇关于2020 ICPC·小米邀请赛 决赛 J. Rikka with Book(状压dp)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/503728

相关文章

百度/小米/滴滴/京东,中台架构比较

小米中台建设实践 01 小米的三大中台建设:业务+数据+技术 业务中台--从业务说起 在中台建设中,需要规范化的服务接口、一致整合化的数据、容器化的技术组件以及弹性的基础设施。并结合业务情况,判定是否真的需要中台。 小米参考了业界优秀的案例包括移动中台、数据中台、业务中台、技术中台等,再结合其业务发展历程及业务现状,整理了中台架构的核心方法论,一是企业如何共享服务,二是如何为业务提供便利。

hdu4826(三维DP)

这是一个百度之星的资格赛第四题 题目链接:http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1004&cid=500 题意:从左上角的点到右上角的点,每个点只能走一遍,走的方向有三个:向上,向下,向右,求最大值。 咋一看像搜索题,先暴搜,TLE,然后剪枝,还是TLE.然后我就改方法,用DP来做,这题和普通dp相比,多个个向上

hdu1011(背包树形DP)

没有完全理解这题, m个人,攻打一个map,map的入口是1,在攻打某个结点之前要先攻打其他一个结点 dp[i][j]表示m个人攻打以第i个结点为根节点的子树得到的最优解 状态转移dp[i][ j ] = max(dp[i][j], dp[i][k]+dp[t][j-k]),其中t是i结点的子节点 代码如下: #include<iostream>#include<algorithm

hdu4865(概率DP)

题意:已知前一天和今天的天气概率,某天的天气概率和叶子的潮湿程度的概率,n天叶子的湿度,求n天最有可能的天气情况。 思路:概率DP,dp[i][j]表示第i天天气为j的概率,状态转移如下:dp[i][j] = max(dp[i][j, dp[i-1][k]*table2[k][j]*table1[j][col] )  代码如下: #include <stdio.h>#include

usaco 1.3 Mixing Milk (结构体排序 qsort) and hdu 2020(sort)

到了这题学会了结构体排序 于是回去修改了 1.2 milking cows 的算法~ 结构体排序核心: 1.结构体定义 struct Milk{int price;int milks;}milk[5000]; 2.自定义的比较函数,若返回值为正,qsort 函数判定a>b ;为负,a<b;为0,a==b; int milkcmp(const void *va,c

usaco 1.1 Broken Necklace(DP)

直接上代码 接触的第一道dp ps.大概的思路就是 先从左往右用一个数组在每个点记下蓝或黑的个数 再从右到左算一遍 最后取出最大的即可 核心语句在于: 如果 str[i] = 'r'  ,   rl[i]=rl[i-1]+1, bl[i]=0 如果 str[i] = 'b' ,  bl[i]=bl[i-1]+1, rl[i]=0 如果 str[i] = 'w',  bl[i]=b

uva 10154 DP 叠乌龟

题意: 给你几只乌龟,每只乌龟有自身的重量和力量。 每只乌龟的力量可以承受自身体重和在其上的几只乌龟的体重和内。 问最多能叠放几只乌龟。 解析: 先将乌龟按力量从小到大排列。 然后dp的时候从前往后叠,状态转移方程: dp[i][j] = dp[i - 1][j];if (dp[i - 1][j - 1] != inf && dp[i - 1][j - 1] <= t[i]

uva 10118 dP

题意: 给4列篮子,每次从某一列开始无放回拿蜡烛放入篮子里,并且篮子最多只能放5支蜡烛,数字代表蜡烛的颜色。 当拿出当前颜色的蜡烛在篮子里存在时,猪脚可以把蜡烛带回家。 问最多拿多少只蜡烛。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cs

uva 10069 DP + 大数加法

代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <queue>#include <map>#include <cl

uva 10029 HASH + DP

题意: 给一个字典,里面有好多单词。单词可以由增加、删除、变换,变成另一个单词,问能变换的最长单词长度。 解析: HASH+dp 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc