雾雨小镇

 找回密码
 注册
搜索

最近看过此主题的会员

楼主: [091]导演A

[091]『无主之城』—— 第一日讨论帖(2019.11.28 9:00 —— 2019.11.29 21:00)

 关闭 [复制链接]

0

荣誉

15

硬币

167

回帖

卸妆演员

发表于 2019-11-28 20:49| 字数 61 | 显示全部楼层

对!
然后让它控制蝙蝠
放大家出去!


我继续研究它的代码了!
回复

使用道具 举报

0

荣誉

575

硬币

104

回帖

卸妆演员

发表于 2019-11-28 21:11| 字数 44 来自手机 | 显示全部楼层
[091]江雪 发表于 2019-11-28 20:49
对!
然后让它控制蝙蝠
放大家出去!

好的喵
回复

使用道具 举报

0

荣誉

485

硬币

228

回帖

卸妆演员

发表于 2019-11-28 21:12| 字数 6 来自手机 | 显示全部楼层
加油^0^~
回复

使用道具 举报

0

荣誉

15

硬币

167

回帖

卸妆演员

发表于 2019-11-28 21:50| 字数 33 | 显示全部楼层

回复

使用道具 举报

0

荣誉

0

硬币

83

回帖

卸妆演员

发表于 2019-11-28 22:02| 字数 7 来自手机 | 显示全部楼层
有没有人加速啊
回复

使用道具 举报

0

荣誉

0

硬币

83

回帖

卸妆演员

发表于 2019-11-28 22:02| 字数 8 来自手机 | 显示全部楼层
晚上真刀真枪的干
回复

使用道具 举报

0

荣誉

0

硬币

83

回帖

卸妆演员

发表于 2019-11-28 22:03| 字数 9 来自手机 | 显示全部楼层
这规则白天就是划水
回复

使用道具 举报

0

荣誉

0

硬币

83

回帖

卸妆演员

发表于 2019-11-28 22:03| 字数 3 来自手机 | 显示全部楼层
好无聊
回复

使用道具 举报

0

荣誉

0

硬币

83

回帖

卸妆演员

发表于 2019-11-28 22:03| 字数 6 来自手机 | 显示全部楼层
戴面具藏身份
回复

使用道具 举报

0

荣誉

0

硬币

83

回帖

卸妆演员

发表于 2019-11-28 22:03| 字数 0 来自手机 | 显示全部楼层
回复

使用道具 举报

0

荣誉

0

硬币

77

回帖

卸妆演员

发表于 2019-11-28 23:41| 字数 41 来自手机 | 显示全部楼层
[091]江雪 发表于 2019-11-28 15:43
Test

#include

小坏雪!
回复

使用道具 举报

0

荣誉

0

硬币

77

回帖

卸妆演员

发表于 2019-11-28 23:43| 字数 19 来自手机 | 显示全部楼层
好多秀恩爱的,我已经彻底脸盲了啊哈哈哈
回复

使用道具 举报

0

荣誉

354

硬币

34

回帖

卸妆演员

发表于 2019-11-28 23:44| 字数 22 来自手机 | 显示全部楼层
可不能加速啊 我还得研究一下 这版比我想的有趣啊
回复

使用道具 举报

0

荣誉

354

硬币

34

回帖

卸妆演员

发表于 2019-11-28 23:45| 字数 9 来自手机 | 显示全部楼层
二二  我被魔鬼迷惑了
回复

使用道具 举报

0

荣誉

15

硬币

167

回帖

卸妆演员

发表于 2019-11-28 23:47| 字数 764 | 显示全部楼层
emm
让我用一个 naive code 来算 π

#include <iostream>
#include <stdlib.h>
#include <time.h>

#ifdef _OPENMP
#include <omp.h>
#endif

using namespace std;
int main(int argc, char *argv[])

{
        clock_t start,end;
        start = clock();
        long nr_tries=100000;
        if(argc >1)
                nr_tries=atol(argv[1]);

        long nr_sucess=0;

#pragma omp parallel
        {
#pragma omp for reduction (+:nr_sucess)
        for (int try_nr=0;try_nr<nr_tries;try_nr++)
        {

         double x = ((double) rand())/((double) RAND_MAX);
         double y=((double) rand())/((double) RAND_MAX);

         if(x*x+y*y<1.0)
         {
                nr_sucess++;
         }
        }
        }

        double pi=4.0*((double)nr_sucess)/((double)nr_tries);

        cout<<pi<<" "<<nr_sucess<<" "<<nr_tries<<endl;
        end = clock();
       
        cout<<"clock_time="<<(double)(end-start)/CLOCKS_PER_SEC<<endl;
   
    return 0;
}

回复

使用道具 举报

0

荣誉

15

硬币

167

回帖

卸妆演员

发表于 2019-11-28 23:48| 字数 29 | 显示全部楼层
回复

使用道具 举报

0

荣誉

0

硬币

9

回帖

卸妆演员

发表于 2019-11-28 23:53| 字数 47 | 显示全部楼层
[091]江雪 发表于 2019-11-28 23:47
emm
让我用一个 naive code 来算 π

魔鬼
回复

使用道具 举报

0

荣誉

0

硬币

77

回帖

卸妆演员

发表于 2019-11-28 23:53| 字数 44 来自手机 | 显示全部楼层
[091]江雪 发表于 2019-11-28 23:48

看写代码这么厉害,就知道是我的小坏雪!
回复

使用道具 举报

0

荣誉

0

硬币

77

回帖

卸妆演员

发表于 2019-11-28 23:57| 字数 283 来自手机 | 显示全部楼层
[091]江雪 发表于 2019-11-28 23:47
emm
让我用一个 naive code 来算 π

ang = 0:0.01:2*pi;

r = zeros(size(ang));
N = length(ang);

r(1) = 1;

syms x y
f = (x^2-1)^3/(x^5) - sin(y)^3 * cos(y)^2;

for i = 2 : N
    f1 = subs(f, y, ang(i));
    rlt = solve(f1);
    rlt = eval(rlt);

    [~, ind] = min(abs(rlt-r(i-1)));
    r(i) = rlt(ind);
end

x = r .* cos(ang);
y = r .* sin(ang);
plot(x,y);
回复

使用道具 举报

0

荣誉

485

硬币

228

回帖

卸妆演员

发表于 2019-11-29 00:39| 字数 19 来自手机 | 显示全部楼层
哇,代码写的666,给大佬们打call
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|手机版|雾雨小镇狼人村

GMT+8, 2025-8-4 16:01 , Processed in 0.076477 second(s), 24 queries .

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表