2012年6月8日星期五

举手之劳

凑个热闹

有个船主,让漆工给船涂漆。漆工涂好船后,顺便将船上的漏洞补好了。过了不久,船主给漆工送了一大笔钱.漆工说:“工钱已经给过了”。船主说:“这是感谢补船漏洞的钱。”漆工说:“那是顺便补的。”船主说:“当得知我的孩子们驾船出海,我就知道他们回不来了。因为船上有漏洞,现在他们却平安归来,所以我感谢你!"

淹神的铁民是不会让自己的船上有漏洞的,所有的漏洞都应该在第一时间被修补上。臭佬就是这样一个漏洞。席恩・葛雷乔伊已经分不清楚到底谁是臭佬,是拉姆斯・波顿还是自己。他终于从史坦尼斯手里逃了出来。隐姓埋名并没有太大的意义,因为已经没有人能够认出他来了。他只是一个给船涂油漆的漆工,由于缺少几根手指,他每次涂漆的动作相当缓慢。不过涂漆这件事,本来就是个细活。慢一点到也没什么坏处。

席恩掂了掂手里的钱袋,笑了。不过他的豁牙吓了船主一跳。船主意识到自己的失态,他黝黑肥硕的脑袋上很快又恢复了优雅的微笑。席恩知道他是谁――札罗・赞旺・达梭斯,十三巨子之一,也知道他的船上为什么有个漏洞――他这么黑乎乎的家伙居然把拥有银色长发紫色眼睛的丹妮莉丝称为自己的孩子。

“铁民拿到报酬必须是铁钱”。突然,一片刀锋从札罗胖胖的前胸中穿出。阿莎从札罗的身后闪了出来,擦擦她的刀。“札罗,你已经在另一处里被风暴丹妮关在自己的金库里饿死了。再死一次也没什么了不起的。”

席恩和阿莎紧紧的拥抱在一起。一次取得十三条战船,这是一个开端。不仅仅是从札罗那里获得的,更是从乔治那里获得的。席恩向淹神祈祷,他们能避开马丁的思绪足够长的时间。

参考:

  1. 淹神
  2. 席恩・葛雷乔伊
  3. 拉姆斯・波顿
  4. 札罗・赞旺・达梭斯
  5. 阿莎・葛雷乔伊
  6. 乔治・R・R・马丁

Posted via email from goldengrape's posterous

如何使用dropbox催稿

step 0. 安装dropbox 到https://www.dropbox.com/下载并安装,至少要1.4版本以上的,这样可以支持实时更新。碰上IT水平稍低的作者,请编辑帮忙下载并安装。不过话说连科学上网都做不到,还怎么做科研,怎么写科普呢。 step 1. 请作者将写作的文稿存储在dropbox的特定目录中。 step 2. 右键点击该目录。dropbox》share this folder [image]http://img1.guokr.com/gkimage/ct/zf/l5/ctzfl5.jpg[/image] 填入编辑的dropbox账户。 step 3. 完成 作者的每一次存盘,都会在编辑的dropbox中同步更新,并且还有提醒。如果不满意一个版本,还可以回溯到上一个版本中去。万一不慎删除,短时间内也是可以恢复的。 下面补充一下写作工具: 写作什么的,最好不要使用word。要知道word是个排版工具,是面向打印机的,而不是面向写作过程的。纯文本才是最好的写作文件格式。如果要增加插图、超链接、标注什么的,推荐使用[Markdown](http://daringfireball.net/projects/markdown/dingus)。 写作软件在小众软件中介绍了许多[无干扰编辑器](www.appinn.com/tag/文本编辑器/) ,Evernote也在可选之列。更喜欢折腾的推荐使用号称[手稿级的编辑器](http://etherpad.org/)能够记录每一个字符的键入和修改过程。对于不折腾不死星人,可以考虑使用Vim或者Emacs,不过我认为那更多的是为... 我认为,最佳的写作工具其实是笔和纸。逃到远离一切显示器的地方去,只用笔在纸上写字,这才是效率最高的。笔推荐使用钢笔,摩擦阻力最小,一口气写个几千字你就知道圆珠笔的痛苦了。纸张切不可选择什么小清新的牛皮笔记本。死重不说,还要腾出一只手压着。用过的A4纸背面最佳,其次是活页本。这些可以“非线性编辑”,但要记着标记页码。 同时应该搭配一台二手扫描仪,200元以下足以。手写过的手稿统统扫描下来,然后把原件保存好。万一日后成名了可以一张一张拍卖。

