結果
問題 | No.1358 [Zelkova 2nd Tune *] 語るなら枚数を... |
ユーザー |
|
提出日時 | 2021-01-22 23:39:26 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 213 ms / 2,000 ms |
コード長 | 1,685 bytes |
コンパイル時間 | 819 ms |
コンパイル使用メモリ | 33,920 KB |
実行使用メモリ | 10,880 KB |
最終ジャッジ日時 | 2024-12-28 08:59:14 |
合計ジャッジ時間 | 2,380 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 17 |
ソースコード
#pragma GCC optimize("Ofast")#pragma GCC target("avx,avx2,fma")#include <cstdio>const long long mod = 1e9+7;long long g0;long long first[1000005];long long gcd(long long a,long long b){ return b?gcd(b,a%b):a; }//pair that ax+by = k, modlong long countPair(long long a,long long b,long long k,int left){if(first[left]==-1) return 0;else{long long s = first[left]*a;long long step = a*b/g0;if(s>k) return 0;else return (k-s)/step+1ll;}}void preWork(long long a,long long b){for(int i = 0; i < b; i++) first[i] = -1;int left = 0;int delta = a%b;int ptr = 1;first[0] = 0;long long maxn = a*b/(gcd(a,b)*a);while(ptr<=maxn){left += delta;if(left>=b) left -= b;if(first[left]==-1) first[left] = ptr;else break;ptr++;}}long long solve(long long n,long long k,long long h,long long y){preWork(k,h);long long step = k*h/gcd(k,h);long long ans = 0;long long left = (y+n)%h;long long delta = n%h;for(int i = 0; i*1ll*n<=y; i++){left -= delta;if(left<0) left += h;long long kk = y-i*1ll*n;if(first[left]!=-1){long long s = first[left]*k;if(s<=kk) ans +=((kk-s)/step+1ll);}}return ans%mod;}int main(){int T;scanf("%d",&T);while(T--){long long n,k,h;long long y;scanf("%lld%lld%lld%lld",&n,&k,&h,&y);if(n>=k && n>=h){printf("%lld\n",solve(n,k,h,y));}else if(k>=n && k>=h){printf("%lld\n",solve(k,n,h,y));}else printf("%lld\n",solve(h,k,n,y));}return 0;}