提出 #63890122

---

ソースコード 拡げる

/**
 *    author:  USERNAME
 *    create-mm-dd hh:mm:ss
 **/

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

// clang-format off
/* accelration */

// 高速バイナリ生成
//  #pragma GCC optimize("O3")
//  #pragma GCC optimize("unroll-loops")

// cin cout の結びつけ解除, stdioと同期しない(入出力非同期化)
// cとstdの入出力を混在させるとバグるので注意
struct Fast {Fast() {std::cin.tie(0); ios::sync_with_stdio(false);}} fast;

/* alias */
using ull = unsigned long long;
using ll = long long;
using vi = vector<int>;
using vd = vector<double>;
using vl = vector<long long>;
using vvi = vector<vi>;
using vvl = vector<vl>;
using vvvl = vector<vvl>;
using vvd = vector<vd>;
using vs = vector<string>;
using pii = pair<int, int>;
using in = int;
using qi = queue<int>;
using qvi = queue<vi>;


/* define short */
#define pb push_back
// #define mp make_pair
#define all(obj) (obj).begin(), (obj).end()
#define allr(obj) (obj).rbegin(), (obj).rend()
#define YESNO(bool) if(bool){cout<<"YES"<<endl;}else{cout<<"NO"<<endl;}
#define yesno(bool) if(bool){cout<<"yes"<<endl;}else{cout<<"no"<<endl;}
#define YesNo(bool) if(bool){cout<<"Yes"<<endl;}else{cout<<"No"<<endl;}
// #define INF INFINITY;
const ll INF = 1e16;

/* REP macro */
#define reps(i, a, n) for (ll i = (a); i < (ll)(n); ++i)
#define rep(i, n) reps(i, 0, n)
#define rrep(i, n) reps(i, 1, n + 1)
#define repd(i,n) for(ll i=n-1;i>=0;i--)
#define rrepd(i,n) for(ll i=n;i>=1;i--)

/* debug */
// 標準エラー出力を含む提出はrejectされる場合もあるので注意
// #define debug(x) cerr << "\033[33m(line:" << __LINE__ << ") " << #x << ": " << x << "\033[m" << endl;

/* func */
inline int in_int() {int x; cin >> x; return x;}
inline ll in_ll() {ll x; cin >> x; return x;}
inline string in_str() {string x; cin >> x; return x;}
inline auto  time() {return chrono::system_clock::now();}
inline auto  msec(auto time) {return  chrono::duration_cast<chrono::milliseconds>(time).count(); }
inline int random_number(int min, int max) {
    std::random_device rd;
    std::mt19937 gen(rd());
    std::uniform_int_distribution<> dis(min, max);
    return dis(gen);
}
inline int rand(int max) {
    std::random_device rd;
    std::mt19937 gen(rd());
    std::uniform_int_distribution<> dis(0, max);
    return dis(gen);
}


// search_length: 走査するベクトル長の上限(先頭から何要素目までを検索対象とするか、1始まりで)
template <typename T> inline bool vector_finder(std::vector<T> vec, T element, unsigned int search_length) {
    auto itr = std::find(vec.begin(), vec.end(), element);
    size_t index = std::distance( vec.begin(), itr );
    if (index == vec.size() || index >= search_length) {return false;} else {return true;}
}
template <typename T> inline void print(const vector<T>& v, string s = " ")
    {rep(i, v.size()) cout << v[i] << (i != (ll)v.size() - 1 ? s : "\n");}
template <typename T, typename S> inline void print(const pair<T, S>& p)
    {cout << p.first << " " << p.second << endl;}
template <typename T> inline void print(const T& x) {cout << x << "\n";}
template <typename T, typename S> inline void print(const vector<pair<T, S>>& v)
    {for (auto&& p : v) print(p);}
// 第一引数と第二引数を比較し、第一引数(a)をより大きい/小さい値に上書き
template <typename T> inline bool chmin(T& a, const T& b) {bool compare = a > b; if (a > b) a = b; return compare;}
template <typename T> inline bool chmax(T& a, const T& b) {bool compare = a < b; if (a < b) a = b; return compare;}

/*********       T I M E R       **********/



const int MAX = 50;
const int MOD = 1000000007;

struct Item {
    int origin;   // 出所：社員 j
    int idx;      // どちらのアイテムか（0: a_j用, 1: b_j用）
    int weight;   // 重み = T[j]
    int assigned; // 割り当て先の社員 i
};

