Alguém no JSPerf abandonou uma implementação incrivelmente rápida para verificar os anos bissextos do calendário ISO (link:Ddd manipulações de bits):

function isLeapYear(year) {
  return !(year & 3 || year & 15 && !(year % 25));
}

Usando o Node.js, verifiquei rapidamente com outras duas implementações de uma linha que conheç

function isLeapClassic(y) { return (y % 4 == 0) && !(y % 100 == 0) || (y % 400 == 0); }
function isLeapXOR(y) { return (y % 4 == 0) ^ (y % 100 == 0) ^ (y % 400 == 0); }
function isLeapBitwise(y) { return !(y & 3 || y & 15 && !(y % 25)); }

//quick'n'dirty test on a small range!
//works with negative integers too
for (var i = 1900; i <= 2100; i++) {
    console.log(
        "year = %d,\t%d%d%d",
        i,
        isLeapClassic(i),
        isLeapXOR(i),
        isLeapBitwise(i)
    );
}

Funciona como esperado, mas meu problema é que não consigo entender como. Eu sei((a % b) == (a & (b-1)) quando b é a potência de dois, então(year % 4) == (year & 3), masyear & 15 && !(year % 25) é bastante difícil de entender. Alguém pode me explicar como isso funciona? Alguma referência sobre esta implementação?