NYOJ219 a problem about date(给出日期计算周几)

给出日期计算周几。

公式留作备用。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstdlib>
using namespace std;
int Hash[7]= {5,6,0,1,2,3,4};
int run[13]= {0,31,29,31,30,31,30,31,31,30,31,30,31};
int pin[13]= {0,31,28,31,30,31,30,31,31,30,31,30,31};
int yy[10000];
//暴力方法，闰年预处理
void init() {
    yy[1600]=1;
    for(int i=1601; i<=9600; i++) {
        if(i%400==0||(i%100&&i%4==0)) {
            yy[i]=yy[i-1]+1;
        } else {
            yy[i]=yy[i-1];
        }
    }
}
int main() {
    int y,m,d;
    init();
    while(~scanf("%d%d%d",&y,&m,&d)) {
        int ans=(y-1600)*365;
        ans+=yy[y-1];

        if(y%400==0||(y%100!=0&&y%4==0)) {
            for(int i=0; i<m; i++) {
                ans+=run[i];
            }
            ans+=d;
        } else {
            for(int i=0; i<m; i++) {
                ans+=pin[i];
            }
            ans+=d;
        }
        printf("%d\n",Hash[ans%7]);
    }
    return 0;
}

//公式<span style="display: none; width: 0px; height: 0px;" id="transmark"></span>
int main(){
    int y, m, d, w;
    while(scanf("%d%d%d", &y, &m, &d) == 3){
        if(m < 3) m += 12, --y;
        w = ( d + 2*m + 3*(m+1)/5 + y + y/4 - y/100 + y/400 + 1) % 7;
        printf("%d\n", w);
    }
    return 0;
}