ll N, L;
// Timer (1.90秒以内で走らせるための簡易実装)
struct Timer {
    double time_limit;
    clock_t start;
    Timer(double t) : time_limit(t) {
        start = clock();
    }
    bool is_time_up() {
        double elapsed = (double)(clock() - start) / CLOCKS_PER_SEC;
        return elapsed >= time_limit;
    }
    double get_elapsed_ratio() {
        double elapsed = (double)(clock() - start) / CLOCKS_PER_SEC;
        return min(1.0, elapsed / time_limit);
    }
};

// // calcScore：以前ご紹介した周期検出付きシミュレーションによるスコア計算関数
// #include <bitset>
// #include <unordered_map>
// #include <string>
// #include <cmath>

struct State {
    int current;
    bitset<100> parity; // 各社員の担当回数が奇数なら1、偶数なら0
    bool operator==(const State &other) const {
        return current == other.current && parity == other.parity;
    }
};

struct StateHash {
    size_t operator()(const State &s) const {
        return hash<int>()(s.current) ^ (hash<string>()(s.parity.to_string()) << 1);
    }
};

ll calcScore(int L, const vl& T, const vl& a, const vl& b, vi& count) {
    int N = T.size();
    // vector<int> count(N, 0);
    count[0] = 1;
    State state;
    state.current = 0;
    state.parity.reset();
    state.parity.set(0, true);
    unordered_map<State, pair<int, vector<int>>, StateHash> seen;
    seen[state] = {1, count};
    int step = 1;
    while(step < L) {
        int x = state.current;
        int nxt = (state.parity.test(x)) ? a[x] : b[x];
        state.parity.flip(x);
        state.current = nxt;
        count[nxt]++;
        step++;
        if(seen.find(state) != seen.end()){
            auto [prevStep, prevCount] = seen[state];
            int cycleLength = step - prevStep;
            vector<int> delta(N);
            for (int i = 0; i < N; i++){
                delta[i] = count[i] - prevCount[i];
            }
            int remaining = L - step;
            int cycles = remaining / cycleLength;
            for (int i = 0; i < N; i++){
                count[i] += delta[i] * cycles;
            }
            step += cycleLength * cycles;
            break;
        } else {
            seen[state] = {step, count};
        }
    }
    while(step < L){
        int x = state.current;
        int nxt = (state.parity.test(x)) ? a[x] : b[x];
        state.parity.flip(x);
        state.current = nxt;
        count[nxt]++;
        step++;
    }
    
    ll error = 0;
    for (int i = 0; i < N; i++){
        error += abs(count[i] - T[i]);
    }
    return error;
}







