y=function(x){s=strftime;b=ISOdate a=seq(b(x,1,1),t=b(x,12,31),b='d') length(a[s(a,'%u')==3&trimws(s(a,'%e'))%in%c(2,3,5,7,11,13,17,19,23,29,31)])} # test cases sapply(c(1912,1914,1984,1996,2063,2150,2199,2233), y)