結果
問題 | No.1358 [Zelkova 2nd Tune *] 語るなら枚数を... |
ユーザー |
![]() |
提出日時 | 2021-02-18 18:33:16 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 3,796 bytes |
コンパイル時間 | 4,812 ms |
コンパイル使用メモリ | 274,000 KB |
実行使用メモリ | 4,480 KB |
最終ジャッジ日時 | 2023-10-13 02:22:16 |
合計ジャッジ時間 | 7,055 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | RE | - |
testcase_01 | AC | 2 ms
4,352 KB |
testcase_02 | AC | 2 ms
4,348 KB |
testcase_03 | AC | 1 ms
4,352 KB |
testcase_04 | AC | 2 ms
4,352 KB |
testcase_05 | AC | 1 ms
4,352 KB |
testcase_06 | AC | 11 ms
4,348 KB |
testcase_07 | AC | 8 ms
4,348 KB |
testcase_08 | AC | 8 ms
4,348 KB |
testcase_09 | AC | 8 ms
4,352 KB |
testcase_10 | AC | 8 ms
4,348 KB |
testcase_11 | RE | - |
testcase_12 | RE | - |
testcase_13 | RE | - |
testcase_14 | RE | - |
testcase_15 | RE | - |
testcase_16 | AC | 828 ms
4,352 KB |
ソースコード
// Enjoy your stay. Code by evima #include <atcoder/all> using namespace atcoder; using mint0 = modint; using mint1 = modint1000000007; #include <bits/stdc++.h> using namespace std; using ld = long double; using ll = long long; using vl = vector<ll>; using LOOPVAR_TYPE = ll; #define all(x) (x).begin(), (x).end() #define sq(x) ((x) * (x)) #define sz(x) ll((x).size()) #define GET_MACRO(_1, _2, _3, NAME, ...) NAME #define rep1(i, n) rep2(i, 0, n) #define rep2(i, a, b) for(LOOPVAR_TYPE i = LOOPVAR_TYPE(a); i < LOOPVAR_TYPE(b); i++) #define rep(...) GET_MACRO(__VA_ARGS__, rep2, rep1)(__VA_ARGS__) #define eb emplace_back #define fir first #define sec second #define mp make_pair #define mt make_tuple const ld EPS = 1e-9; const ld PI = 3.14159265358979323846L; const ll INF = 1070000000LL; const ll MOD = 998244353LL;// 1000000007LL; void fast_io() { cin.tie(NULL); ios_base::sync_with_stdio(false); } ll ll_input() { ll x; cin >> x; return x; } string input() { string s; cin >> s; return s; } vl vl_input(int n) { vector<ll> v(n); rep(i, n) cin >> v[i]; return v; } template<class T> bool chmin(T& a, const T& b) { return (b < a) ? (a = b, true) : false; } template<class T> bool chmax(T& a, const T& b) { return (a < b) ? (a = b, true) : false; } template<class T> map<T, ll> freq(const vector<T>& v){ map<T, ll> ret; for(T x: v) ++ret[x]; return ret; } template<class T> vector<T> reversed(vector<T> v) { reverse(all(v)); return v; } template<class T> vector<T> sorted(vector<T> v) { sort(all(v)); return v; } template<class T> vector<T> sub(const vector<T>& v, int from, int to) { vector<T> ret; copy(&v[from], &v[to], back_inserter(ret)); return ret; } template<class T> string str(const T& x, int prec = 20) { stringstream ss; ss << setprecision(prec) << x; return ss.str(); } template<class T> string str(const vector<T>& v, int prec = 20) { stringstream ss; ss << setprecision(prec); rep(i, sz(v)) ss << v[i] << (i < sz(v) - 1 ? " " : ""); return ss.str(); } template<class T> void print1(T&& x, const string& end) { cout << str(x) << end; } void print() { print1("", "\n"); } template<class T, class... U> void print(T&& head, U&&... tail) { print1(head, " "); print(forward<U>(tail)...); } template<class T> void eprint1(T&& x, const string& end) { cout << str(x) << end; } void eprint() { eprint1("", "\n"); } template<class T, class... U> void eprint(T&& head, U&&... tail) { eprint1(head, " "); eprint(forward<U>(tail)...); } template<class... T> void quit(T&&... x) { print(forward<T>(x)...); exit(0); } void yn(bool cnd, const string& yes = "Yes", const string& no = "No") { if(cnd) print(yes); else print(no); } void zip(vl& v) { vl w = v; sort(all(w)); int n = unique(all(w)) - w.begin(); for(ll& x: v) x = lower_bound(&w[0], &w[n], x) - &w[0]; } vl zipped(vl v) { zip(v); return v; } map<ll, ll> zipmap(vl v) { map<ll, ll> ret; sort(all(v)); v.erase(unique(all(v)), v.end()); rep(i, sz(v)) ret[v[i]] = i; return ret; } void solveOne(){ vl v = sorted(vl_input(3)); ll Y = ll_input(); ll g = gcd(v[0], v[1]); ll a = v[0] / g, b = v[1] / g; mint0::set_mod(b); ll inv_a = mint0(a).inv().val(); mint1 ans = 0; rep(i, Y / v[2] + 1){ ll y = Y - i * v[2]; if(y % g) continue; ll c = y / g; // ax + by = c ll x0 = inv_a * c; ll y0 = (c - a * x0) / b; // (x0 - kb, y0 + ka) ll y_lo = (a - -y0 % a) % a; ll x_lo = x0 % b; ll y_hi = y0 + (x0 - x_lo) * a / b; ans += (y_hi - y_lo) / a + 1; assert((x0 - x_lo) * a <= 1e18); } print(ans.val()); } int main(){ fast_io(); int num_tc = 1; cin >> num_tc; rep(tc, 1, num_tc + 1){ // cout << "Case #" << tc << ": " ;// << endl; solveOne(); } }