七个7问题问题的编程解法
七个7问题见http://bbs.emath.ac.cn/thread-73-1-1.html
现在要求使用编程的办法解决,要求求解的时间尽可能的短,
程序尽可能写的比较清晰.
可以使用mathematica\matlab之类的编程语言!
请管理员保留这个帖子,删除另外的那个重复的! 实在不是我想发两遍!而是我的浏览器或者论坛反映太慢,以致于我以为自己没有把帖子发出去 /*此程序使用pari/gp运行*/
/*程序编制时间:2012年4月26日*/
GetDigits(n,d)=
{
/*获得整数n从右向左数的第d位数字,
比如GetDigits(7458,3)=4 */
(n%10^d-n%10^(d-1))/10^(d-1)
}
{
for(a=1,9,
for(b=0,9,
for(c=0,9,
for(d=0,9,
for(e=0,9,
for(f=1,9,
for(g=0,9,
for(h=0,9,
for(i=0,9,
chushu=a*10^5+b*10^4+c*10^3+d*10^2+70+e;\\形成除数
shang=f*10^4+g*10^3+700+h*10+i;\\形成商
beichushu=chushu*shang;\\得到被除数,也就是两者的成绩
if((beichushu<10^9)||(beichushu>10^10),next);\\如果被除数不是10位数,则继续下一个循环
if(GetDigits(beichushu,8)!=7,next);\\如果被除数倒数第八位不是7,那么继续下一个循环
x4=GetDigits(beichushu,4);\\得到被除数的千位
x3=GetDigits(beichushu,3);\\得到被除数的百位
x2=GetDigits(beichushu,2);\\得到被除数的十位
x1=GetDigits(beichushu,1);\\得到被除数的个位
num=(beichushu-beichushu%10000)/10000;\\被除数舍掉后四位
\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
b0=chushu*f;
if(b0<10^5||b0>10^6,next);
\\第01次除法
a1=num-chushu*f;
if(a1<10^5||a1>10^6,next);
if(GetDigits(a1,1)!=7,next);
a12=a1*10+x4;
\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
b1=chushu*g;
if(b1<10^6||b0>10^7,next);
\\第02次除法
a2=a12-b1;
if(a2<10^4||a2>10^5,next);
if(GetDigits(a2,4)!=7,next);
a22=a2*10+x3;
\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
b2=chushu*7;
if(b2<10^5||b2>10^6,next);
\\第03次除法
if(GetDigits(chushu*7,5)!=7,next);
a3=a22-chushu*7;
if(a3<10^5||a2>10^6,next);
a32=a3*10+x2;
\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
b3=chushu*h;
if(b3<10^6||b3>10^7,next);
\\第04次除法
if(GetDigits(chushu*h,3)!=7,next);
a4=a32-chushu*h;
if(a4<10^4||a2>10^5,next);
a42=a4*10+x1;
if((a42-chushu*i)!=0,next);
\\结束上面的所有的操作,并且打印结果
print()
)
)
)
)
)
)
)
)
)
}
这是一个马虎的结果,其实我觉得代码还可以简单地优化一下,这样才更方便阅读! 虽然能得到结果,但是运行的时间实在是太长了,连我自己都受不了了! http://bbs.emath.ac.cn/viewthread.php?tid=73&page=4&fromuid=865#pid42933 3# mathematica
一年多后,自己写的代码自己都不怎么看明白了,这真是让人无语了! mathematica 发表于 2012-4-26 14:05
/*此程序使用pari/gp运行*/
/*程序编制时间:2012年4月26日*/
GetDigits(n,d)=
我现在自己都不明白为什么代码外面要加一层大括号，不加大括号都没办法运行！
页:
[1]