Posted via email from goldengrape's posterous

2012年5月25日星期五

贡献评价函数

贡献评价函数.xlsx Download this file

现在遇到一个假想问题,如果一群人一起做一件事情,挣了钱,那么如何分呢?事先股份制确认,然后按照股权分账,显然是一种标准的商业方式。但是这种方式也不一定是好的,尤其是在创业早期,比如一次创意头脑风暴,产生了一个惊人的发明,有的人提出最初的内核,有的人进行了大幅度的扩展和提升,还有的人只做了少量修补,这个头脑风暴会议可以录制下来,大家可以播放、讨论,但是如何判定各自的贡献呢?而且如果会议之前事先确定专利的股权,头脑风暴的时候,可能就有人偷懒,没那么积极了。

于是我希望设计出一种评价方法,使每个人尽可能诚实的对自己和他人的贡献作出评价。这样组织头脑风暴会议讨论发明的时候,也可以清晰的分割知识产权了。

==问题描述==
现有若干人完成了一项任务,假定每个人都能够清楚的了解其他人所作出的贡献事件。这个条件其实大企业并不满足,员工通常并不知道老板在做什么。员工做错事顶多扣掉三个月工资,再找一份工作。老板做错事可能30年的积蓄就没了。但是微型企业应该是满足的。

每个人(i)对完成任务作出了若干贡献C(i),但C(i)不可测,只能是估计了。

现在每个人(i)对自己和其他人评价贡献。X(i,j),j=1:N。对自己贡献的评价就是X(i,i)了。
那么对C(j)的估计E(i)=E(X),比如Average(X(i,j),i=1:N)

为了要求使每个人的评价务必客观准确,所以引入罚分
P(i)=P(X,E)
比如P(i)=sqrt(sum((X(i,j)-E(j))^2,j=1:N)),相当于均方根

如何设计每个人的收益R(E(i),P(i))?
使Sum(R)=100%,
并且每个人(i)给出的X(i,j)=E(j),j=1:N,ji是纳什均衡点?

也就是说,
如果(i)对自己的评价过高,X(i,i)>>E(i)的时候,应该是R(i)去很小的值,
而X(i,i)=E(i)的时候,R(i)应该=E(i)/sum(E)

==一个方案==
现在每个人(i)对自己和其他人评价贡献。X(i,j),j=1:N
那么对C(j)的估计E0(i)=(Sum(X(i,j),i=1:N)-X(i,i))/N
修正E(i)=E0(i)/sum(E0)

为了要求使每个人的评价务必客观准确,所以引入罚分
P0(i)=sqrt(sum((X(i,j)-C(j))^2,j=1:N))
修正P1(i)=exp(-P0(i))
修正P(i)=P1(i)/sum(P1)

R0(i)=E(i)*P(i)
修正R(i)=R0(i)/sum(R)

==举个栗子==

A的评价:X(A,A)=0.30 X(A,B)=0.10 X(A,C)=0.60
B的评价:X(B,A)=0.30 X(B,B)=0.20 X(B,C)=0.50
C的评价:X(C,A)=0.10 X(C,B)=0.40 X(C,C)=0.50

对每个人的贡献估计E(i)
E(A)=0.20 E(B)=0.25 E(C)=0.55
大家评价误差产生的惩罚P(i),注意P(i)越小,惩罚越大
P(A)=0.33 P(B)=0.35 P(C)=0.33
最后结果:
R(i) R(A)=0.20 R(B)=0.26 R(C)=0.54

