4 ms·

In C-lang the algorithm is more mysterious: dow(m,d,y) { y-=m<3; return(y+y/4-y/100+y/400+"-bed=pen+mad."[m]+d)%7; }