24小时热门版块排行榜    

查看: 1698  |  回复: 16
本帖产生 4 个 程序强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

holmescn

金虫 (正式写手)

[交流] Euler 工程 第三十五题:循环质数 已有4人参与

197 这个质数很特别,因为1, 9, 7这三个数的循环排列也是质数。(197,719,971)

100以下有这样性质的质数有13个:2, 3, 5, 7, 11, 13, 17, 31, 37, 71, 73, 79, 97.

那么100百万以下这样的数有多少个呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
xzhdty(金币+1): 欢迎常来程序语言版讨论 2011-07-07 19:51:23
余泽成(金币+2, 程序强帖+1): 很全面! 2011-07-11 20:34:48
引用回帖:
Originally posted by libralibra at 2011-07-07 14:50:40:
多谢指教哦,我刚开始看c++,c也忘记的差不多了,还得多学习.感觉c和c++的那些个函数,库晕晕乎乎分不清楚,

那些库很容易分的啦,规则有三:
C++标准库都是不带.h的,如
C++继承C标准库是带C头的,比较容易弄混的就是这两个了:,前面的那个是c版本的,后面那个是string类
C++的头文件可以是hpp,据说是h+cpp的意思,这个多半是模板,因为模板需要编译时实例化,实现代码也在头文件里面
C和C++的操作在某些方面不同,主要有:
a.IO库是不能混合用的,他们之间不存在共享或者同步关系,scanf/printf和istream/ostream要分开使用
b.C没有名字空间,而C++有,于是C里面的一些struct被封装到了std里面,比如struct tm,C++里面它是一个类,并且引用它需要打开名字空间std
c.C的字符串处理统统集中在了string.h里面,字符串是char*和const char*,字符串的IO应该使用sscanf/sprintf,而C++的字符串处理集中在了string类中,这意味着C++使用字符串都应该用string而非char *,字符串的IO则应该使用stringstream,不过上述做法在某些时候有性能问题,操作起来很费事又挨批。
d.C没有而C++有的一个就是容器库和算法库,C的数据法宝是指针+结构体+数组,C++的数据法宝则是迭代器+类+容器,C++的容器库跟算法库对接很麻烦,迭代器固然是好东西,别人实现的拿过来就用着爽,自己实现接口就想骂人,哎……
再多的也说不上了,反正C++用起来就是很麻烦,天生麻烦,纯种没杂种优势啊,看起来好看,难养活。
漩涡的中心有一块空地,空空的。
16楼2011-07-07 16:09:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

huycwork

金虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-07 14:43:38
这个题目带有强烈的多线程暗示呀~
漩涡的中心有一块空地,空空的。
2楼2011-07-06 23:31:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-07 14:43:45
matlab使用了primes函数函数不快哦
CODE:
% Elapsed time is 49.161032 seconds.
% ans =
%     55
function result = euler35()
tic;
pm = primes(1000000);
result = 0;
for i=1:length(pm)
    x = pm(i);
    si = num2str(pm(i));
    for j=2:length(si)
        x(end+1) = str2double([si(j:end),si(1:j-1)]);
    end
    x = unique(x);
    if any(~ismember(x,pm)) % use any() is twice faster than all()
        continue;
    end
    result = result+1;
end
toc;
end

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
3楼2011-07-07 00:11:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-07 14:43:53
余泽成(金币+2, 程序强帖+1): 鼓励交流! 2011-07-11 20:35:50
c++的效率还是高
不过初学c++,谁知道还有什么数字,字符串转换的高效方法吗?
第二个子函数中计算rotation number的办法好丑陋
CODE:
#include
using namespace std;

// 素数筛: 返回0-stop的数组,素数位置为1,其余位置为0
// 从第一个非素数开始,直到大于stop的平方根
// 方法是挨着划去当前数的倍数
int *primes(int stop)
{
    int *numlist = new int[stop+1];
    int i;

    for(i=0;i<=stop;++i)
    {
        if(i<2)
            numlist[i] = 0;
        else
            numlist[i] = i;
    }

    i = 2;
    while(i*i<=stop)
    {
        if(numlist[i]>0)
        {
            for(int j=i+1;j<=stop;++j)
                if(numlist[j]%i==0)
                    numlist[j] = 0;
        }
        i++;
    }

    return numlist;
}

// 判断所有循环是否都是素数
bool checkRotation(int n, int *pmlist)
{
    bool flag = true;
    if(n<10) return flag; // 小于10,直接返回true

    char strnum[7] = ""; // <=1,000,000,7位就够
    sprintf(strnum,"%d",n); // 转为字符串

    int len = strlen(strnum); // 长度
    int curNum,i,j;

    char sNewNum[7] = ""; // 旋转数的字符串形式

    for(i=1;i     {
        // 从位置i开始的旋转数
        for(j=0;j         {
            if(i+j                 sNewNum[j] = strnum[i+j];
            else
                sNewNum[j] = strnum[i+j-len];
        }

        sscanf(sNewNum,"%d",&curNum); // 转为数字
        //cout<
        // 判断是否是素数
        if(pmlist[curNum]==0)
        {
            flag = false;
            break;
        }
    }

    return flag;
}

// euler35
int main(int args, char* argv[])
{
    int stop = 1000000;
    int i,num = 0;
    int *pmlist = primes(stop);

    for(i=0;i<=stop;++i)
        if(pmlist[i]>0 && checkRotation(i,pmlist))
            num++;

    cout<<"小于 "<     return 0;
}

结果
CODE:
//小于 1000000 满足条件的数有: 55 个.
//Process returned 0 (0x0)   execution time : 1.703 s
//Press any key to continue.

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
4楼2011-07-07 03:10:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见