===再举个栗子===
如果B同学极度自恋,
A 0.30 0.10 0.60
B 0.05 0.90 0.05
C 0.30 0.10 0.60
那么
R(A)=0.32 R(B)=0.08 R(C)=0.60

Posted via email from goldengrape's posterous

2012年5月9日星期三

新五毛入推指南

1)到糊涂佛的页面https://twitter.com/#!/whotofo/following找些有趣的人交谈发些日常生活推;
2)加关注@mranti @xiaomi2020 @aiww了解时局动态;
3)加关注@baozuitun @sanwolfy学习如何攻击『反动言论』;
4)贴个头像
5)有些言论在推上使你的五毛身份非常容易被识别,典型的,包括但不限于称赞GFW;
6)请珍惜政府给予的翻墙机会,抓紧时间通过youtube以及其他不易获得的资源学习英语、学习公开课等等;
7)与推上的程序员多交流,学习一技之长,早日换个工作

拿工资做老板要求的工作是天经地义的,但工作也是学习和提高的珍贵机会。一个人从工作中获得如果不仅仅是工资,还有其他的更多的收获,就好像在进行利滚利的存钱。
做事呐,没有必要要求做到最好,那样是不现实的,但是应该每一次比上一次有所进步,或者有所优化,让收益与付出的比值不断增大。
开始的工作可以是为了养家糊口,慢慢的能力积累到一定,就要考虑自己的擅长和喜好。做自己最擅长的事情,可以在单位时间内拿到最多的钱。做自己最喜欢的事情,可以获得除了工资以外的幸福感。幸运的是,往往我们最喜欢做的事情也恰恰是最擅长做的事情。

Posted via email from goldengrape's posterous

2012年4月26日星期四

