| ||||||||||
| 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 | |||||||||
没话说了,兄弟太仗义了,谢谢In Reply To:这个是我的迭代 Posted by:frkstyc at 2005-07-06 01:41:22 > 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: Post your reply here: |
All Rights Reserved 2003-2013 Ying Fuchen,Xu Pengcheng,Xie Di
Any problem, Please Contact Administrator