結果
問題 | No.679 不思議マーケット |
ユーザー |
![]() |
提出日時 | 2017-09-22 08:53:17 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 6 ms / 2,000 ms |
コード長 | 11,262 bytes |
コンパイル時間 | 2,174 ms |
コンパイル使用メモリ | 195,896 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-08 19:58:36 |
合計ジャッジ時間 | 2,927 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
外部呼び出し有り |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 18 |
コンパイルメッセージ
main.cpp: In function ‘void vizGraph(Graph&, int, std::string)’: main.cpp:170:11: warning: ignoring return value of ‘int system(const char*)’ declared with attribute ‘warn_unused_result’ [-Wunused-result] 170 | system(((string)"dot -T png out.dot >" + filename).c_str()); | ~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
ソースコード
#include <bits/stdc++.h>#include <sys/time.h>using namespace std;#define rep(i,n) for(long long i = 0; i < (long long)(n); i++)#define repi(i,a,b) for(long long i = (long long)(a); i < (long long)(b); i++)#define pb push_back#define all(x) (x).begin(), (x).end()#define fi first#define se second#define mt make_tuple#define mp make_pair#define ZERO(a) memset(a,0,sizeof(a))template<class T1, class T2> bool chmin(T1 &a, T2 b) { return b < a && (a = b, true); }template<class T1, class T2> bool chmax(T1 &a, T2 b) { return a < b && (a = b, true); }#define exists find_if#define forall all_ofusing ll = long long; using vll = vector<ll>; using vvll = vector<vll>; using P = pair<ll, ll>;using ld = long double; using vld = vector<ld>;using vi = vector<int>; using vvi = vector<vi>; vll conv(vi& v) { vll r(v.size()); rep(i, v.size()) r[i] = v[i]; return r; }using Pos = complex<double>;template <typename T, typename U> ostream &operator<<(ostream &o, const pair<T, U> &v) { o << "(" << v.first << ", " << v.second << ")"; return o; }template<size_t...> struct seq{}; template<size_t N, size_t... Is> struct gen_seq : gen_seq<N-1, N-1, Is...>{}; template<size_t... Is> struct gen_seq<0, Is...> : seq<Is...>{};template<class Ch, class Tr, class Tuple, size_t... Is>void print_tuple(basic_ostream<Ch,Tr>& os, Tuple const& t, seq<Is...>){ using s = int[]; (void)s{0, (void(os << (Is == 0? "" : ", ") << get<Is>(t)),0)...}; }template<class Ch, class Tr, class... Args>auto operator<<(basic_ostream<Ch, Tr>& os, tuple<Args...> const& t) -> basic_ostream<Ch, Tr>& { os << "("; print_tuple(os, t, gen_seq<sizeof...(Args)>()); return os << ")"; }ostream &operator<<(ostream &o, const vvll &v) { rep(i, v.size()) { rep(j, v[i].size()) o << v[i][j] << " "; o << endl; } return o; }template <typename T> ostream &operator<<(ostream &o, const vector<T> &v) { o << '['; rep(i, v.size()) o << v[i] << (i != v.size()-1 ? ", " : ""); o<< "]"; return o; }template <typename T> ostream &operator<<(ostream &o, const set<T> &m) { o << '['; for (auto it = m.begin(); it != m.end(); it++) o << *it << (next(it) != m.end() ? ", " : ""); o << "]"; return o; }template <typename T> ostream &operator<<(ostream &o, const unordered_set<T> &m) { o << '['; for (auto it = m.begin(); it != m.end(); it++) o << *it<< (next(it) != m.end() ? ", " : ""); o << "]"; return o; }template <typename T, typename U> ostream &operator<<(ostream &o, const map<T, U> &m) { o << '['; for (auto it = m.begin(); it != m.end(); it++) o<< *it << (next(it) != m.end() ? ", " : ""); o << "]"; return o; }template <typename T, typename U, typename V> ostream &operator<<(ostream &o, const unordered_map<T, U, V> &m) { o << '['; for (auto it = m.begin();it != m.end(); it++) o << *it; o << "]"; return o; }vector<int> range(const int x, const int y) { vector<int> v(y - x + 1); iota(v.begin(), v.end(), x); return v; }template <typename T> istream& operator>>(istream& i, vector<T>& o) { rep(j, o.size()) i >> o[j]; return i;}string bits_to_string(ll input, ll n=64) { string s; rep(i, n) s += '0' + !!(input & (1ll << i)); reverse(all(s)); return s; }template <typename T> ostream &operator<<(ostream &o, const priority_queue<T> &v) { auto tmp = v; while (tmp.size()) { auto x = tmp.top(); tmp.pop();o << x << " ";} o << endl; return o; }template <typename T> unordered_map<T, ll> counter(vector<T> vec){unordered_map<T, ll> ret; for (auto&& x : vec) ret[x]++; return ret;};string substr(string s, P x) {return s.substr(x.fi, x.se - x.fi); }struct ci : public iterator<forward_iterator_tag, ll> { ll n; ci(const ll n) : n(n) { } bool operator==(const ci& x) { return n == x.n; } booloperator!=(const ci& x) { return !(*this == x); } ci &operator++() { n++; return *this; } ll operator*() const { return n; } };size_t random_seed; namespace std { using argument_type = P; template<> struct hash<argument_type> { size_t operator()(argument_type const& x) const{ size_t seed = random_seed; seed ^= hash<ll>{}(x.fi); seed ^= (hash<ll>{}(x.se) << 1); return seed; } }; }; // hash for various classnamespace myhash{ const int Bsizes[]={3,9,13,17,21,25,29,33,37,41,45,49,53,57,61,65,69,73,77,81}; const int xor_nums[]={0x100007d1,0x5ff049c9,0x14560859,0x07087fef,0x3e277d49,0x4dba1f17,0x709c5988,0x05904258,0x1aa71872,0x238819b3,0x7b002bb7,0x1cf91302,0x0012290a,0x1083576b,0x76473e49,0x3d86295b,0x20536814,0x08634f4d,0x115405e8,0x0e6359f2}; const int hash_key=xor_nums[rand()%20]; const int mod_key=xor_nums[rand()%20]; template<typename T> struct myhash{ std::size_t operator()(const T& val) const { return (hash<T>{}(val)%mod_key)^hash_key; } }; };template <typename T> class uset:public std::unordered_set<T,myhash::myhash<T>> { using SET=std::unordered_set<T,myhash::myhash<T>>; public: uset():SET(){SET::rehash(myhash::Bsizes[rand()%20]);} };uint32_t randxor() { static uint32_t x=1+(uint32_t)random_seed,y=362436069,z=521288629,w=88675123; uint32_t t; t=(x^(x<<11));x=y;y=z;z=w; return( w=(w^(w>>19))^(t^(t>>8)) ); }struct timeval start; double sec() { struct timeval tv; gettimeofday(&tv, NULL); return (tv.tv_sec - start.tv_sec) + (tv.tv_usec - start.tv_usec) *1e-6; }struct init_{init_(){ gettimeofday(&start, NULL); ios::sync_with_stdio(false); cin.tie(0); struct timeval myTime; struct tm *time_st; gettimeofday(&myTime, NULL); time_st = localtime(&myTime.tv_sec); srand(myTime.tv_usec); random_seed = RAND_MAX / 2 + rand() / 2; }} init__;#define rand randxorstatic const long long mo = 1e9+7;#define ldout fixed << setprecision(40)/***********************/// 共通部分/***********************/#define REP(i,n) for(int i=0;i<(int)n;++i)#define FOR(i,c) for(__typeof((c).begin())i=(c).begin();i!=(c).end();++i)#define ALL(c) (c).begin(), (c).end()typedef long long Weight;const Weight INF = 1e18;const Weight EPS = 0; // 浮動小数点なら1e-14struct Edge {ll src, dst;Weight weight; // 最小費用流ではcostの役割ll cap;ll rev; // 残余グラフの対応用bool rev_flag = false; // revなら1Edge(int src = 0, int dst = 0, int weight = 0, Weight cap = 1):src(src), dst(dst), weight(weight), cap(cap){ }};bool operator < (const Edge &e, const Edge &f) {return e.weight != f.weight ? e.weight > f.weight : // !!INVERSE!!e.src != f.src ? e.src < f.src : e.dst < f.dst;}typedef vector<Edge> Edges;typedef vector<Edges> Graph;typedef vector<Weight> Array;typedef vector<Array> Matrix;// 最大流と最小費用流の有向// 無向は自分でひっくり返して追加して下さいvoid addDirected(Graph& g, ll src, ll dst, Weight weight = 0, ll cap = 0) {assert(src < g.size() && src >= 0 && dst < g.size() && dst >= 0);g[src].push_back(Edge(src, dst, weight, cap));}vector<string> names;unordered_map<string, int> name_server;void clearNameServer(void) {names.clear();name_server.clear();}int constructNameServer(vector<string>& names_) {names = names_;for (int i = 0; i < names.size(); i++) {name_server[names[i]] = i;}return name_server.size();}void addDirected(Graph& g, string src, string dst, Weight weight = 0, ll cap = 0) {if (!name_server.count(src)) { cerr << "No such node : " << src << endl; exit(1); }if (!name_server.count(dst)) { cerr << "No such node : " << dst << endl; exit(1); }addDirected(g, name_server[src], name_server[dst], weight, cap);}void addUndirected(Graph& g, ll src, ll dst, Weight weight = 0, ll cap = 0) {addDirected(g, src, dst, weight, cap);addDirected(g, dst, src, weight, cap);}void transformFromMatrixToGraph(Graph& g, Matrix& m) {ll n = m.size();g.resize(n);rep(i, n) rep(j, n) if (m[i][j] != INF) {addDirected(g, i, j, m[i][j]);}}void printGraph(Graph& g) {rep(i, g.size()) {if (!g[i].size())continue;rep(j, g[i].size())cout << "(" << i << ", " << g[i][j].dst << "), ";cout << endl;}}void printGraphCap(Graph& g) {rep(i, g.size()) {if (!g[i].size())continue;rep(j, g[i].size())cout << "(" << i << ", " << g[i][j].dst << " : " << g[i][j].cap << ", " << (g[i][j].rev_flag ? "rev" : "for") << "), ";cout << endl;}}// mode//// 0 : デフォルト。辺だけ表示// 1 : 重みweightだけ表示// 2 : フロー用。weight, capを表示。void vizGraph(Graph& g, int mode = 0, string filename = "out.png") {ofstream ofs("./out.dot");ofs << "digraph graph_name {" << endl;rep(i, g.size()) {if (!g[i].size())continue;rep(j, g[i].size()) {if (name_server.size()) {ofs << " " << names[i] << " -> " << names[g[i][j].dst];} else {ofs << " " << i << " -> " << g[i][j].dst;}if (mode == 1) {ofs << " [ label = \"" << g[i][j].weight << "\"];";} else if (mode == 2) {ofs << " [ label = \"" << g[i][j].weight << "/" << (g[i][j].cap == INF ? "inf" : to_string(g[i][j].cap)) << "\"];";}ofs << endl;}}ofs << "}" << endl;ofs.close();system(((string)"dot -T png out.dot >" + filename).c_str());}// 強連結分解// O(V+E)// const Graph& g 有向グラフ// vector< vector<ll> >& scc グラフの強連結成分void visitStronglyConnectedComponents(const Graph &g, ll v, vector< vector<ll> >& scc,stack<ll> &S, vector<bool> &inS,vector<ll> &low, vector<ll> &num, ll& time) {low[v] = num[v] = ++time;S.push(v); inS[v] = true;FOR(e, g[v]) {ll w = e->dst;if (num[w] == 0) {visitStronglyConnectedComponents(g, w, scc, S, inS, low, num, time);low[v] = min(low[v], low[w]);} else if (inS[w])low[v] = min(low[v], num[w]);}if (low[v] == num[v]) {scc.push_back(vector<ll>());while (1) {ll w = S.top(); S.pop(); inS[w] = false;scc.back().push_back(w);if (v == w) break;}}}void stronglyConnectedComponents(const Graph& g, vector< vector<ll> >& scc) {const ll n = g.size();vector<ll> num(n), low(n);stack<ll> S;vector<bool> inS(n);ll time = 0;rep(u, n) if (num[u] == 0)visitStronglyConnectedComponents(g, u, scc, S, inS, low, num, time);}int main(void) {ll n, m; cin >> n >> m;Graph g = Graph(n);rep(i, m) {ll x, r; cin >> x >> r; x--;rep(i, r) {ll u; cin >> u; u--;addDirected(g, u, x);}}vizGraph(g);vector<vector<ll>> scc;stronglyConnectedComponents(g, scc);vll colored(n);rep(i, scc.size()) if (scc[i].size() != 1)for (auto x : scc[i])colored[x] = 1;function<void(ll)> dfs = [&](ll u) {for (auto x : g[u]) if (!colored[x.dst])colored[x.dst] = 1, dfs(x.dst);};rep(i, n) if (colored[i])dfs(i);ll ret = 0;rep(i, n)ret += colored[i] == 0;cout << ret << endl;return 0;}