結果

問題 No.3107 Listening
ユーザー harurunharurun
提出日時 2024-03-07 05:18:00
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 5,160 bytes
コンパイル時間 7,687 ms
コンパイル使用メモリ 428,176 KB
実行使用メモリ 68,748 KB
最終ジャッジ日時 2024-09-30 22:05:17
合計ジャッジ時間 41,569 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
13,764 KB
testcase_01 AC 3 ms
6,812 KB
testcase_02 AC 3 ms
6,816 KB
testcase_03 AC 938 ms
38,852 KB
testcase_04 AC 149 ms
10,248 KB
testcase_05 AC 439 ms
14,848 KB
testcase_06 AC 161 ms
12,164 KB
testcase_07 AC 1,000 ms
27,496 KB
testcase_08 AC 636 ms
19,316 KB
testcase_09 AC 491 ms
24,572 KB
testcase_10 AC 625 ms
39,876 KB
testcase_11 AC 966 ms
48,544 KB
testcase_12 AC 630 ms
19,300 KB
testcase_13 AC 964 ms
38,608 KB
testcase_14 AC 836 ms
23,548 KB
testcase_15 AC 246 ms
19,084 KB
testcase_16 AC 308 ms
11,800 KB
testcase_17 AC 672 ms
14,832 KB
testcase_18 AC 93 ms
7,660 KB
testcase_19 AC 693 ms
30,340 KB
testcase_20 AC 584 ms
35,244 KB
testcase_21 AC 1,620 ms
57,012 KB
testcase_22 AC 883 ms
44,028 KB
testcase_23 AC 1,589 ms
54,460 KB
testcase_24 AC 870 ms
49,140 KB
testcase_25 AC 1,209 ms
32,864 KB
testcase_26 AC 495 ms
14,708 KB
testcase_27 AC 443 ms
27,080 KB
testcase_28 AC 1,002 ms
47,244 KB
testcase_29 AC 1,387 ms
49,528 KB
testcase_30 AC 821 ms
23,720 KB
testcase_31 AC 880 ms
27,512 KB
testcase_32 AC 1,127 ms
42,116 KB
testcase_33 TLE -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
testcase_41 -- -
testcase_42 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <algorithm>
#include <deque>
#include <queue>
#include <string>
#include <iomanip>
#include <set>
#include <unordered_set>
#include <map>
#include <unordered_map>
#include <utility>
#include <stack>
#include <random>
#include <complex>
#include <functional>
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>
#include <assert.h>
#if __has_include(<atcoder/all>)
#include <atcoder/all>
using namespace atcoder;
#endif
#if __has_include(<boost/multiprecision/cpp_int.hpp>)
#include <boost/multiprecision/cpp_int.hpp>
using cpp_int=boost::multiprecision::cpp_int;
#endif
#if __has_include(<boost/multiprecision/cpp_dec_float.hpp>)
#include <boost/multiprecision/cpp_dec_float.hpp>
template<unsigned size>using cpp_float=boost::multiprecision::number<boost::multiprecision::cpp_dec_float<size>>;
template<unsigned size>using cpp_double=boost::multiprecision::number<boost::multiprecision::cpp_dec_float<size,long long>>;
#endif
using namespace std;
using ll=long long;
#define read(x) cin>>(x);
#define readll(x) ll (x);cin>>(x);
#define readS(x) string (x);cin>>(x);
#define readvll(x,N) vector<ll> (x)((N));for(int i=0;i<(N);i++){cin>>(x)[i];}
#define rep(i,N) for(ll (i)=0;(i)<(N);(i)++)
#define rep2d(i,j,H,W) for(ll (i)=0;(i)<(H);(i)++)for(ll (j)=0;(j)<(W);j++)
#define is_in(x,y) (0<=(x) && (x)<H && 0<=(y) && (y)<W)
#define yn {cout<<"Yes"<<endl;}else{cout<<"No"<<endl;}
#define double_out(x) fixed << setprecision(x)
template<class T> inline void erase_duplicate(vector<T>& A){sort(A.begin(),A.end());A.erase(unique(A.begin(),A.end()),A.end());}
inline ll powll(ll x,ll n){ll r=1;while(n>0){if(n&1){r*=x;};x*=x;n>>=1;};return r;}