讨论-决策系统

  • 突然想起来的,有人帮忙做了么?要不我自己写?

    讨论-决策系统

    2006-11-11

    版权声明:转载时请以超链接形式标明文章原始出处和作者信息及本声明
    http://goldengrape.blogbus.com/logs/3808246.html

    这是很久以前设计的一种讨论-决策系统。那篇文章不一定能找到了,现在重新发布一下。这个创意是在GNU自由文档许可下发布。

    简单起见,以纸笔作为说明,实际操作中必须使用计算机进行。

    每一个参与者各自在一张纸上写下3个意见
    随机交换这些纸
    每人划掉一个最反对的意见,再写下一个意见。
    再次随机交换,重复这一过程。
    直到某一个意见在意见库达到一定比例后停止。
    讨论过程中,保持匿名、随机。

    可以进行一些修改,比如不止3条,为了输入方便,可以保留待删除的语句,可以直接复制/粘贴,为了达到中止条件,可以比较字符串的相似程度,如果近似程度比较高,算作一类意见。。等等

    上面这个过程,近似于一个遗传算法的过程,包含了选择、复制、交换、突变的各个步骤,系统会自进化。讨论过程中,保持匿名、随机。那么开会中常见的人与人之间的摩擦就会减少。讨论以及投票是统一的过程。也许会很有效率。

    因为没有实际测试过,所以不知道效果会怎样。比如如果重复的次数太多,特别是接近收敛的时候,会比较无聊。为了在开始快速收敛,可以调整“观察窗宽“,也就是说,一开始每个人写3条意见,交换以后,每张纸上出现6个或者9个意见,让每个人能够看到更多的意见。而到了末尾,可以只是二选一,甚至是二变一,也就是对两个不同意见选择,被否认的意见自动变成被赞同的意见。

    Technorati Tags:,,


    历史上的今天:


    收藏到:Del.icio.us



    评论

  • "如果剩下来的意见满足了最低回复数,那么就显示给用户。"

    主要是担心不够回复数,那么只能全部展现给用户,而这是我们不想看到的,因此我们可以设置最低回复数,例如2倍的参与人数量,这样可以保证不要把所有的意见都呈现出来。

    金色葡萄回复Shenxinaz说:
    期待中。
    最初的设计是作为一个病历讨论的系统的,相当于一个delphi专家系统和一个投票系统,后来某次参加创业大赛,把这个改成了管理系统中的一部分。虽然过了初赛,不过恐怕这部份没有引起什么重视。
    某日被问起产品经理是否可以兼职。又重新想了起来。
    2006-11-23 04:38:42
  • 这个设计够细了吧?

    还有什么补充的不?可以开始写软件了。呵呵

    匿名讨论系统概要设计:

    1、基于 Web的数据库

    2、发布制度,由某个人发起讨论,同时记录编号、日期、内容等信息,还有一些相关参数设置,例如,收敛的结果在那一个范围内,同时参与最低人数等。

    3、参与制度,用户可以自由选择参与讨论的内容,参与以后,就会获得一个许可。

    4、限制制度(为了防止机器人作弊,可能要加入某些限制)

    5、通知机制,作为非立刻表决的结果,需时较长,可能涉及到不同时区的工作时段,需要有通知机制,帮助大家同步。

    6、投票机制:

    6.1、第一次提交,(是否需要先确认看完了对内容的描述?)

    6.2、第一次通知,必须在满足条件(最低人数、最低回复数)下,系统自动通知(邮件或者界面方式),可以是用户登录以后。

    6.3、筛选条件:不看自己的,不看看过的人的同一条意见,如果剩下来的意见满足了最低回复数,那么就显示给用户。

    6.4、其他,根据总参与人数比例,决定显示窗口的大小。

    7、严格投票和松散投票,严格意味着必须每个人都提交了才可以交换意见,松散意味着可以随时提交。

    7.1、严格投票时,由发起讨论人决定是否把某人踢出,以防止由于他的低效率造成延误

    7.2、松散投票时,可以多次登录以获得更大的窗口,因为多次登录表示对该事件的重视,理应鼓励。

    7.3、问题,松散投票时,如何保证后来者的意见不被这个机制忽视?

    8、每次的新的意见,建议限制在30或100个字以内,但是可以添加链接到某个地方,用辅助文字说明。

    金色葡萄回复Shenxinaz说:
    "如果剩下来的意见满足了最低回复数,那么就显示给用户。"

    把所有的意见都呈现出来么?

    投票的过程是我说的那种方式么?因为相当与一个开放选项的投票,项目个数可能在开始几轮内暴涨。这个阶段更类似一个头脑风暴的讨论阶段,只呈现方案。

    很想看看实际效果如何?终于有人开始做了。太好了。 另外,很激动,很感动,多谢多谢!

    2006-11-15 13:12:35
  • 貌似挺有意思的,应该可以写个软件来做得把?

Posted via email from goldengrape's posterous

2012年4月14日星期六

转贴:医患大战之一--中国人为什么看不起病

医患大战之一--中国人为什么看不起病

阿鼠

阶级斗争的幽灵始终在中华大地游荡。医患大战只是其中精心导演的一幕。

几年来,我一直在准备一部医疗剧,希望为越演越烈,水火不容的医患关系,搭起一座沟通的桥梁。目前的医疗剧,和所有国产影视剧一样,除了拙劣的写作技术,就是医生们高大全的虚假形象。连真相的1%都没有触及。

在层层剥茧似的准备过程中,我发现,医患冲突,远远不是13.3亿百姓和190万执业医生以及233万执业助理医师之间的冲突。医患冲突的的背后,是党,政,警,司法,新闻,企业,利益集团交互纠缠的身影,甚至黑社会也占据一席之地。

最后,我颓然放弃。

在呼吸不到自由空气的环境里,所有你读到的新闻资讯,所有你了解的"真相",都是经过精心选择和过滤的,是故意让你看到的,称之为洗脑。真相的出炉和文明的传播,原本需要不同观点不同立场的激烈碰撞。可是,今天,体制内的医生全体被禁言。只要为医生说话,无论是我写的文章,还是同事在医院内部网上发表的观点,都因为涉及敏感话题迅速被删。

