結果

問題 No.317 辺の追加
ユーザー paruki
提出日時 2016-06-27 00:50:06
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 458 ms / 2,000 ms
コード長 2,446 bytes
コンパイル時間 2,152 ms
コンパイル使用メモリ 184,172 KB
実行使用メモリ 5,888 KB
最終ジャッジ日時 2024-10-11 22:20:58
合計ジャッジ時間 10,197 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 38
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
#define FOR(i,j,k) for(int (i)=(j);(i)<(int)(k);++(i))
#define rep(i,j) FOR(i,0,j)
#define each(x,y) for(auto &(x):(y))
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define debug(x) cout<<#x<<": "<<(x)<<endl
#define smax(x,y) (x)=max((x),(y))
#define smin(x,y) (x)=min((x),(y))
#define MEM(x,y) memset((x),(y),sizeof (x))
#define sz(x) (int)(x).size()
typedef long long ll;
typedef pair<int, int> pii;
typedef vector<int> vi;
typedef vector<ll> vll;

class UnionFind{
public:
    UnionFind(int _n):n(_n), cnt(_n), par(_n), rank(_n), size(_n, 1){
        for(int i=0;i<n;++i) par[i]=i;
    }
    int find(int k){
        return (k==par[k])?k:(par[k]=find(par[k]));
    }
    int operator[](int k){
        return find(k);
    }
    int getSize(int k){
        return size[find(k)];
    }
    void unite(int x, int y){
        x = find(x); y = find(y);
        if(x==y) return;
        --cnt;
        if(rank[x] < rank[y]){
            par[x] = y;
            size[y] += size[x];
        }else{
            par[y] = x;
            size[x] += size[y];
            if(rank[y] == rank[x]) ++rank[y];
        }
    }
    int count(){
        return cnt;
    }
private:
    int n, cnt;
    vector<int> par, rank, size;
};

const int INF = INT_MAX / 100;
int dp[2][100001];

int main(){
    int N, M;
    while(cin >> N >> M){
        UnionFind uft(N + 1);
        rep(i, M){
            int u, v;
            scanf("%d%d", &u, &v);
            uft.unite(u, v);
        }

        vi cnt(N + 1);
        FOR(i, 1, N + 1)cnt[uft[i]]++;
        map<int, int> cc;
        FOR(i, 1, N + 1)if(cnt[i])cc[cnt[i]]++;
        rep(i, N + 1)dp[0][i] = INF;
        dp[0][0] = -1;
        int cur = 0, nex = 1;
        each(p, cc){
            int a, b;
            tie(a, b) = p;
            memcpy(dp[nex], dp[cur], sizeof dp[cur]);
            rep(c, a){
                deque<pii> dq;
                for(int k = 0; c + k*a <= N; ++k){
                    if(sz(dq))smin(dp[nex][c + k*a], dq[0].second + k);
                    int val = dp[cur][c + k*a] - k;
                    while(sz(dq) && dq.back().second > val)dq.pop_back();
                    dq.emplace_back(k, val);
                    if(dq[0].first == k - b)dq.pop_front();
                }
            }
            swap(nex, cur);
        }
        FOR(i, 1, N+1)
            printf("%d\n", dp[cur][i] < INF ? dp[cur][i] : -1);
    }
}
0