#include using namespace std; #define REP(i,a,n) for(int i=(a); i<(int)(n); i++) #define rep(i,n) REP(i,0,n) #define FOR(it,c) for(__typeof((c).begin()) it=(c).begin(); it!=(c).end(); ++it) #define ALLOF(c) (c).begin(), (c).end() typedef long long ll; typedef unsigned long long ull; class Dinic { int MAX_V; int INF; struct edge{ int to, cap, rev, icap, flow; }; vector< vector > G; vector level; //sからの距離 vector iter; //どこまで調べたか void max_flow_bfs(int s){ fill(level.begin(), level.end(), -1); queue que; level[s] = 0; que.push(s); while(!que.empty()){ int v = que.front(); que.pop(); for(int i=0; i0 && level[e.to]<0){ level[e.to] = level[v] + 1; que.push(e.to); } } } } int max_flow_dfs(int v, int t, int f){ if(v==t) return f; for(int &i=iter[v]; i0 && level[v]0){ e.cap -= d; G[e.to][e.rev].cap += d; e.flow += d; return d; } } } return 0; } public: Dinic(int N):MAX_V(N),G(N),level(N),iter(N){ INF = 99999999; } void add_edge(int from, int to, int cap){ G[from].push_back((edge){to, cap, (int)G[to].size(), cap, 0}); G[to].push_back((edge){from, 0, (int)G[from].size()-1, 0, 0}); } int get_flow(int from, int to){ //untried rep(i,G[from].size()){ if(G[from][i].to == to){ return G[from][i].flow; } } return -1; } int max_flow(int s, int t){ int flow = 0; while(true){ max_flow_bfs(s); if(level[t]<0) return flow; fill(iter.begin(), iter.end(), 0); int f; while((f = max_flow_dfs(s, t, INF))>0){ flow += f; } } } int min_cut(int s, int t, vector& S, vector& T){ S.clear(); T.clear(); int maxf = max_flow(s, t); for(int i=0; i= 0) S.push_back(i); else T.push_back(i); } return maxf; } }; template class MinCostFlow { struct Edge { int rev, from, to; F cap, icap; C cost; Edge(int rev, int from, int to, F cap, C cost): rev(rev), from(from), to(to), cap(cap), icap(cap), cost(cost){} }; int N; vector> G; const C INF; public: MinCostFlow(int N):N(N),G(N),INF(numeric_limits::max()){} void add_edge(int from, int to, F cap, C cost){ G[from].emplace_back((int)(G[to].size()), from, to, cap, cost); G[to].emplace_back((int)(G[from].size()) - 1, to, from, 0, -cost); } C solve(int s, int t, F init_f){ vector dist(N); vector prevv(N); vector preve(N); C ret = 0; F f = init_f; while(f > 0){ fill(dist.begin(), dist.end(), INF); dist[s] = 0; while(true){ bool update = false; for(int v=0; v 0 && dist[e.to] > dist[v] + e.cost) { dist[e.to] = dist[v] + e.cost; prevv[e.to] = v; preve[e.to] = i; update = true; } } } if(!update) break; } if(dist[t] == INF) return 0; F d = f; for(int v=t; v!=s; v=prevv[v]){ d = min(d, G[prevv[v]][preve[v]].cap); } f -= d; ret += dist[t] * d; for(int v=t; v!=s; v=prevv[v]){ Edge& e = G[prevv[v]][preve[v]]; Edge& re = (e.from != e.to) ? G[e.to][e.rev] : G[e.to][e.rev+1]; e.cap -= d; re.cap += d; } } return ret; } vector get_selected_edges(int i){ vector ret; for(int j=0; j> N; string S; cin >> S; vector v; rep(i,N){ ll a; cin >> a; v.push_back(a); } Dinic dinic(2*N+2); MinCostFlow mcf(2*N+2); int s = 2*N; int t = s+1; rep(i,N){ if(S[i] == 'y'){ mcf.add_edge(s, i, 1, -v[i]); dinic.add_edge(s, i, 1); } mcf.add_edge(i, i+N, 1, 0); dinic.add_edge(i, i+N, 1); if(S[i] == 'i'){ mcf.add_edge(i+N, t, 1, 0); dinic.add_edge(i+N, t, 1); } rep(j,N){ if(S[i] == 'y' && S[j] == 'u'){ if(i