作为一个医生,同时作为一个挣扎在死亡线上九年的罕见病病人,作为曾经在美国生活过,又视写作为生命的人,我相信很少有人比我更有资格,站在公允的角度来写中国医疗的真相和缺陷。

1, 改革开放前的全民免费医疗:

1983年,世界卫生组织(WHO)将中国列为发展中国家全民免费医疗的典范。

1983年的中国十分贫穷。政府却在医疗方面投入巨资。组建了健全的全民免费医疗网。超过96%的中国人。享受到了免费医疗的保障。医疗保障当时由三个系统组成,劳保医疗覆盖职工及家属达2亿多人,公费医疗覆盖5千多万人口,合作医疗覆盖超过90%的行政村和农民。

另一方面,所有医院由国家财政全额拨款。尽管医生态度很差。过度医疗占医疗总费用的20%。但是百姓给了医生很高的评价。

2, 医改目的:

1984年,轰轰烈烈的国企改革和公费医疗制度改革拉开序幕。从此政府甩开一切包袱,大干快上经济,冲高GDP。同期,大批下岗职工和农村人口被抛出传统医疗体制之外。

到1998年,参与农村合作医疗的,降到农村居民的6.5%。2002年,全国城镇职工医疗保险只覆盖了9千多万人。超过88%的中国人,被国家摒弃在医疗保险体制外,看病全自费。因病返穷大批出现。

2000年,同样是WHO报告,在全世界191个国家和地区中,中国卫生服务的公正性排名倒数第四。比阿富汗还低。医患大战开始上演。

3, 医院现状:

截至2008年,我国共有医院19244家,财政拨款仅占医院收入的2%到8%,其中大多数省份的财政补助占医院收入在3%-5%之间。连退休人员工资都不够。最低的县医院,财政补助连续三年只有17.5万元,占医院收入的0.01%。

2008年,中国医院总资产10191亿元,其中负债约3250亿元。资产负债率约31.9%,目前全国90%的县医院负债经营。按83年前的标准衡量,中国已经没有一家公立医院。

以北大医院为例,为满足日益增多的病人需求,买地建第二住院部时,贷款1个多亿5年还清,每年两千多万元的还款,全部来自于对病人的收费。医院,在裹胁而来的市场经济浪潮中,在政府的数次医改中,早已失去了公益性质。留给中国医院的选择只有两个,要么靠医生卖药做检查勉强维持,要么关门大吉。

4, 中国医疗费用贵不贵?

即便所有的医院大楼,昂贵的仪器设备,以及医务人员的工资奖金都来自于病人收费。

2010年,中国卫生总费用仍然只占GDP的4.3%,而WHO的最低标准是5%。美国是16.1%。今天,美国几乎所有物价都比中国便宜很多,里维斯牌牛仔裤是中国的1/7,汽车是1/3,猪肉是1/2到1/3,别墅是中国公寓的几分之一。汽油便宜一小半而且质量好。

唯有医疗价格,美国却是中国的6-10倍。同样一个感冒,美国平均400美元。一台阑尾炎手术7000多美元。急诊什么药也没有开,做了些检查,账单寄过来是5000多美元。一台宫外孕手术,2.5万美元。远远比中国贵。

中国医疗总费用占GDP的比例,几乎全世界最低。中国的医务人员,以极其艰辛的努力和低廉的收入,中国的百姓,以最大的意志力忍受着疾病的折磨,他们共同实现了全世界最低的医疗费用,实现了达到发达国家的人均寿命。中国看病真的不贵。

5,为什么中国百姓觉得看病贵?

很多人指责,既然中国医疗费用如此低廉,为什么百姓看病要倾家荡产?为什么百姓感受到的实际情况和数据截然相反?确实。在中国,无数挣扎在疾病中的人,承受着肉体的折磨,心灵的痛苦外,还要面对经济的巨大压力,心理频临崩溃。然而,不仅仅是病人,很多医生得病后,同样面临看不起病的困难,尤其退休后的医生们。

且看一组数据。

