結果

問題 No.177 制作進行の宮森あおいです！
コンテスト
ユーザー tubo28
提出日時 2015-09-24 18:26:52
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 77 ms / 2,000 ms
コード長 3,995 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,055 ms
コンパイル使用メモリ 174,040 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-19 09:00:54
合計ジャッジ時間 2,389 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 13
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#define _CRT_SECURE_NO_WARNINGS
// #define _GLIBCXX_DEBUG
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
// #define int ll
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef pair<int,int> pii;
#define all(c) begin(c), end(c)
#define range(i,a,b) for(ll i=a; i<ll(b); i++)
#define rep(i,b) range(i,0,b)
#define rangei(i,a,b) for(ll a=a;i<=ll(b);i++)
#define repi(i,b) rangei(i,1,b)
#define pb push_back
#define eb emplace_back
#define mp make_pair
#define mt make_tuple
template<class T> ostream & operator << (ostream &os, vector<T> const &);
template<int n, class...T>
typename enable_if<(n>=sizeof...(T))>::type
_ot(ostream &, tuple<T...> const &){}
template<int n, class...T>
typename enable_if<(n< sizeof...(T))>::type
_ot(ostream &os, tuple<T...> const &t){
    os << (n==0?"":" ") << get<n>(t); _ot<n+1>(os, t);
}
template<class...T>
ostream & operator << (ostream &os, tuple<T...> const &t){
    _ot<0>(os, t); return os;
}
template<class T, class U>
ostream & operator<<(ostream &os, pair<T,U> const &p){
    return os << "(" << p.first << ", " << p.second << ") ";
}
template<class T>
ostream & operator<<(ostream &os, vector<T> const &v){
    rep(i,v.size()) os << v[i] << (i+1==(int)v.size()?"":" "); return os;
}
#ifdef DEBUG
#define dump(...) (cerr << #__VA_ARGS__ << " = " << mt(__VA_ARGS__) \
                   << " [" << __LINE__ << "]" << endl)
#else
#define dump(...)
#endif
void fastios(){
    ios_base::sync_with_stdio(0);
    cin.tie(0);
    // #define endl '\n'
}
template<class T>
size_t uniq(vector<T> &v){
    sort(v.begin(), v.end());
    v.erase(unique(v.begin(), v.end()), v.end());
    return v.size();
}
template<class T>
size_t uniq(T *l, size_t n){
    sort(l,l+n);
    return unique(l,l+n) - l;
}
#define mems(arr,val) memset(arr,val,sizeof(arr));
int const mod = 1000000007;
int const inf = numeric_limits<int>::max()/8;

struct Edge {
    int src, dst;
    int cap;
    Edge(int src_, int dst_, int cap_) :
        src(src_), dst(dst_), cap(cap_) { }
};
typedef vector<Edge> Edges;
typedef vector<Edges> Graph;
typedef vector<int> Array;
typedef vector<Array> Matrix;

int dfs(int v, int t, int lim, Matrix &cap, vector<bool> &vis){
    vis[v] = true;
    const size_t n = cap.size();
    if(v == t) return lim;
    for(size_t i = 0; i < n; i++){
        if(vis[i] || cap[v][i] == 0) continue;
        int f = dfs(i, t, min(lim, cap[v][i]), cap, vis);
        if(f > 0){
            cap[v][i] -= f;
            cap[i][v] += f;
            return f;
        }
    }
    return 0;
}

int max_flow(int s, int t, Matrix &cap){
    vector<bool> vis;
    int n = cap.size();
    int maxflow = 0;
    while(1){
        vis.assign(n,false);
        int f = dfs(s,t,inf,cap,vis);
        if(f == 0) break;
        maxflow += f;
    }
    return maxflow;
}

int max_flow(int s, int t, const Graph &g){
    int n = g.size();
    Matrix cap(n, Array(n));
    for(size_t i = 0; i < g.size(); i++){
        for(size_t j = 0; j < g[i].size(); j++){
            cap[g[i][j].src][g[i][j].dst] = g[i][j].cap;
        }
    }
    rep(i,n) dump(cap[i]);
    int res = max_flow(s, t, cap);
    rep(i,n) dump(cap[i]);
    return res;
}

ll W,N;
ll J[60];
ll M;
ll C[60];
bool ok[60][60];

int main(){
    while(cin >> W){
        cin >> N;
        rep(i,N) cin >> J[i];
        cin >> M;
        rep(i,M) cin >> C[i];
        rep(i,N)rep(j,M) ok[i][j] = true;
        rep(i,M){
            int Q; cin >> Q;
            rep(j,Q){
                int X; cin >> X;
                X--;
                ok[X][i] = false;
            }
        }

        int S = N+M, T = N+M+1;
        Graph g(N+M+2);
        rep(i,N){
            g[S].eb(S,i,J[i]);
        }
        rep(i,N)rep(j,M){
            if(ok[i][j]){
                g[i].eb(i,N+j,min(J[i],C[j]));
            }
        }
        rep(i,M){
            g[i+N].eb(i+N,T,C[i]);
        }
        ll f = max_flow(S,T,g);
        dump(f);
        puts(f >= W ? "SHIROBAKO" : "BANSAKUTSUKITA");
    }
}
0