关于求质数的dc程序翻译

如果我是疯

上次有位哥们发了个dc程序，是求质数的，分析了一下，现转换成C程序放在下面：
原dc程序：
echo '2p3p[dl!d2+s!%0=@l!l^!<#]s#[s/0ds^]s@[p]s&[ddvs^3s!l#x0<&2+l.x]ds.x'|dc

该程序使用的算法我不讨论了，翻译的c程序如下：
#include <stdio.h>
#include <math.h>

int stack[100];
int sflag=0;

int a;        /*^*/
int b;        /*!*/

void check_doing(void); /*#*/
void check_do2(void);   /*@*/
void check_prime(void); /*.*/
void xadd(void);
void xmod(void);
int  get(void);
int pop(void);
void push(int);
void dupli(void);
void xprint(void);    /*&*/

int get(void)
{
        return stack[sflag];
}

int pop(void)
{
        if(sflag<0) return -1000;

        sflag--;       
        return stack[sflag+1];
}

void push(int i)
{
        sflag++;
        stack[sflag]=i;
}

void dupli(void)
{
        push(stack[sflag]);
}

void xprint(void)
{
        printf("%d\n",stack[sflag]);       
}

void check_prime(void)
{

/*now begin the main */
        dupli();
        dupli();
        a=sqrt(pop());

        b=3;
       
        check_doing();
/*end the main*/
       
        if(pop()>0) xprint();

        push(2+pop());

        if(get()>50) return;
       
        check_prime();       
}

void check_doing(void)
{
        int c;

        dupli();
        push(b);
        dupli();
        push(2);
        xadd();
       
        b=pop();
        xmod();

        if(pop()==0) check_do2();
       
        if(!(a<b)) check_doing();
        else return;
}

void check_do2(void)
{
        pop();
        push(0);
        a=0;
}

void xadd(void)
{
        int k;

        k=stack[sflag]+stack[sflag-1];
       
        pop();
        pop();
        push(k);
}

void xmod(void)
{
        int k;
        k=stack[sflag-1]%stack[sflag];
        pop();
        pop();
        push(k);
}

void main(void)
{
        push(2);
        xprint();
        push(3);
        xprint();
        check_prime();
}

如果存为x1.c
请用gcc -g -o x x1.c -lm编译，然后用gdb跟踪就可以了。
原dc程序一直循环，我改成取出50之前的质数，如果使用dc，则改动如下：

echo '2p3p[dl!d2+s!%0=@l!l^!<#]s#[s/0ds^]s@[p]s&[ddvs^3s!l#x0<&2+d6000>.]ds.x'|dc
这样只执行到6000以内；

r2007

好!翻译思路清晰、严谨。推荐！
建议论坛的字体该改一下了，最近总惦记着这段dc代码，直到读了楼主的帖子，才反推出来原dc代码的意思。我一直把其中的l(load)看成管道符，刚解读到.宏命令[ddvs^3s!l#x0<&2+l.x]就一头雾水，以为要做x的y次方再模除z，以为是什么高深的求解质数的方法^_^

hoxide

应该是用了欧拉定理吧,如果p是质数,则对任意2<=r<p有(r^p)%p = p

另dc是什么东西阿?

zhiwood

segment fault on calculating primes within 100000, the last one is 65257.

li2002

谁有dc的中文文档？？
想不到这个东东竟然这么强！

寂寞烈火

dc是一个任意精度的计算器.dc从标准输入读入命令,如果向dc提供参数,那么这个参数就必需是文件名,dc从文件读入并执行,所有的输出将被送到标准输出,所有的错误被送到标准错误设备,输入的数据将存放在堆栈中,每输入一个数据,dc将它放入栈中,如果输入的是算术运算符,那么dc就从栈中读取数据,并计算,将结果返回堆栈
PS:抄自:<<linux命令打全>>

li2002

有dc应用的例子吗

寂寞烈火

[quote]原帖由 "li2002"]有dc应用的例子吗[/quote 发表：

/home/lee#dc
1
2
3
f    //显示栈中所有数据
3
2
1
p    //显示栈顶数据
3
*    //栈顶头2个数据相乘
p   
6
f  //显示当前栈中数据
6
1
- //栈顶头2个数据相减
p
-5
f //显示栈中当前数据
-5
ctr+D //退出
PS:对dc的了解止于此

ylbaby2003

火兄，llinux命令打全，你知道哪里有下吗？

寂寞烈火

[quote]原帖由 "ylbaby2003"]火兄，llinux命令打全，你知道哪里有下吗？[/quote 发表：

这个大全,是我手头的一本纸书,网上应该有的,google一下吧!