// https://ei1333.github.io/library/graph/flow/gabow-edmonds.hpp.html
// https://qiita.com/Kutimoti_T/items/5b579773e0a24d650bdf

/**
 * @brief Gabow Edmonds(一般グラフの最大マッチング)
 * @docs docs/gabow-edmonds.md
 */
struct GabowEdmonds {

  struct edge {
    int to, idx;
  };

  vector< vector< edge > > g;
  vector< pair< int, int > > edges;
  vector< int > mate, label, first;
  queue< int > que;

  GabowEdmonds(int n) : g(n + 1), mate(n + 1), label(n + 1, -1), first(n + 1) {}

  void add_edge(int u, int v) {
    ++u, ++v;
    g[u].push_back(edge{v, (int) (edges.size() + g.size())});
    g[v].push_back(edge{u, (int) (edges.size() + g.size())});
    edges.emplace_back(u, v);
  }

  int find(int x) {
    if(label[first[x]] < 0) return first[x];
    first[x] = find(first[x]);
    return first[x];
  }

  void rematch(int v, int w) {
    int t = mate[v];
    mate[v] = w;
    if(mate[t] != v) return;
    if(label[v] < (int)g.size()) {
      mate[t] = label[v];
      rematch(label[v], t);
    } else {
      int x = edges[label[v] - g.size()].first;
      int y = edges[label[v] - g.size()].second;
      rematch(x, y);
      rematch(y, x);
    }
  }

  void assign_label(int x, int y, int num) {
    int r = find(x);
    int s = find(y);
    int join = 0;
    if(r == s) return;
    label[r] = -num;
    label[s] = -num;
    while(true) {
      if(s != 0) swap(r, s);
      r = find(label[mate[r]]);
      if(label[r] == -num) {
        join = r;
        break;
      }
      label[r] = -num;
    }
    int v = first[x];
    while(v != join) {
      que.push(v);
      label[v] = num;
      first[v] = join;
      v = first[label[mate[v]]];
    }
    v = first[y];
    while(v != join) {
      que.push(v);
      label[v] = num;
      first[v] = join;
      v = first[label[mate[v]]];
    }
  }

  bool augment_check(int u) {
    que = queue< int >();
    first[u] = 0;
    label[u] = 0;
    que.push(u);
    while(!que.empty()) {
      int x = que.front();
      que.pop();
      for(auto e : g[x]) {
        int y = e.to;
        if(mate[y] == 0 && y != u) {
          mate[y] = x;
          rematch(x, y);
          return true;
        } else if(label[y] >= 0) {
          assign_label(x, y, e.idx);
        } else if(label[mate[y]] < 0) {
          label[mate[y]] = x;
          first[mate[y]] = y;
          que.push(mate[y]);
        }
      }
    }
    return false;
  }

  vector< pair< int, int > > max_matching() {
    for(int i = 1; i < (int)g.size(); i++) {
      if(mate[i] != 0) continue;
      if(augment_check(i)) label.assign(g.size(), -1);
    }
    vector< pair< int, int > > ret;
    for(int i = 1; i < (int)g.size(); i++) {
      if(i < mate[i]) ret.emplace_back(i - 1, mate[i] - 1);
    }
    return ret;
  }
};

struct pair_hash{
  inline size_t operator()(const pair<ll,ll>& p)const{
    return 10'000'000LL*p.first+p.second;
  }
};

int main(){
  ll N,M;
  cin>>N>>M;
  GabowEdmonds GE(N);
  unordered_map<ll,ll> index;
  vector<pair<ll,ll>> Edges(M);
  for(ll i=0;i<M;i++){
    cin>>Edges[i].first>>Edges[i].second;
    Edges[i].first--;
    Edges[i].second--;
    GE.add_edge(Edges[i].first,Edges[i].second);
  }
  auto res=GE.max_matching();
  unordered_set<pair<ll,ll>,pair_hash> ans(res.begin(),res.end());
  cout<<ans.size()<<endl;
  for(ll i=0;i<M;i++){
    if(ans.find(Edges[i])!=ans.end() || ans.find(make_pair(Edges[i].second,Edges[i].first))!=ans.end()){
      cout<<i+1<<endl;
    }
  }
}
0