2007年,美国政府的财政支出中,社会性支出(含教育、社会保障、公共医疗卫生、住房支出)比重高达58.4%。其中医疗支出又是重中之重。而加拿大和欧洲更高,丹麦为71.6%,德国为70.8%。

美国7800万穷人,62岁以上的老人,残疾人看病全免费。医院照常收费,只是钱由美国政府支付给医院。目前美国医疗产业的总值1.4万亿美元,超过教育和国防的总开支。

而中国,以2005年为例。当年卫生总费用为8659.91亿元。占GDP4.73%,其中政府预算投资1552.53亿元,占卫生总费用的17.93%。政府投入的卫生总费用一共包含十项,遗憾的是,去除官方的各项事业费基建费,基本医疗保险补助只占37.13亿元,卫生事业费占593.22亿元。合计占政府预算投资的40%。而在卫生事业费里,属于医院经费的又只有36.30%,加上社区服务中心和卫生院补助,也只占卫生事业费里的51%左右。

计算一下,政府投入的17%医疗支出中,用于13亿百姓的,只有40%×51%=20.4%。

卫生部前副部长殷大奎给出了相同的数据。800万党政干部花掉了医疗投入的80%。

换言之,中国百姓每花一百块钱看病,来自于政府投入的只有100×17%×20.4%=3.47块钱。职工医保是扣你自己工资,新农合国家补助每人一两百块钱。仅此而已。

2010年,撇开800万各级党政干部的医疗经费,中国政府对于13亿百姓的实际医疗投入,只占GDP的4.3%×17%×20.4%=0.149%。而同年,美国政府为医疗投入约占GDP的16%×60%=9.6%。甩掉医疗这个包袱,美国政府立刻富的流油,倒过来成为债权国。

下面是各国政府补贴占卫生总费用的比例。充分说明了为什么中国的医疗费用明明很低,百姓却看不起病。补充一下,2010年美国的人均卫生费用已达8000多美元。

国家 各国人均卫生费用 人均政府卫生支出:       政府补助占卫生总费用比例

中国: 91美元,   39美元(其中20%用于医院和病人)8.57%

美国: 6719美元,  3076美元             45.7%(加上政府减税的补贴, 实际美国政府承担了美国百姓超过60%的卫生费用)

意大利:2819美元 2172美元 77.0%

法国: 3937美元,  3139美元             79.7%

日本: 2357美元, 2242美元 95.1%

十几年来,中国的医生被杀或被重伤致残的案例,有几百个。大多数被杀医生的血被都悄无声息的和谐了。公众所见到的医生被杀,是近两年微博兴起后医生群体的努力。自哈医大一院王浩医生被杀,超过60%的网民拍手叫好,我开始查阅资料,准备写医患大战的系列文章,然而数次修改,却总是难以完成。医患冲突背后有着太多深层次的原因,无数部门都牵涉其中。

昨天,北京再次发生两起杀医生的血案。然而,一天后,相关新闻从新浪的头版里消失。无数相关微博都被删除。和以前一样,医生被杀再次被和谐了。

作为医生,我们理解,现今的中国,各行各业都充斥着仇恨的情绪,医患冲突不过是其中的一个缩影。我们也理解,生病的人尤其脆弱,他们渴望哪怕一丁点的情感温暖和经济帮助,然而就是这一丁点,在这个社会是如此奢求。很多病人举债而疾病不愈,很多病人遭到亲人的冷漠和抛弃,很多病人面对失业后的经济窘迫,很多病人看到社会大众的歧视和欺压,这一丁点温暖,在最后寄希望于医生而不得时,就化成了对医生的仇恨。

和你们每个人一样,医生只是一种工作,当我们每天要面对一百多个病人时,当我们每天超时工作透支体力脑力时,留给你们每个人的关心和照顾,和留给我们自己的一样,微乎其微。

看病难看病贵的体制,你们改变不了,我们也改变不了。情感的温暖和帮助,无论我们多么倾心给予,对于病人的巨大痛苦,都不过是杯水车薪。比医生关怀更重要的,是亲人的情感支持和心理医生的心理辅导。而最终的平静,只能来自于内心的力量,以及上帝以及神灵的信仰。

