#include #include #include #include #include using namespace std; namespace mp = boost::multiprecision; using Real1024 = mp::number>; Real1024 eps = mp::pow((Real1024)10,(Real1024)-12); #define rep(i, n) for (int i = 0; i < (int)(n); i++) #define length size() //#define int long long #define ll long long constexpr ll inf = 1001001001001001001ll; constexpr ll mod = /* 1000000007*/ 998244353; constexpr double pi = 3.141592653589793; //小数出力 //cout << std::setprecision(12); //imos void to_ruiseki(vector &vec){ for(int i=1;i<(int)vec.size();i++){ vec[i] += vec[i-1]; } } void to_2druiseki(vector> &vec){ rep(i,vec.size()){ rep(j,vec[0].size()){ int s = 0; if(i!=0) s += vec[i-1][j]; if(j!=0) s += vec[i][j-1]; if(i!=0 && j!=0) s -= vec[i-1][j-1]; s += vec[i][j]; vec[i][j] = s; } } } //modpow(a,n,mod) int modpow(int a, int n, int mod) { long long res = 1; while (n > 0) { if (n & 1) res = res * a % mod; a = a * a % mod; n >>= 1; } return res; } //素因数分解 vector> p_fact(int N) { vector> res; for (int a = 2; a * a <= N; ++a) { if (N % a != 0) continue; int ex = 0; // 指数 // 割れる限り割り続ける while (N % a == 0) { ++ex; N /= a; } // その結果を push res.push_back({a, ex}); } // 最後に残った数について if (N != 1) res.push_back({N, 1}); return res; } //ctoi int ctoi(const char c){ switch(c){ case '0': return 0; case '1': return 1; case '2': return 2; case '3': return 3; case '4': return 4; case '5': return 5; case '6': return 6; case '7': return 7; case '8': return 8; case '9': return 9; default : throw runtime_error("hoge"); } } //max template< typename T1, typename T2 > inline bool chmax(T1 &a, T2 b) { return a < b && (a = b, true); } //min template< typename T1, typename T2 > inline bool chmin(T1 &a, T2 b) { return a > b && (a = b, true); } //print void print() { cout << '\n'; } template void print(const T &t) { cout << t << '\n'; } template void print(const Head &head, const Tail &... tail) { cout << head << ' '; print(tail...); } //join template string join(vector &vec ,const string &sp){ int si = vec.length; if(si==0){ return ""; }else{ stringstream ss; rep(i,si-1){ ss << vec[i] << sp; } ss << vec[si - 1]; return ss.str(); } } //木の直径 int tree_diameter(vector> &vec){ //0-indexed queue> que; vector is_searched(vec.size(),false); que.push(make_pair(0,0)); int max_d = 0; int max_e = 0; while(!que.empty()){ int e = que.front().first; int distance = que.front().second; is_searched[e] = true; if(max_d> &vec,int st){ queue> que; vector is_searched(vec.size(),false); is_searched[st] = true; que.push(make_pair(st,0)); int max_d = 0; int max_e = 0; while(!que.empty()){ int e = que.front().first; int distance = que.front().second; if(max_d par, rank, siz; // 構造体の初期化 UnionFind(int n) : par(n,-1), rank(n,0), siz(n,1) { } // 根を求める int root(int x) { if (par[x]==-1) return x; // x が根の場合は x を返す else return par[x] = root(par[x]); // 経路圧縮 } // x と y が同じグループに属するか (= 根が一致するか) bool issame(int x, int y) { return root(x)==root(y); } // x を含むグループと y を含むグループを併合する bool unite(int x, int y) { int rx = root(x), ry = root(y); // x 側と y 側の根を取得する if (rx==ry) return false; // すでに同じグループのときは何もしない // union by rank if (rank[rx] > prime_factorize(long long N) { vector > res; for (long long a = 2; a * a <= N; ++a) { if (N % a != 0) continue; long long ex = 0; // 指数 // 割れる限り割り続ける while (N % a == 0) { ++ex; N /= a; } // その結果を push res.push_back({a, ex}); } // 最後に残った数について if (N != 1) res.push_back({N, 1}); return res; } class BIT { public: BIT() = default; // 長さ size の数列で初期化 explicit BIT(size_t size) : m_bit(size + 1) {} // 数列で初期化 explicit BIT(const std::vector& v) : BIT(v.size()) { for (int i = 0; i < (int)v.size(); ++i) { add((i + 1), v[i]); } } // 閉区間 [1, r] の合計を返す (1-based indexing) long long sum(int r) const { long long ret = 0; for (; 0 < r; r -= (r & -r)) { ret += m_bit[r]; } return ret; } // 閉区間 [l, r] の合計を返す (1-based indexing) long long sum(int l, int r) const { return (sum(r) - sum(l - 1)); } // 数列の i 番目の要素を加算 (1-based indexing) void add(int i, long long value) { for (; i < (int)m_bit.size(); i += (i & -i)) { m_bit[i] += value; } } private: std::vector m_bit; }; // Binary Indexed Tree (1.2 区間加算対応) // 1-based indexing class BIT_RAQ { public: BIT_RAQ() = default; explicit BIT_RAQ(size_t size) : m_bit0(size) , m_bit1(size) {} explicit BIT_RAQ(const std::vector& v) : m_bit0(v) , m_bit1(v.size()) {} // 閉区間 [1, r] の合計を返す (1-based indexing) long long sum(int r) const { return (m_bit0.sum(r) + m_bit1.sum(r) * r); } // 閉区間 [l, r] の合計を返す (1-based indexing) long long sum(int l, int r) const { return (sum(r) - sum(l - 1)); } // 数列の i 番目の要素を加算 (1-based indexing) void add(int i, long long value) { m_bit0.add(i, value); } // 閉区間 [l, r] の要素を加算 (1-based indexing) void add(int l, int r, long long value) { m_bit0.add(l, (-value * (l - 1))); m_bit0.add((r + 1), (value * r)); m_bit1.add(l, value); m_bit1.add((r + 1), -value); } private: BIT m_bit0; BIT m_bit1; }; int bs(vector> &vec,int n){ int left = -1; //「index = 0」が条件を満たすこともあるので、初期値は -1 int right = (int)vec.size(); // 「index = a.size()-1」が条件を満たさないこともあるので、初期値は a.size() /* どんな二分探索でもここの書き方を変えずにできる! */ while (right - left > 1) { int mid = left + (right - left) / 2; if (n <= vec[mid].second) right = mid; else left = mid; } /* left は条件を満たさない最大の値、right は条件を満たす最小の値になっている */ return (int)(right < (int)vec.size() && vec[right].first <= n); } signed main(void){ int n; cin >> n; rep(p,n){ Real1024 a,b,c; cin >> a >> b; c = mp::sqrt(a)+mp::sqrt(b); if(mp::abs(c-mp::ceil(c))