#542118 (C++17) No.1200 お菓子配り-3

• 問題: No.1200 お菓子配り-3
• ユーザー: chocopuu
• 提出日時: 2020-08-28 22:05:50
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: WA
• 実行時間: -
• コード長: 2,579 bytes
• コンパイル時間: 2,477 ms
• コンパイル使用メモリ: 195,364 KB
• 実行使用メモリ: 5,376 KB
• 最終ジャッジ日時: 2024-04-26 23:02:02
• 合計ジャッジ時間: 27,746 ms
• ジャッジサーバーID （参考情報）: judge5 / judge4

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 5,248 KB
testcase_01	AC	2 ms; 5,248 KB
testcase_02	AC	4 ms; 5,248 KB
testcase_03	AC	4 ms; 5,376 KB
testcase_04	AC	4 ms; 5,376 KB
testcase_05	AC	4 ms; 5,376 KB
testcase_06	AC	4 ms; 5,376 KB
testcase_07	AC	21 ms; 5,376 KB
testcase_08	AC	22 ms; 5,376 KB
testcase_09	AC	21 ms; 5,376 KB
testcase_10	AC	21 ms; 5,376 KB
testcase_11	AC	22 ms; 5,376 KB
testcase_12	AC	185 ms; 5,376 KB
testcase_13	AC	193 ms; 5,376 KB
testcase_14	AC	191 ms; 5,376 KB
testcase_15	AC	193 ms; 5,376 KB
testcase_16	AC	191 ms; 5,376 KB
testcase_17	AC	689 ms; 5,376 KB
testcase_18	AC	990 ms; 5,376 KB
testcase_19	AC	228 ms; 5,376 KB
testcase_20	AC	1,905 ms; 5,376 KB
testcase_21	AC	1,900 ms; 5,376 KB
testcase_22	AC	2,249 ms; 5,376 KB
testcase_23	AC	1,892 ms; 5,376 KB
testcase_24	AC	1,904 ms; 5,376 KB
testcase_25	AC	1,892 ms; 5,376 KB
testcase_26	AC	1,891 ms; 5,376 KB
testcase_27	AC	2 ms; 5,376 KB
testcase_28	AC	3,143 ms; 5,376 KB
testcase_29	AC	2,301 ms; 5,376 KB
testcase_30	AC	2,307 ms; 5,376 KB
testcase_31	WA	-
testcase_32	WA	-

ソースコード

#include "bits/stdc++.h"
using namespace std;
#define int long long
#define REP(i, n) for (int i = 0; i < (int)n; ++i)
#define RREP(i, n) for (int i = (int)n - 1; i >= 0; --i)
#define FOR(i, s, n) for (int i = s; i < (int)n; ++i)
#define RFOR(i, s, n) for (int i = (int)n - 1; i >= s; --i)
#define ALL(a) a.begin(), a.end()
#define IN(a, x, b) (a <= x && x < b)
template<class T>istream&operator >>(istream&is,vector<T>&vec){for(T&x:vec)is>>x;return is;}
template<class T>inline void out(T t){cout << t << "\n";}
template<class T,class... Ts>inline void out(T t,Ts... ts){cout << t << " ";out(ts...);}
template<class T>inline bool CHMIN(T&a,T b){if(a > b){a = b;return true;}return false;}
template<class T>inline bool CHMAX(T&a,T b){if(a < b){a = b;return true;}return false;}
constexpr int INF = 1e18;

#define endl '\n'
#define IOS() ios_base::sync_with_stdio(0);cin.tie(0)

// 負の数にも対応した mod
// 例えば -17 を 5 で割った余りは本当は 3 (-17 ≡ 3 (mod. 5))
// しかし単に -17 % 5 では -2 になってしまう
inline long long mod(long long a, long long m) {
    return (a % m + m) % m;
}

// 拡張 Euclid の互除法
// ap + bq = gcd(a, b) となる (p, q) を求め、d = gcd(a, b) をリターンします
long long extGcd(long long a, long long b, long long &p, long long &q) {  
    if (b == 0) { p = 1; q = 0; return a; }  
    long long d = extGcd(b, a%b, q, p);  
    q -= a/b * p;  
    return d;  
}

// 中国剰余定理
// リターン値を (r, m) とすると解は x ≡ r (mod. m)
// 解なしの場合は (0, -1) をリターン
pair<long long, long long> ChineseRem(long long b1, long long m1, long long b2, long long m2) {
  long long p, q;
  long long d = extGcd(m1, m2, p, q); // p is inv of m1/d (mod. m2/d)
  if ((b2 - b1) % d != 0) return make_pair(0, -1);
  long long m = m1 * (m2/d); // lcm of (m1, m2)
  long long tmp = (b2 - b1) / d * p % (m2/d);
  long long r = mod(b1 + m1 * tmp, m);
  return make_pair(r, m);
}

void solve(){
	//配列は必要な分だけ毎回確保する！
	int X, Y;
	cin >> X >> Y;

	int ans = 0;
	auto f = [&](int a) -> int {
		--a;
		if(a == 0) return 0;
		if(1 - a * a == 0) return 0;
		if ((X - a * Y) % (1 - a * a)) return 0;
		int c = (X - a * Y) / (1 - a * a);
		if((X - c) % a) return 0;
		int b = (X - c) / a;
		if(b == 0) return 0;
		return c > 0 && b > 0;
	};
	for(int i = 1; i * i <= X + Y; ++i) {
		if((X + Y) % i) continue;
		ans += f(i);
		if(i != (X + Y) / i) ans += f((X + Y) / i);
	}
	out(ans);
}

signed main(){
	IOS();
	int Q = 1;
	cin >> Q;
	while(Q--)solve();
}