#line 2 "/home/sakflat/CP/_library/cpp/template/template.cpp" //yukicoder@cpp17 #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using ll = long long; using P = pair; const ll MOD = 998244353; const ll MODx = 1000000007; const int INF = (1<<30)-1; const ll LINF = (1LL<<62LL)-1; const double EPS = (1e-10); P ar4[4] = {{0,1},{0,-1},{1,0},{-1,0}}; P ar8[8] = {{-1,-1},{-1,0},{-1,1},{0,-1},{0,1},{1,-1},{1,0},{1,1}}; template vector make_vector(size_t a, T b) { return vector(a, b); } template auto make_vector(size_t a, Ts... ts) { return vector(a, make_vector(ts...)); } /* 確認ポイント cout << fixed << setprecision(n) << 小数計算//n桁の小数表記になる 計算量は変わらないが楽できるシリーズ min(max)_element(iter,iter)で一番小さい(大きい)値のポインタが帰ってくる count(iter,iter,int)でintがiterからiterの間にいくつあったかを取得できる */ /* comment outed because can cause bugs __attribute__((constructor)) void initial() { cin.tie(0); ios::sync_with_stdio(false); } */ #line 2 "/home/sakflat/CP/_library/cpp/template/basic.cpp" ll calc(ll d) { return (ll)((-1 + sqrt(1 + 4*d)) / 2.0); } ll calcAC(ll d){ long long mn = 0; long long mx = 1000000000; while(mx-mn > 1){ long long ce = (mn+mx)/2; if(ce*ce+ce <= d){ mn = ce; }else{ mx = ce; } } return mn; } int main(){ long long mx = 1000000000000000000LL; int c = 0; for(long long i = 500000000 - 900; 500000000 >= i; i++){ for(long long j = -2000; 2000 > j; j++){ long long x = 4*i*i-1000000000+j; if(calc(x) != calcAC(x)){ cout << x << endl; c++; if(c == 100000)return 0; } } } return 0; }