#pragma GCC optimize("Ofast") #pragma GCC target("avx2") #define rd_init() char*rp=({char*mmap();mmap(0l,1l<<25,1,2,0,0ll);}) #define rd() ({int _v=0,_c;while(_c=*rp++-48,_c>=0)_v=_v*10+_c;_v;}) #define wt(v) ({ulong _z=v;do*--wp=_z%10+48;while(_z/=10);}) #define wt1(v) ({char wbuf[64],*wp=wbuf+sizeof wbuf;wt(v);write(1,wp,wbuf+sizeof wbuf-wp);}) #define rep(v,e) for(typeof(e) v=0;v>6]|=1ul<>=1){ if(k&1){ b10=b10*a10%p; b9=b9*a9%p; b7=b7*a7%p; b5=b5*a5%p; } a10=a10*a10%p; a9=a9*a9%p; a7=a7*a7%p; a5=a5*a5%p; } m10[i]=b10; m9[i]=b9; m7[i]=b7; m5[i]=b5; } rrep(ci,n){ unsigned cv=m7[ci]; rrep(bi,ci){ unsigned bv=m9[bi]; rrep(ai,bi){ unsigned av=m10[ai]; z+=d[(p*3+q-(cv+bv+av))%p]; } } d[m5[ci]]+=1; } wt1(z); } _exit(0); }