本帖最后由 plp626 于 2012-4-11 21:02 编辑
C测试，40万年以内测试用时秒杀- //&cls&@type %~fs0|tcc -run -
-
- // 把tcc.exe 和 lib\msvcrt.def 放在当前目录下；双击本脚本解释执行下面C代码
- int i2date(int i, int *, int *, int *);
- int date2i(int, int, int);
-
- /*
- * 把0-3-1作为参考日历{0-0-0}的0号索引，后面依次类推;
- * 参考日历月份数为 0,1,2,...,11;
- * 11月分闰月和平月；当所在年份+1为闰年时为闰月；
- * 参考日历日期数为 0,1,...,30;
- * 大月最大偏移30,小月最大偏移29;
- * 11月中，闰月最大偏移28；平月最大偏移27；
- */
- int main(){
-
- int i,j;
- int y,t,m,d;
-
- for (i=0; i<1461*100000; i++){ // 0 ~ 40万年索引号
- i2date(i,&y,&m,&d);
- j=date2i(y,m,d);
- if (i!=j) { // 测试，寻找不对称转换的具体值。。
- printf("/%5d:%4d/%2d/%2d\n",i,y,m,d);
- getchar();
- }
- }
-
- return 0;
- }
-
- int i2date(int i, int *year, int *month, int *day){ // 索引转日期
-
- int t,y,m,d;
-
- y=(i*4+999)/1461; // 1/365.2425的前6位小数最佳有理逼近
- // y=(i*99+145)/36159;
- // y=i*33/12053; //1/365.2425的前9位小数最佳有理逼近
-
- t=i-y*365-y/4+y/100-y/400;
- y+=t>>9; // 获得参考年份；若t为负数,将参考年减1
- t=i-y*365-y/4+y/100-y/400; // 获得参考年0月0日的偏移数
-
- m=(t*5+2)/153; // 获得参考月份
- d=t-(m*153+2)/5; // 获得参考日期
-
- d=d+1; // 获得实际日期
- y+=(m+2)/12; // 获得实际年份
- m=(m+2)%12+1; // 获得实际月份
-
- *year=y;
- *month=m;
- *day=d;
- return 0;
- }
-
- int date2i(int y, int m, int d){// 日期转索引
- int i;
-
- m+=9;
- m%=12;
- y-=m/10;
- //上面三句做平移取模，3月作为当年的0月；2月作为上1年的11月；
-
- i=365*y+y/4-y/100+y/400+(m*153+2)/5+d-1;
- // 11月初的偏移为337；337/11≈153/5;
- // (m*153+2)/5 将 {0,1,2,...,11}映射为{0,31,61,...,337}
-
- return i;
- }
可是，这样的对称测试 有意义吗？ 即使在date2index正确的条件下，也存在index2date出错；
即使date2index正确，date2index它也只是保证正确的日期生成正确的索引；
错误的日期也有可能生成正确的索引；2002-2-29号是错的，但它在意义上式2002-3-1；生成的索引也是2002-3-1日的索引；
如此以来，还需验证index2date生成正确的日期。。。 | 
发表于 2012-4-11 18:36
| 
发表于 2012-4-11 21:59
| 
发表于 2012-4-11 22:15
|