#pragma GCC optimize("Ofast") #pragma GCC target("avx2") char*mmap(); #define rd_skip() while(*rp++>=48) #define rd(v) long v=0;{int _c;while(_c=*rp++-48,_c>=0)v=v*10+_c;} #define wt(v) {long _z=v;do*--wp=_z%10+48;while(_z/=10);} long gcd(long a,long b){ if(a==0||b==0){ return a+b; } long x=__builtin_ctzl(a); long y=__builtin_ctzl(b); a>>=x; b>>=y; long c; while(c=a-b){ c>>=__builtin_ctzl(c); if(c<0){ b=-c; }else{ a=c; } } return a<<(x<=y?x:y); } main(){ char*rp=mmap(0l,1l<<25,1,2,0,0ll); rd(p); rd(q); rd_skip(); long r=0; if(p+q==0){ while(*rp){ long x=*rp==48; rd_skip(); long y=*rp==48; rd_skip(); r+=x&y; } }else{ long g=gcd(p,q); long b=__builtin_ctzl(g); unsigned long t,e; { long x=g>>b; t=~0ul/x; e=x; e*=2-x*e; e*=2-x*e; e*=2-x*e; e*=2-x*e; e*=2-x*e; } p=(p>>b)*e; q=(q>>b)*e; long m=(1<>b)*e; y=(y>>b)*e; if(x<=t&&y<=t){ r+=f|~(x^y)&1; } } } } char wbuf[64],*wp=wbuf+sizeof wbuf; wt(r); write(1,wp,wbuf+sizeof wbuf-wp); _exit(0); }