1526.cpp

这是哈尔滨工业大学acmOJ的源代码

/*  This Code is Submitted by wywcgs for Problem 1526 on 2005-10-24 at 16:54:49 */ 
#include <cstdio>
#include <cmath>

const double CW = 4.19;
const double CI = 2.09;
const int EM = 335;
const double eps = 1e-2;

int main()
{
	double mw, mi, tw, ti;
	double Qm, tm, mm, dm;
	
	while(scanf("%lf %lf %lf %lf", &mw, &mi, &tw, &ti) == 4) {
		if(fabs(mw) < eps && fabs(mi) < eps && fabs(tw) < eps && fabs(ti) < eps) {
			return 0;
		} else {
			Qm = CW * mw * tw + CI * mi * ti;
			mm = mi + mw;
			if(Qm - EM * mi > eps) {
				tm = (Qm - EM * mi) / (CW * mm);
				printf("0.0 g of ice and %.1lf g of water at %.1lf C\n", mm, tm);
			} else if(Qm + EM * mw < eps) {
				tm = (Qm + EM * mw) / (CI * mm);
				printf("%.1lf g of ice and 0.0 g of water at %.1lf C\n", mm, tm);
			} else {
				dm = Qm / EM;
				mw += dm;
				mi -= dm;
				printf("%.1lf g of ice and %.1lf g of water at 0.0 C\n", mi, mw);
			}
		}
	}
	
	return 0;
}