洛谷P3959 宝藏(模拟退火乱搞)

题意

题目链接

题面好长啊。。。自己看吧。。

Sol

自己想了一个退火的思路，没想到第一次交85，多退了几次就A了哈哈哈

首先把没用的边去掉，然后剩下的边从小到大排序

这样我们就得到了一个选边的序列，我们要求答案强制按照这个序列选

每次退火的时候选两个点交换。

枚举每个点，判断是否能更新答案，

时间复杂度$O(200 * 1000 * N * M)$

/*
*/
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<ctime>
#include<cstring>
#include<algorithm>
#include<vector>
#define Pair pair<int, int>
#define MP(x, y) make_pair(x, y)
#define fi first
#define se second
using namespace std;
const int MAXN = ;
const double eps = 1e-, Dlt = 0.97, INF = 1e9 + ;
inline int read() {
    char c = getchar(); int x = , f = ;
    while(c < '' || c > '') {if(c == '-') f = -; c = getchar();}
    while(c >= '' && c <= '') x = x *  + c - '', c = getchar();
    return x * f;
}
int N, M;
struct Edge {
    int u, v, w;
    bool operator < (const Edge &rhs) const {
        return w < rhs.w;
    }
}E[MAXN];
int link[MAXN][MAXN], num, fa[MAXN];
void unionn(int x, int y) {
    fa[x] = y;
}
int find(int x) {
    if(fa[x] == x) return fa[x];
    else return fa[x] = find(fa[x]);
}
vector<Pair> v[MAXN];
int dfs(int x, int cnt, int fa) {
    int ans = ;
    for(int i = ; i < v[x].size(); i++) {
        int to = v[x][i].fi, w = v[x][i].se;
        if(to != fa) ans += dfs(to, cnt + , x) + w * cnt;
    }
    return ans;
}
int solve() {
    int cur = INF, tot = , base = ;
    for(int i = ; i <= N; i++) fa[i] = i, v[i].clear();
    for(int i = ; i <= M; i++) {
        int x = E[i].u, y = E[i].v;
        int fx = find(x), fy = find(y);
        if(fx == fy) continue;
        tot++; unionn(fx, fy);
        v[x].push_back(MP(y, E[i].w));
        v[y].push_back(MP(x, E[i].w));
    }
    if(tot != N - ) return INF;
    for(int i = ; i <= N; i++)
        cur = min(cur, dfs(i, , ));
    return cur;
}
int main() {
//    freopen("testdata.in", "r", stdin);
    srand((unsigned)time(NULL));
    memset(link, 0x7f, sizeof(link));
    N = read(); M = read();
    if(N == ) {
        puts(""); return ;
    }
    for(int i = ; i <= M; i++) {
        int x = read(), y = read(), w = read();
        link[x][y] = min(link[x][y], w);
        link[y][x] = min(link[y][x], w);
    }
    for(int i = ; i <= N; i++)
        for(int j = i + ; j <= N; j++)
            if(link[i][j] <= INF)
                E[++num] = (Edge) {i, j, link[i][j]};
    sort(E + , E + num + );
    int ans = solve();
    int times = ;
    while(times--) {
        int now = INF;
        for(double T = ; T > eps; T *= Dlt) {
            int x = rand() % num + , y = rand() % num + ;
            //check(x, y);
            swap(E[x], E[y]);
            int nxt = solve();
            if(nxt < ans) {ans = nxt; continue;}
            if(nxt < now || ((exp(now - nxt / T) < rand() / RAND_MAX))) {now = nxt; continue;}
            swap(E[x], E[y]);
        }
    }
    printf("%d", ans);
    return ;
}
/*
4
0 0
0 5000
2354 10000
8787 0
*/