技术宅社区 | Z站»家 › 技术区 › 程序设计 › 技术分享 › 踢馆帝原创小测试，第二题！

查看: 822|回复: 10

[普通教程] 踢馆帝原创小测试，第二题！

[复制链接]

踢馆帝

签到天数: 1 天

连续签到: 1 天

[LV.1]初来乍到

9 主题	11 好友	1118 积分

Continue

积分: 1118

发表于 2012-4-12 19:50:48 | 显示全部楼层 |阅读模式

╱人◕‿‿◕人╲定下契约

您需要登录才可以下载或查看，没有账号？╱人◕‿‿◕人╲订下契约（注册新用户）

本帖最后由踢馆帝于 2012-4-12 20:22 编辑

开始吧！计算结果要发出来哦
简述你的算法！

测试, 原创

评分

参与人数 1	宅币 +10	收起理由
风音洛洛	+ 10

查看全部评分

相关帖子

签名被小宅喵吞掉了~~~~(>_<)~~~~

使用道具举报

风音洛洛

该用户从未签到

81 主题	92 好友	2万积分

第一章

去哪里找比我更懒得

积分: 28955

发表于 2012-4-12 19:52:19 | 显示全部楼层

竟然不能用sqrt。。。

喂，不要戳奇怪的地方啦~要戳这里！

回复支持反对

使用道具举报

ZeroWinte

该用户从未签到

5 主题	7 好友	1万积分

第一章

积分: 10526

发表于 2012-4-12 20:15:51 | 显示全部楼层

本帖最后由 ZeroWinte 于 2012-4-12 20:33 编辑

using System;
namespace ConsoleApplication2
{
class Program
{
static void Main(string[] args)
{
double e = 0.0;
e=Math.Pow((double)(1 + 1 / 100000000.0), 100000000.0);
Console.WriteLine("e={0}",e);
double t = 0.0;
t = 1.3;
double a = 0.0;
while (Math.Abs(t - a) > 0.0001)
{
a = t;
t=t + 0.5 * (e / t - t);
}
Console.WriteLine("sqrt(e)={0}", t);
Console.ReadKey();
}
}
}

复制代码

e=lim0->∞(1+1/x)^x 精度要求也不是很高所以取个大一点的数就可以了。
开方：
计算公式
X(n+1)=Xn+( (A/Xn) - Xn ) * 1/2

A是要开方的数。
1，给x一个初始值，（我用变量t代x）
2，记录x的值，（变量a）
3，将x代入公式计算，并把结果赋给x
4，比较x和a（差值越小就越接近精确值）
5，若差值较大，返回第2步
6，输出结果x
（一开始代入错误的值也没关系）。

其实我公式是从度娘来的：http://baike.baidu.com/view/1188043.htm

签名被小宅喵吞掉了~~~~(>_<)~~~~

回复支持反对

使用道具举报

踢馆帝

签到天数: 1 天

连续签到: 1 天

[LV.1]初来乍到

9 主题	11 好友	1118 积分

Continue

积分: 1118

楼主| 发表于 2012-4-12 20:22:48 | 显示全部楼层

ZeroWinte 发表于 2012-4-12 20:15

登录/注册后可看大图

简述一下算法，方便大家交流

签名被小宅喵吞掉了~~~~(>_<)~~~~

回复支持反对

使用道具举报

秋声赋

该用户从未签到

28 主题	103 好友	9783 积分

技术宅认证程序员

苦逼,受剥削人民

积分: 9783

发表于 2012-4-12 20:26:12 | 显示全部楼层

本帖最后由秋声赋于 2012-4-12 20:33 编辑

C#

e=1+1+1/2!+1/3!+…+1/n!+
开方是网上看的...

private void button2_Click(object sender, EventArgs ef)
{
double e = 1f;
int i = 1;
do
{
double temp=1;
for (int j = 1; j <= i; j++)
{
temp *= j;
}
e += 1.0/temp;
i++;
}
while(!(i==10));
double x1 = e/2;
double x2 = (x1+e/x1)/2;
do
{
x1 = x2;
x2 = (x1 + e / x1) / 2;
} while (Math.Abs(x1 - x2) >= 0.01);
MessageBox.Show(x1.ToString());
}

复制代码

签名被小宅喵吞掉了~~~~(>_<)~~~~

回复支持 1 反对 0

使用道具举报

Flagalac

该用户从未签到

8 主题	11 好友	3239 积分

序章

积分: 3239

发表于 2012-4-12 21:10:07 | 显示全部楼层

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
namespace ConsoleApplication3
{
class Program
{
static void Main(string[] args)
{
float es = 1, b = 1, c = 1,d=0,k=0,s=0;
int i;
for (i = 0; i < 50;i++ )
{
c = c * b;
es = es + (1 / c);
b++;
}
Console.WriteLine("e=" + es.ToString() + "\n");
if (es > 1)
{ k = 1; d = 0.00001F; }
if (es > 0 & es<1)
{ k = 1; d = -0.00001F; }
if (es == 0)
{ k = 0; d = 0; }
if (es < 0 & es>-1)
{ k = -1; d = 0.00001F; }
if (es < -1)
{ k = -1; d = -0.00001F; }
s = k;
for (i = 0; i < 100000; i++)
{
s = s + d;
if (Math.Abs(s * s - es) < 0.001)
{break;}
}
Console.WriteLine("x=" + s.ToString() + "\n");
Console.ReadLine();
}
}
}

复制代码

签名被小宅喵吞掉了~~~~(>_<)~~~~

回复支持反对

使用道具举报