int main(){
    int N, L;
    cin >> N >> L;
    // 入力：各社員の目標掃除回数 T[i]
    vvl T_arr(N, vl(2));
    vl t(N);
    rep(i, N) {
        cin >> T_arr[i][0];
        t[i] = T_arr[i][0];
        T_arr[i][1] = i; // 社員番号を保持
    }
    
    // ※「スコアの高い順」にソート（ここでは T_arr[i][0] が大きいほど良いと仮定）
    sort(all(T_arr), [](const vl &a, const vl &b){
        return a[0] > b[0]; // 降順ソート
    });
    
    // 初期解の作成：b はサイクル構造、a は基本的に各社員自身（ただし最後は例外）
    vl a(N), b(N);
    rep(i, N-1) {
        b[T_arr[i][1]] = T_arr[i+1][1]; // 上位から下位へ
    }
    b[T_arr[N-1][1]] = T_arr[0][1]; // サイクルの閉じる部分
    
    rep(i, N-1) {
        a[T_arr[i][1]] = T_arr[i][1]; // 初期は自己ループ
    }
    a[T_arr[N-1][1]] = T_arr[0][1]; // 一部例外
    
    // 初期スコア計算

    vi best_count(N, 0);
    ll current_score = calcScore(L, t,a, b, best_count);
    ll best_score = current_score;
    vl best_a = a;
    vl best_b = b;
    
    

    // 近いものの結合
    int diff = 0; //差が50以下なら結合する
    int connect = 0;
    vl connected(N, 0);
    rep(i,N) connected[i] = i;

    rep(i, N-1){
        if(abs(T_arr[i][0] - T_arr[i+1][0]) < diff){
            a[T_arr[i][1]] = T_arr[i+1][1];
            connect++;
            connected[T_arr[i][1]] = 1;
        }
    }  


    Timer timer(1.80); // 1.90秒以内で探索
    // 焼きなましパラメータ
    double T0 = 1000.0;
    double T_end = 1e-4;
    int cnt = 0;
    
    // a_bの組み合わせを変更する.
    vector<bool> changed(N,false);

    
    while (timer.get_elapsed_ratio() < 1.00) {
        double ratio = timer.get_elapsed_ratio();
        double temp = T0 * pow(T_end / T0, ratio);
        
        // 局所操作：ランダムに T_arr 内のインデックス i (1 <= i < N) を選び、
        // a[T_arr[i][1]] を T_arr 内のより上位 (j < i) の社員 T_arr[j][1] に変更
        
        
        int i = 1 + rand() % (N - 1);  // i = 1,...,N-1
        int pos = T_arr[i][1];  

        // 1から10までをランダムに選ぶ
        if(rand() % 10 < 11){
                        int loop_count = 0;
                    while(connected[i]!=0 && loop_count < 100){
                        i = 1 + rand() % (N - 1);  // i = 1,...,N-1
                        loop_count++;
                    }
                    
                         // 対象社員 
                    
                    int j = rand() % i; 
                    loop_count = 0;
                    while(connected[j]!=0 && loop_count < 100){
                        j = rand() % i; //
                        loop_count++;
                    }
                    int new_val = T_arr[j][1];
                    
                    // 変更前の値を保存
                    int old_val;
                    if (changed[pos]){
                        old_val = b[pos];
                        b[pos] = new_val;
                        
                    } else {
                        
                        old_val = a[pos];
                        a[pos] = new_val;
                    }
                

                    vi tmp_count(N, 0);
                    ll new_score = calcScore(L, t,a, b, tmp_count);
                    int delta = new_score - current_score;
                    
                    // 改善するか、温度に応じた確率で受容. スコアが良いほど
                    if(delta < 0 || exp(-delta / temp) > ((double)rand() / RAND_MAX)) {
                        current_score = new_score;
                        if(new_score < best_score) {
                            best_score = new_score;
                            best_a = a;
                            best_b = b;
                            best_count = tmp_count;
                    // cnt++;
                        }
                    } else {
                        // 不採用の場合、元に戻す
                        // a[pos] = old_val;
                        if (changed[pos]){
                            
                            
                            b[pos] = old_val;
                        } else {
                            a[pos] = old_val;
                        }
                    }
                    
        } else {

            // iのa , bを入れ替える
            int old_a = a[pos];
            int old_b = b[pos];
            a[pos] = old_b;
            b[pos] = old_a;
            changed[pos] = !changed[pos];

             vi tmp_count(N, 0);
            ll new_score = calcScore(L, t,a, b, tmp_count);
            int delta = new_score - current_score;
            
            // 改善するか、温度に応じた確率で受容. スコアが良いほど
            if(delta < 0 || exp(-delta / temp) > ((double)rand() / RAND_MAX)) {
                current_score = new_score;
                if(new_score < best_score) {
                    best_score = new_score;
                    best_a = a;
                    best_b = b;
                    best_count = tmp_count;
            // cnt++;
                }
            } else {
                // 不採用の場合、元に戻す
                // a[pos] = old_val;
                a[pos] = old_a;
                b[pos] = old_b;
            }

          
        }

        
        // cnt++;
        
    }

    while(!timer.is_time_up()){
        double ratio = timer.get_elapsed_ratio();
        double temp = T0 * pow(T_end / T0, ratio);
        //  a , bをランダムに入れ替える
        int i = 1 + rand() % (N - 1);  // i = 1,...,N-1
        int pos = T_arr[i][1];    // 対象社員
        int a_pos = a[pos];
        int b_pos = b[pos];
        a[pos] = b_pos;
        b[pos] = a_pos;

        vi tmp_count(N, 0);
        ll new_score = calcScore(L, t,a, b, tmp_count);
        int delta = new_score - current_score;
        // 改善するか、温度に応じた確率で受容. スコアが良いほど
        if(delta < 0 || exp(-delta / temp) > ((double)rand() / RAND_MAX)) {
            current_score = new_score;
            if(new_score < best_score) {
                best_score = new_score;
                best_a = a;
                best_b = b;
                best_count = tmp_count;
            cnt++;
            }
        } else {
            // 不採用の場合、元に戻す
            a[pos] = a_pos;
            b_pos = b_pos;
        }
        // cnt++;
    }
    
    // 最終スコアと解を出力
    rep(i, N) {
        cout << best_a[i] << " " << best_b[i] << endl;
    }
    // cout << 1e6 - best_score << endl;
    // cout<<cnt<<endl;
    
    return 0;
}

提出情報

コンパイルエラー

Main.cpp:64:19: warning: use of ‘auto’ in parameter declaration only available with ‘-std=c++20’ or ‘-fconcepts’
   64 | inline auto  msec(auto time) {return  chrono::duration_cast<chrono::milliseconds>(time).count(); }
      |                   ^~~~

ジャッジ結果