【9309】求Y=X1/3

Time Limit: 1 second

Memory Limit: 2 MB

问题描述

求Y=X1/3次方的值。X由键盘输入(x不等于0，在整型范围内）。利用下列迭代公式计算：

yn + 1=2/3*yn+x/(3*yn2)（这里的2是平方），初始值y0=x,误差要求小于ε=10-4次方。

Input

输入只有一行，一个整数x

Output

输出只有一行，输出x和y的值，x按照整数输出，不含小数部分，y的值紧跟其后，加"：3x="字样，按照实数形式输出。

Sample Input1

Sample Output2

8：3X=2.000000011E+00 （冒号后面无空格）

Sample Input2

Sample Output2

error！

【题解】

这题用C++根本写不出来。。

用的pascal。就是一个迭代过程。

把前面算出来的值an+1赋给原来参与运算的an.

然后再用新的an来算an+1即可。

精度就是an+1-an的值。让他的绝对值小于1e-4就可以了

记得x==0时判错。

【代码1】pascal

program cheng;

var

	x:longint;

	y,y0:real;

begin

	readln(x);

	if x = 0 then

		begin

			writeln('error!');

			halt;

		end;

	y0:=x;

	repeat

		y:=y0;

		y0:=2/3*y0+x/(3*y0*y0);

	until abs(y-y0) < 1e-4;

	writeln(x,':3x=',y0);

end.

【代码2】C++版（但是不能过测评）

#include <cstdio>

int t;

double temp,temp2,x;

int main()

{

	scanf("%d",&t);

	x = t;

	temp = x;

	temp2 = ((2.0/3.0)*temp) + (x/(3.0*temp*temp));

	double tt = temp2-temp;

	if (tt < 0) tt=-tt;

	while (tt >= 1e-4)

		{

			temp = temp2;

			temp2 = (2.0/3.0)*temp + x/(3.0*temp*temp);

			tt = temp2-temp;

			if (tt <0 ) tt=-tt;

		}

	printf("%d:3X=%E",t,temp2);

	return 0;

}

巴特西

【9309】求Y=X1/3

Input

Output

Sample Input1

Sample Output2

Sample Input2

Sample Output2

最新文章

热门文章