結果

問題 No.1358 [Zelkova 2nd Tune *] 語るなら枚数を...
ユーザー nikkukunnikkukun
提出日時 2021-01-22 22:25:46
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,267 bytes
コンパイル時間 2,558 ms
コンパイル使用メモリ 167,188 KB
実行使用メモリ 10,624 KB
最終ジャッジ日時 2024-12-28 02:39:54
合計ジャッジ時間 17,423 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 12 TLE * 5
権限があれば一括ダウンロードができます

ソースコード

diff # 

//
#include <bits/stdc++.h>
using namespace std;

#define fi first
#define se second
#define all(x) x.begin(), x.end()
#define lch (o << 1)
#define rch (o << 1 | 1)

typedef double db;
typedef long long ll;
typedef unsigned int ui;
typedef pair<int, int> pint;
typedef tuple<int, int, int> tint;

const int MOD = 1e9 + 7;
const int INF = 0x3f3f3f3f;
const ll INF_LL = 0x3f3f3f3f3f3f3f3f;

void ExGCD(ll a, ll b, ll &x, ll &y, ll &g) {
	if (!b)
		x = 1, y = 0, g = a;
	else
		ExGCD(b, a % b, y, x, g), y -= x * (a / b);
}

ll divCeil(ll upp, ll low) {
	if (upp < 0)
		return upp / low;
	else
		return (upp + low - 1) / low;
}

ll divFloor(ll upp, ll low) {
	if (upp > 0)
		return upp / low;
	else
		return (upp - low + 1) / low;
}

void solve() {
	ll a, b, c, s;
	cin >> a >> b >> c >> s;

	if (s / a < s / c) swap(a, c);
	if (s / b < s / c) swap(b, c);

	ll sum = 0;
	for (ll z = 0; z * c <= s; z++) {
		ll p = s - z * c;
		if (p == 0) sum++;
		if (p <= 0) continue;

		ll x, y, g;
		ExGCD(a, b, x, y, g);
		if (p % g) continue;

		ll L = divCeil(-x * (p / g), b / g);
		ll R = divFloor(y * (p / g), a / g);
		sum += max(0LL, R - L + 1);
	}
	cout << sum % MOD << endl;
}

int main() {
	ios::sync_with_stdio(0);

	int t;
	cin >> t;
	while (t--) solve();

	return 0;
}
0