C语言根据日期判断星期几(使用基姆拉尔森计算公式)加解释

C语言根据日期判断星期几（使用基姆拉尔森计算公式）算法如下：基姆拉尔森计算公式W=(d+2*m+3*(m+1)/5+y+y/4-y/100+y/400)mod7在公式中d表示日期中的日数，m表示月份数，y表示年数。注意：在公式中有个与其他公式不同的地方：把一月和二月看成是上一年的十三月和十四月，例：如果是2004-1-10则换算成：2003-13-10来代入公式计算。以公元元年为参考，公元元年1月1日为星期一/PREPRE程序如下：/*利用基姆拉尔森计算日期公式w=(d+2*m+3*(m+1)/5+y+y/4-y/100+y/400)*/#includestdio.hvoidCaculateWeekday(inty,intm,intd){intw;//代表星期几if(m==1||m==2){m+=12;y--;}w=(d+1+2*m+3*(m+1)/5+y+y/4-y/100+y/400)%7;switch(w){case0:printf(星期天！\n);break;case1:printf(星期一\n);break;case2:printf(星期二\n);break;case3:printf(星期三\n);break;case4:printf(星期四\n);break;case5:printf(星期五\n);break;case6:printf(星期六\n);break;}}intmain(){intyear,month,day;charch='1';while(ch!='\033'){printf(\n请输入日期：\n格式为：1900,1,1\n);scanf(%d,%d,%d,&year,&month,&day);CaculateWeekday(year,month,day);flushall();printf(按Esc键退出,其他键继续!);ch=getch();printf(\n);flushall();}}运行效果：请输入日期：格式为：1900,1,12008,4,29星期二输入0退出,其他继续:d2008,1,1星期二输入0退出,其他继续:l2008,8,8星期五输入0退出,其他继续:0请按任意键继续...编者注：用来算现在真实日期的星期是没有问题的。原理是根据已知公元1年1月1日的星期数来推算。如果在你的题目中约定了某天是星期几，你要注意那天的星期是否跟真实的星期相同，如果不同，需要考虑相差几天！如果大家觉得不够过瘾，可以看看以下该公式的推导过程，让大家对历法有个更深刻的认识下面我们完全按自己的思路由简单到复杂一步步进行推导……推导之前，先作两项规定：①用y,m,d,w分别表示年月日星期(w=0-6代表星期日-星期六②我们从公元0年1月1日星期日开始一、只考虑最开始的7天，即d=1---7变换到w=0---6很直观的得到:w=d-1二、扩展到整个1月份模7的概念大家都知道了，也没什么好多说的。不过也可以从我们平常用的日历中看出来，在周历里边每列都是一个按7增长的等差数列，如1、8、15、22的星期都是相同的。所以得到整个1月的公式如下：w=(d-1)%7---------公式⑴三、按年扩展由于按月扩展比较麻烦，所以将年扩展放在前面说①我们不考虑闰年，假设每一年都是365天。由于365是7的52倍多1天，所以每一年的第一天和最后一天星期是相同的。也就是说下一年的第一天与上一年的第一天星期滞后一天。这是个重要的结论，每过一年，公式⑴会有一天的误差，由于我们是从0年开始的，所以只须要简单的加上年就可以修正扩展年引起的误差，得到公式如下：w=(d-1+y)%7②将闰年考虑进去每个闰年会多出一天，会使后面的年份产生一天的误差。如我们要计算2005年1月1日星期几，就要考虑前面的已经过的2004年中有多少个闰年，将这个误差加上就可以正确的计算了。根据闰年的定义(能被4整但不能被100整除或能被400整)，得到计算闰年的个数的算式：y/4-y/100+y/400。由于我们要计算的是当前要计算的年之前的闰年数，所以要将年减1，得到了如下的公式w=[d-1+y+(y-1)/4-(y-1)/100+(y-1)/400]%7-----公式⑵现在，我们得到了按年扩展的公式⑵，用这个公式可以计算任一年的1月份的星期四、扩展到其它月考虑这个问题颇费了一翻脑筋，后来还是按前面的方法大胆假才找到突破口。①现在我们假设每个月都是28天，且不考虑闰年有了这个假设，计算星期就太简单了，因为28正好是7的整数倍，每个月的星期都是一样的，公式⑵对任一个月都适用：)②但假设终究是假设，首先1月就不是28天，这将会造成2月份的计算误差。1月份比28天要多出3天，就是说公式⑵的基础上，2月份的星期应该推后3天。而对3月份来说，推后也是3天(2月正好28天，对3月的计算没有影响)。依此类推，每个月的计算要将前面几个月的累计误差加上。要注意的是误差只影响后面月的计算，因为12月已是最后一个月，所以不用考虑12月的误差天数，同理，1月份的误差天数是0，因为前面没有月份影响它。由此，想到建立一个误差表来修正每个月的计算。==================================================月误差累计模7130020333333426653816211473136831629219510321011224312—265(闰年时2月会有一天的误差，但我们现在不考虑)==================================================我们将最后的误差表用一个数组存放在公式⑵的基础上可以得到扩展到其它月的公式e[]={0,3,3,6,1,4,6,2,5,0,3,5}w=[d-1+y+e[m-1]+(y-1)/4-(y-1)/100+(y-1)/400]%7--公式⑶③上面的误差表我们没有考虑闰年，如果是闰年，2月会一天的误差，会对后面的3-12月的计算产生影响，对此，我们暂时在编程时来修正这种情况，增加的限定条件是如果当年是闰年，且计算的月在2月以后，需要加上一天的误差。大概代码是这样的：w=(d-1+y+e[m-1]+(y-1)/4-(y-1)/100+(y-1)/400);if(m2&&(y%4==0&&y%100!=0||y%400==0)&&y!=0)++w;w%=7;现在，已经可以正确的计算任一天的星期了。注意：0年不是闰年，虽然现在大都不用这个条件，但我们因从公元0年开始计算，所以这个条件是不能少的。④改进公式⑶中，计算闰年数的子项(y-1)/4-(y-1)/100+(y-1)/400没有包含当年，如果将当年包含进去，则实现了如果当年是闰年，w自动加1。由此带来的影响是如果当年是闰年，1,2月份的计算会多一天误差，我们同样在编程时修正。则代码如下w=(d-1+y+e[m-1]+y/4-y/100+y/400);----公式⑷if(m3&&(y%4==0&&y%100!=0||y%400==0)&&y!=0)--w;w%=7;与前一段代码相比，我们简化了w的计算部分。实际上还可以进一步将常数-1合并到误差表中，但我们暂时先不这样做。至此，我们得到了一个阶段性的算法，可以计算任一天的星期了。publicclassWeek{publicstaticvoidmain(String[]args){inty=2005;intm=4;intd=25;inte[]=newint[]{0,3,3,6,1,4,6,2,5,0,3,5};intw=(d-1+e[m-1]+y+(y2)-y/100+y/400);if(m3&&((y&3)==0&&y%100!=0||y%400==0)&&y!=0){--w;}w%=7;System.out.println(w);}}五、简化现在我们推导出了自己的计算星期的算法了，但还不能称之为公式。所谓公式，应该给定年月日后可以手工算出星期几的，但我们现在的算法需要记住一个误差表才能进行计算，所以只能称为一种算法，还不是公式。下面，我们试图消掉这个误差表……=============================消除闰年判断的条件表达式=============================由于闰年在2月份产生的误差，影响的是后面的月份计算。如果2月是排在一年的最后的话，它就不能对其它月份的计算产生影响了。可能已经有人联想到了文章开头的公式中为什么1,2月转换为上年的13,14月计算了吧：)就是这个思想了，我们也将1,2月当作上一年的13,14月来看待。由此会产生两个问题需要解决：1一年的第一天是3月1日了，我们要对w的计算公式重新推导2误差表也发生了变化，需要得新计算①推导w计算式1用前面的算法算出0年3月1日是星期3前7天,d=1---7===w=3----2得到w=(d+2)%7此式同样适用于整个三月份2扩展到每一年的三月份[d+2+y+(y-1)/4-(y-1)/100+(y-1)/400]%7②误差表==================================================月误差累计模7330042335355628173103831369216210318411221012323213326514-291==================================================③得到扩展到其它月的公式e[]={0,3,5,1,3,6,2,4,0,2,5,1}w=[d+2+e[m-3]+y+(y-1)/4-(y-1)/100+(y-1)/400]%7(3=m=14)我们还是将y-1的式子进行简化w=[d+2+e[m-3]+y+y/4-y/100+y/400]%7(3=m=14)这个式子如果当年是闰年，会告成多1的误差但我们将1,2月变换到上一年的13,14月，年份要减1，所以这个误差会自动消除，所以得到下面的算法：inte[]=newint[]{0,3,5,1,3,6,2,4,0,2,5,1};if(m3){m+=12;--y;}intw=(d+2+e[m-3]+y+(y/4)-y/100+y/400)%7;-----公式⑸我们可以看到公式⑸与公式⑷几乎是一样的，仅仅是误差天和一个常数的差别常数的区别是由起始日期的星期不同引起的，0年1月1日星期日，0年3日1日星期三，有三天的差别，所以常数也从-1变成了2。现在，我们成功的消除了繁琐的闰年条件判断。=============================消除误差表=============================假如存在一种m到e的函数映射关系，使得e[m-3]=f(m)则我们就可以用f(m)取代公式⑸中的子项e[m-3]，也就消除了误差表。由于误差表只有12个项，且每一项都可以加减7n进行调整，这个函数关系是可以拼凑出来的。但是这个过程可能是极其枯燥无味的，我现在不想自己去推导它，我要利用前人的成果。所谓前人栽树，后人乘凉嘛：)文章开头开出的公式中的2*m+3*(m+1)/5这个子项引起了我的兴趣经过多次试试验，我运行下面的代码for(m=1;m=14;++m)System.out.print((-1+2*m+3*(m+1)/5)%7+);System.out.println();天哪，输出结果与我的误差表不谋而合，成功了，哈哈24035136240251Pressanykeytocontinue...上面就是输出结果，看它后面的12项，与我的误差表完全吻合!!!现在就简单的，将f(m)=-1+2*m+3*(m+1)/5代入公式⑸，得到w=(d+1+2*m+3*(m+1)/5+y+(y/4)-y/100+y/400)%7----公式6约束条件:m=1,m=2时m=m+12,y=y-1;现在，我们得到了通用的计算星期的公式，并且“完全”是按自己的思想推导出来的(那个函数映射关系不算)，只要理解了这个推导