結果

問題 No.5021 Addition Pyramid
コンテスト
ユーザー kaede2020
提出日時 2025-02-25 21:40:44
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,903 ms / 2,000 ms
コード長 3,590 bytes
コンパイル時間 5,264 ms
コンパイル使用メモリ 281,280 KB
実行使用メモリ 6,820 KB
スコア 45,126,461
最終ジャッジ日時 2025-02-25 21:42:29
合計ジャッジ時間 104,705 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 50
権限があれば一括ダウンロードができます

ソースコード

diff # 

//#pragma GCC optimize("Ofast")
//#pragma GCC optimize "O3,omit-frame-pointer,inline"
#include <iostream> // cout, endl, cin
#include <string> // string, to_string, stoi
#include <vector> // vector
#include <algorithm> // min, max, swap, sort, reverse, lower_bound, upper_bound
#include <utility> // pair, make_pair
#include <tuple> // tuple, make_tuple
#include <cstdint> // int64_t, int*_t
#include <cstdio> // printf
#include <map> // map
#include <queue> // queue, priority_queue
#include <set> // set
#include <stack> // stack
#include <deque> // deque
#include <unordered_map> // unordered_map
#include <unordered_set> // unordered_set
#include <bitset> // bitset
#include <cctype> // isupper, islower, isdigit, toupper, tolower
#include <iomanip>//fixed,setprecision
#include <limits.h>//INT_MAX
#include <math.h>//M_PI
#include <random>
#include <regex> // 正規表現
#include <time.h>
#include <fstream>
#include <array>
#include <bit>
#include <chrono>
#include <span>
#include <cmath>
#include <complex>//複素数
//#include <bits/stdc++.h>
using namespace std;
#include <atcoder/all>
using namespace atcoder;
//using mint = modint998244353;
using ll = long long;
using ull = unsigned long long;
 
#define rep(i, n) for (ll i = 0; i < (ll)(n); i++)
const int dx[4] = {0, 1, 0, -1};
const int dy[4] = {1, 0, -1, 0};
random_device rnd;
mt19937 mt_gen(rnd());
int RandInt(int a, int b) {
    return a + mt_gen() % (b - a + 1);
}
const int MAX = 100000000;

// スコア計算
int score_function(int N, vector<vector<int>> A, vector<int> C) {
    vector<vector<int>> B(N, vector<int>(N, 0));
    for (int i = 0; i < N; i++) B[N - 1][i] = C[i];
    for (int i = N - 2; i >= 0; i--) {
        for (int j = 0; j <= i; j++) B[i][j] = (B[i + 1][j] + B[i + 1][j + 1]) % MAX;
    }
    int max_error = 0;
    for (int i = 0; i < N; i++) {
        for (int j = 0; j <= i; j++) max_error = max(max_error, min(abs(A[i][j] - B[i][j]), MAX - abs(A[i][j] - B[i][j])));
    }
    return MAX / 2 - max_error;
}
int main() {
    // 入力
    int N; cin >> N;
    vector<vector<int>> A(N, vector<int>(N, 0));
    for (int i = 0; i < N; i++) {
        for (int j = 0; j <= i; j++) cin >> A[i][j];
    }
    // SA
    double timeLimit = 1.9; 
    auto t0 = chrono::steady_clock::now();
    double T0 = 1000, T_end = 1e-9;

    random_device rnd;
    mt19937 mt_gen(rnd());
    auto RandDouble = [&]() -> double {
        uniform_real_distribution<double> dist(0.0, 1.0);
        return dist(mt_gen);
    };

    vector<int> C(N);
    rep(i, N) C[i] = RandInt(0, MAX-1);
    
    int current_score = score_function(N, A, C);
    int best_score = current_score;
    vector<int> best_C = C;

    while (true) {
        double elapsed = chrono::duration<double>(chrono::steady_clock::now() - t0).count();
        if (elapsed > timeLimit) break;
        
        double T = T0 + (T_end - T0) * (elapsed / timeLimit);
        
        int r = RandInt(0, N - 1);
        int old_val = C[r];
        int new_val = RandInt(MAX/4, MAX/2);
        C[r] = new_val;
        
        int new_score = score_function(N, A, C);
        int delta = new_score - current_score;
        
        if (delta >= 0 || exp(delta / T) > RandDouble()) {
            current_score = new_score;
            if (new_score > best_score) {
                best_score = new_score;
                best_C = C;
            }
        } else {
            C[r] = old_val; 
        }
    }
    cerr << best_score << endl;
    rep(i, N) {
        if (i >= 1) cout << " ";
        cout << best_C[i];
    }
    cout << endl;
    return 0;
}
0