Detail of message

Online Judge

Problem Set

Authors

Online Contests

User

Web Board
Home Page
F.A.Qs
Statistical Charts

Current Contest
Past Contests
Scheduled Contests
Award Contest

这个是我的迭代

Posted by frkstyc at 2005-07-06 01:41:22 on Problem 1905
In Reply To:谢谢styc，凭感觉写了牛顿法，感觉自己太差了，用什么方法都WA，郁闷。 Posted by:sunmoonstar_love at 2005-07-06 01:39:30

		lp = (1 + n * c) * l;
		r = lp * lp * 0.25 / sqrt(lp * lp - l * l);
		do
		{
			double temp = lp / (2 * r);
			rp = r;
			r += r * (l - 2 * r * sin(temp)) / (l - lp * cos(temp));
		}
		while(fabs(r - rp) > r * FLT_EPSILON);
		printf("%.3lf\n", r - sqrt(r * r - l * l * 0.25));

Followed by:

没话说了，兄弟太仗义了，谢谢 (309) sunmoonstar_love 2005-07-06 01:48:46
是不是应该对r迭代而不是alpha? (0) HatTricker 2006-08-28 10:57:04
请问那个r = lp * lp * 0.25 / sqrt(lp * lp - l * l);中的r是怎么确定的呢？这个式子不太懂。 (311) miaomiaomiaomiao 2007-08-25 20:30:30

Post your reply here:

Home Page Go Back To top