20# wayne
你的C代码未必比下面的快啊，:D
#include<stdio.h>
int main(){
int i,j,k;
for(i=1;i<10;i++)
for(j=0;j<10;j++)
for(k=0;k<10;k++)
printf("%d\t",100001*i+10010*j+1100*k);
}
我总感觉你的方法这次又涉及到了取余运算，而取余的效率应该不及纯粹的加法、乘法
21# chyanog
如果要进一步追求效率，可以将所有的乘法(*)换成递增(+=)。
22# gxqcn
多谢郭大侠关注。这我可不会了，我对C还算不上了解皮毛，业余的很，只是课外看些书了解下简单语法而已。站长对回文数通项公式求法可有什么建议？
建议谈不上，只是如果为应用的话，应以简洁为美，用你在 21# 所给的表达式足矣。
23# chyanog #include<stdio.h>
int main(){
int i,j,k,a=100001;
for(i=1;i<10;i++){
for(j=0;j<10;j++){
for(k=0;k<10;k++){
printf("%d\t",a);
a+=1100;}
a-=990;}
a-=99;}
}不过，for太多了，俺不喜欢
建议谈不上，只是如果为应用的话，应以简洁为美，用你在 21# 所给的表达式足矣。
gxqcn 发表于 2010-8-27 16:37 http://bbs.emath.ac.cn/images/common/back.gif
这个优化编译器能够自动做到
现在常用的IDE是C-Free，以前用的是Dev-C++，相比之下前者更方便点儿，不过根据我的测试，同一个C程序，前者总是比后者稍慢。似乎他们都是MinGW
发现一个恒等式：
http://zdu6vg.bay.livefilestore.com/y1pP0g3nnyLvNJxAzWoazatSjlrZnBBYri-eA8EoUxSEnCgiSfuEMcU2ZWSW15rzUxYBcQ7RxDbzRCYTyLOo41NYzbq9RYsyfQO/QQ%E6%88%AA%E5%9B%BE%E6%9C%AA%E5%91%BD%E5%90%8D.png?psid=1
前提是 x 属于自然数,算是Floor[(x+9)/10]的展开。
不知道正确否，大家可以验证。
Table -
1/100 (45 + 10 x - 5 Cos (2 - 9 x)] -
5 Sqrt Cos (2 - 9 x)] +
5 Cos (1 - 7 x)] -
5 Sqrt Cos (1 - 7 x)] - 5 Cos[\ x] -
5 Cos (2 + x)] - 5 Sqrt Cos (2 + x)] +
5 Cos (1 + 3 x)] -
5 Sqrt Cos (1 + 3 x)] -
Sqrt)] Sin (1 - 2 x)] -
Sqrt)] Sin (1 - 2 x)] -
Sqrt)] Sin (1 + 3 x)] -
Sqrt)] Sin (1 + 8 x)]), {x, 1000}] //
FullSimplify // Union
我在一楼贴的式子在Mathematica里已经不能再化简了，后来我又自定义三角函数和差化积规则来化简，终于看上去更简洁一些，尽管式子依旧有些长。
http://zdu6vg.bay.livefilestore.com/y1pga6os2MQXJdT_qcjGAW3_rBYkgve4Uokm93E4D6CFpD02ERHmKzwVDvgOpzYnHY0AwluLlgT9DU58Wu2wdsZM_HMH9FBrpOD/QQ%E6%88%AA%E5%9B%BE%E6%9C%AA%E5%91%BD%E5%90%8D.png?psid=1
http://public.bay.livefilestore.com/y1pF3SbnxusFUJfnn2i9DDEAz0W4g471rIsdMwtpzAItqmdg8tYRYhvaBsJaSYeykHG4ejdQSFig99FSYzgtP8HrQ/QQ%E6%88%AA%E5%9B%BE%E6%9C%AA%E5%91%BD%E5%90%8D.png?psid=1