LonghronShen

该用户从未签到

10 主题	5 好友	1441 积分

Continue

积分: 1441

发表于 2012-4-12 23:53:58 | 显示全部楼层

本帖最后由 LonghronShen 于 2012-4-13 00:06 编辑

// Learn more about F# at http://fsharp.net

#light

open System

let sqrtfloat( x:double ):double =

let mutable g0 = 0.0

let mutable g1 = 0.0

if x = 0.0 then

0.0

else

g0 <- x / 2.0

g1 <- (g0 + x / g0) / 2.0

while Math.Abs(g1 - g0) > 0.01 do

g0 <- g1

g1 <- (g0 + (x / g0)) / 2.0

let e =

let mutable sum:double = 0.0

let multiplicativeInverse(x:double) = 1.0 / x

let rec factorial(n:double):double =

if n = 0.0 then

1.0

else

n * factorial(n - 1.0)

for i = 0 to 10 do

sum <- sum + ((double i) |> factorial |>multiplicativeInverse)

sum

Console.WriteLine(e)

Console.WriteLine(sqrtfloat(e))

let ignore = System.Console.ReadKey(false)

签名被小宅喵吞掉了~~~~(>_<)~~~~

回复支持反对

使用道具举报

romanosovsky

该用户从未签到

1 主题	11 好友	1715 积分

Continue

积分: 1715

发表于 2012-4-15 00:52:26 | 显示全部楼层

本帖最后由 romanosovsky 于 2012-4-15 18:01 编辑

我也来一发。Ruby代码

pow=sqrt_e=1.0
e=(1..10).inject{|sum, x|sum + (pow/=x)}+1
sqrt_e=(e/sqrt_e+sqrt_e)/2 until (sqrt_e*sqrt_e-e).abs<Float::EPSILON
p e, sqrt_e

复制代码

泰勒公式求e，牛顿迭代开根号。没有用(1+1/x)^x的办法来求e，因为这样不太准确而且其实内部应该还是调用了exp来实现
结果。

note: 原来的第2行：e=(1..10).map{|x|pow *= 1.0/x}.inject{|sum, x|sum + x}+1 多了一次迭代。已修改（@花晴兰）

签名被小宅喵吞掉了~~~~(>_<)~~~~

回复支持反对

使用道具举报

秋声赋

该用户从未签到

28 主题	103 好友	9783 积分

技术宅认证程序员

苦逼,受剥削人民

积分: 9783

发表于 2012-4-15 18:12:44 | 显示全部楼层

romanosovsky 发表于 2012-4-15 00:52

登录/注册后可看大图

我也来一发。Ruby代码泰勒公式求e，牛顿迭代开根号。没有用(1+1/x)^x的办法来求e，因为这样不太准确而且其 ...

这Ruby的代码也略短了点吧。。

签名被小宅喵吞掉了~~~~(>_<)~~~~

回复支持反对

使用道具举报

LonghronShen

该用户从未签到

10 主题	5 好友	1441 积分

Continue

积分: 1441

发表于 2012-4-15 20:47:00 | 显示全部楼层

秋声赋发表于 2012-4-15 18:12

登录/注册后可看大图

这Ruby的代码也略短了点吧。。

拼代码长度，Ruby笑而不语……

签名被小宅喵吞掉了~~~~(>_<)~~~~

回复支持反对

使用道具举报

johnnyworm

该用户从未签到

0 主题	1 好友	3402 积分

序章

积分: 3402

发表于 2012-5-14 08:47:09 | 显示全部楼层

本帖最后由 johnnyworm 于 2012-5-14 09:32 编辑

LonghronShen 发表于 2012-4-15 20:47

登录/注册后可看大图

拼代码长度，Ruby笑而不语……

膜拜下Ruby……看到如此短的代码实在是忍不住想比一下……奈何熟悉的java完全不在一个档次上= =只能祭出才学没多久的python了……
计算方法是一样的……

[mw_shl_code=python,true]f = lambda x: x and x * f(x - 1) or 1
e = reduce(lambda x,y:x + 1.0/f(y), range(1,10)) + 1
sqrt_e = lambda x,y: abs((x + e/x)/2-x)>y and sqrt_e((x + e/x)/2, y) or x
print e, sqrt_e(1.0, 0.0001)[/mw_shl_code]

结果：2.71828152557 1.64876387193
懒得弄图了……

附：java版
[mw_shl_code=java,true]public class Test{
  public static void main(String[] args){
double e = 1.0, sqrt_e = 1.0;
int  f = 1;
for (int i = 1; i <= 10; i++){
   f *= i;
   e += 1.0 / f;
}
while(Math.abs(sqrt_e-(sqrt_e=(sqrt_e+e/sqrt_e)/2))>0.0001){
}
System.out.println(e+" "+ sqrt_e);
  }
}[/mw_shl_code]

结果：2.7182818011463845 1.648721262969893

再附：python看上去比较正常版

[mw_shl_code=python,true]e = sqrt_e = f = 1.0
for i in range(1,10):
  f *= i
  e += 1.0/f
while abs(sqrt_e - (sqrt_e + e/sqrt_e)/2)>0.0001:
  sqrt_e = (sqrt_e + e/sqrt_e)/2
print e, sqrt_e[/mw_shl_code]

签名被小宅喵吞掉了~~~~(>_<)~~~~

回复支持反对

使用道具举报

返回列表

发新帖

合作站点账号登陆

[普通教程] 踢馆帝原创小测试，第二题！

╱人◕‿‿◕人╲定下契约

评分

相关帖子