首页 > 代码库 > POJ 3296 Rinse
POJ 3296 Rinse
果然是杂题,都没有人做啊,这解题报告独一份~~
题目大意:
Granny有一个罐子里有Vw的酒,她有Vb的雨水来冲这个罐子,由于表面张力的原因当罐子往外倒液体时,会有Vr的液体倒不出来。罐子的容量是Vc。
她最多有K次拿水冲洗罐子的机会,每一次它可以取一些雨水倒入罐子充分混合之后在将罐子里的液体倒出来。
问怎样冲洗可以使罐子内残留的酒的数量最少。
解题思路:
对于操作来说,第一次是最重要的,后来k-1次倒入水量是相同的。对于第一次来说可以用三分法来求极值。
下面是代码:
#include <set> #include <map> #include <queue> #include <math.h> #include <vector> #include <string> #include <stdio.h> #include <string.h> #include <stdlib.h> #include <iostream> #include <algorithm> #define eps 1e-8 #define pi acos(-1.0) #define inf 107374182 #define inf64 1152921504606846976 #define lc l,m,tr<<1 #define rc m + 1,r,tr<<1|1 #define iabs(x) ((x) > 0 ? (x) : -(x)) #define clear1(A, X, SIZE) memset(A, X, sizeof(A[0]) * (SIZE)) #define clearall(A, X) memset(A, X, sizeof(A)) #define memcopy1(A , X, SIZE) memcpy(A , X ,sizeof(X[0])*(SIZE)) #define memcopyall(A, X) memcpy(A , X ,sizeof(X)) #define max( x, y ) ( ((x) > (y)) ? (x) : (y) ) #define min( x, y ) ( ((x) < (y)) ? (x) : (y) ) using namespace std; double vb,vw, vr, vc; int k; double does(double first) { double temp=min(vr,vw+first)*(vw/(vw+first)); double v=vb-first,last=vw+first; if(k<=1)return temp; v/=(k-1); if(vr+v>vc)v=vc-vr; for(int i=1; i<k; i++) { temp=min(vr,last+v)*(temp/(min(vr,last)+v)); last=min(vr,last+v); } return temp; } int main() { double l,r,lmid,rmid,leftover,ans1,ans2,ansa; while(scanf("%d",&k),k) { scanf("%lf%lf%lf%lf",&vb,&vw,&vr,&vc); if(vb + vw < vr) { puts("0"); continue; } l=0.0; r=vb; leftover=vb-(k-1)*(vc-vr); if(l<leftover)l=leftover; if(r+vw>vc)r=vc-vw; while(r-l>eps) { lmid=(l+r)/2; rmid=(lmid+r)/2; if(does(lmid)<=does(rmid))r=rmid; else l=lmid; } printf("%d %.2f",k,r); if(k>1) { ansa=(vb-r)/(k-1); ansa=min(vc-vr,ansa); for(int i=1; i<k; i++) { printf(" %.2f",ansa); } } puts(""); } return 0; }
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。