結果
問題 | No.1203 お菓子ゲーム |
ユーザー |
![]() |
提出日時 | 2020-08-28 23:41:33 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 60 ms / 2,000 ms |
コード長 | 4,813 bytes |
コンパイル時間 | 5,669 ms |
コンパイル使用メモリ | 215,240 KB |
最終ジャッジ日時 | 2025-01-13 19:44:59 |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 51 |
ソースコード
#pragma GCC optimize ("Ofast")#include<bits/stdc++.h>using namespace std;void *wmem;char memarr[96000000];template<class S, class T> inline S min_L(S a,T b){return a<=b?a:b;}template<class T> inline void walloc1d(T **arr, int x, void **mem = &wmem){static int skip[16] = {0, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1};(*mem) = (void*)( ((char*)(*mem)) + skip[((unsigned long long)(*mem)) & 15] );(*arr)=(T*)(*mem);(*mem)=((*arr)+x);}inline int my_getchar_unlocked(){static char buf[1048576];static int s = 1048576;static int e = 1048576;if(s == e && e == 1048576){e = fread_unlocked(buf, 1, 1048576, stdin);s = 0;}if(s == e){return EOF;}return buf[s++];}inline void rd(int &x){int k;int m=0;x=0;for(;;){k = my_getchar_unlocked();if(k=='-'){m=1;break;}if('0'<=k&&k<='9'){x=k-'0';break;}}for(;;){k = my_getchar_unlocked();if(k<'0'||k>'9'){break;}x=x*10+k-'0';}if(m){x=-x;}}inline void rd(long long &x){int k;int m=0;x=0;for(;;){k = my_getchar_unlocked();if(k=='-'){m=1;break;}if('0'<=k&&k<='9'){x=k-'0';break;}}for(;;){k = my_getchar_unlocked();if(k<'0'||k>'9'){break;}x=x*10+k-'0';}if(m){x=-x;}}inline int rd_int(void){int x;rd(x);return x;}struct MY_WRITER{char buf[1048576];int s;int e;MY_WRITER(){s = 0;e = 1048576;}~MY_WRITER(){if(s){fwrite_unlocked(buf, 1, s, stdout);}}};MY_WRITER MY_WRITER_VAR;void my_putchar_unlocked(int a){if(MY_WRITER_VAR.s == MY_WRITER_VAR.e){fwrite_unlocked(MY_WRITER_VAR.buf, 1, MY_WRITER_VAR.s, stdout);MY_WRITER_VAR.s = 0;}MY_WRITER_VAR.buf[MY_WRITER_VAR.s++] = a;}inline void wt_L(char a){my_putchar_unlocked(a);}inline void wt_L(int x){int s=0;int m=0;char f[10];if(x<0){m=1;x=-x;}while(x){f[s++]=x%10;x/=10;}if(!s){f[s++]=0;}if(m){my_putchar_unlocked('-');}while(s--){my_putchar_unlocked(f[s]+'0');}}inline void wt_L(long long x){int s=0;int m=0;char f[20];if(x<0){m=1;x=-x;}while(x){f[s++]=x%10;x/=10;}if(!s){f[s++]=0;}if(m){my_putchar_unlocked('-');}while(s--){my_putchar_unlocked(f[s]+'0');}}template<class T> int Factor_L(T N, T fac[], int fs[]){T i;int sz = 0;if(N%2==0){fac[sz] = 2;fs[sz] = 1;N /= 2;while(N%2==0){N /= 2;fs[sz]++;}sz++;}for(i=3;i*i<=N;i+=2){if(N%i==0){fac[sz] = i;fs[sz] = 1;N /= i;while(N%i==0){N /= i;fs[sz]++;}sz++;}}if(N > 1){fac[sz] = N;fs[sz] = 1;sz++;}return sz;}template<class T> int Divisor_L(T N, T res[], void *mem = wmem){int i;int j;int k;int s;int sz = 0;T *fc;int *fs;int fsz;walloc1d(&fc, 100, &mem);walloc1d(&fs, 100, &mem);fsz =Factor_L(N, fc, fs);res[sz++] = 1;for(i=(0);i<(fsz);i++){s = sz;k = s * fs[i];for(j=(0);j<(k);j++){res[sz++] = res[j] * fc[i];}}sort(res, res+sz);return sz;}int ds;long long d[1000000];int main(){int Lj4PdHRW;wmem = memarr;int KL2GvlyY = rd_int();for(Lj4PdHRW=(0);Lj4PdHRW<(KL2GvlyY);Lj4PdHRW++){int i;long long X;long long Y;long long res = 0;__int128_t a;__int128_t b;__int128_t k;__int128_t up;__int128_t dw;rd(X);rd(Y);if(X==Y-X){wt_L(0);wt_L('\n');continue;}X =min_L(X, Y-X);res += 100000000 / Y;ds =Divisor_L(Y, d);for(i=(0);i<(ds);i++){if(d[i]%2==0){continue;}k = d[i] / 2;up = (2*k*k+2*k) * Y;dw = d[i] * X;if(up % dw){continue;}up /= dw;if(up == d[i]){continue;}if(up < 1 || up > 100000000){continue;}res++;}wt_L(res);wt_L('\n');}return 0;}// cLay varsion 20200813-1 [beta]// --- original code ---// int ds; ll d[1d6];// {// REP(rd_int()){// ll X, Y, res = 0;// __int128_t a, b, k, up, dw;// rd(X, Y);//// if(X==Y-X) wt(0), continue;// X = min(X, Y-X);// res += 1d8 / Y;//// ds = Divisor(Y, d);// rep(i,ds){// if(d[i]%2==0) continue;// k = d[i] / 2;// up = (2*k*k+2*k) * Y;// dw = d[i] * X;// if(up % dw) continue;// up /= dw;// if(up == d[i]) continue;// if(up < 1 || up > 1d8) continue;// res++;// }//// wt(res);// }// }