我们唯一不希望的是,你们和我们斗的你死我活时,却放过了我们共同的敌人--疾病。

根据常识,本文将很快被和谐。

欢迎扩散,转载时请注明作者,来自阿鼠医生的文章--医患大战之一。

稍后写医德(红包,回扣和态度),医术,医院收费组成。

以及中国人为什么看病难。

http://blog.sina.com.cn/s/blog_48e8d37001014ks1.html

Posted via email from goldengrape's posterous

2012年4月12日星期四

有理数的葡萄判定

倒时差,失眠中。半夜三更,思维不清楚,如果错了,见谅。
之前@双料黑客 提出过一个问题,[如何判定一个数是无理数](http://www.guokr.com/question/130064/)

今有不知何数x,注意要产生一个可能是无理数的数字,肯定是要无限位的,所以是需要一个程序来生成这个x的,这是如何产生一个无理数的问题,不在此处讨论。请参考[解读求pi的怪异代码](http://wangcong.org/articles/puzzle.html)。

不失一般性,以十进制来表示x,不失一般性,只考虑小数部分,也就是0<x<1

伪代码如下:

x(0)=floor[x*2]
  x=x*2-x(0)
  y=x
x(1)=floor[x*2]
  x=x*2-x(0)
u=x(0)*2+x(1)

Do
n++
x=y
x(0)=floor[x*2]
  x=x*2-x(0)
  y=x
x(1)=floor[x*2]
  x=x*2-x(0)
v=x(0)*2+x(1)

// 将x转化成二进制小数,取前3位,第1位和第2位组成一个数u,第2位和第3位组成一个数v,(u,v)是平面上一个点。(u,v)可能有四个(0,0), (0,1), (1,0), (1,1),其数量分别记为a,b,c,d

if u=0 and v=0 then a++
if u=0 and v=1 then b++
if u=1 and v=0 then c++
if u=1 and v=1 then d++

SD=(a-n/4)^2+(b-n/4)^2+(c-n/4)^2+(d-n/4)^2
// 如果是无理数,那么(u,v)将会等可能的分布在这四个点上,也就是说当重复次数足够多的时候,a,b,c,d是非常接近的,那么其方差Sqrt(SD)是一个非常小的数,这个数是一个小于F(n)的数字,其中F(n)是一个与重复次数有关的函数,具体是什么,我今天也不知道,恐怕要留着下次失眠的时候来算了。

if SD>F(n) then return(x是有理数), break

u=v
Loop

-EOF-

这是根据前几天看过的一个“通过Pi求Pi”的文章来写的,名字不记得了,好像原文的方法是取Pi的小数部分,每4位算一个点,比如3.14159265358....,那么(0.1415, 0.9265)就是一个点,(0.3589,0.7932)是下一个点。如果是无理数,那么应该布满单位正方形,于是u^2+v^2<=1的点的数量,应该总数的Pi/4。

其实不一定是要来求圆,只要是这些点均匀的布满单位正方形之内就可以,那么随便找个图形画上,只要图形内的点数量的比例=图形面积即可。

而且我不喜欢用Pi,这种方法有两个问题,第一,用一个无理数去判定一个无理数,这数值计算烦死了。第二,即使是4位小数,也是平面上离散的有限个点,跟Pi总是有很大差距的。

所以,既然是平面上的离散点,那么干脆取出最小的值,也就是4个点来做就好了,如果只有4个数的话,还算什么单位圆,直接求方差就够了。然后,我不喜欢间断的求点,我觉得那样会漏掉东西,(0.1415, 0.9265),(0.3589,0.7932)我觉得不好,我更喜欢(0.1415,0.4159),(0.1592,0.5926)这样的点阵。

同样无聊的人们,谁愿意来求解一下最大能够满足F(n)判别作用的表达式呢?

Posted via email from goldengrape's posterous