結果

問題 No.1288 yuki collection
ユーザー mugen_1337mugen_1337
提出日時 2020-11-28 20:44:20
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 412 ms / 5,000 ms
コード長 4,479 bytes
コンパイル時間 2,745 ms
コンパイル使用メモリ 218,544 KB
実行使用メモリ 4,592 KB
最終ジャッジ日時 2023-10-11 01:59:29
合計ジャッジ時間 10,421 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,352 KB
testcase_01 AC 2 ms
4,352 KB
testcase_02 AC 1 ms
4,352 KB
testcase_03 AC 2 ms
4,352 KB
testcase_04 AC 2 ms
4,356 KB
testcase_05 AC 1 ms
4,352 KB
testcase_06 AC 1 ms
4,352 KB
testcase_07 AC 1 ms
4,352 KB
testcase_08 AC 2 ms
4,352 KB
testcase_09 AC 2 ms
4,352 KB
testcase_10 AC 2 ms
4,352 KB
testcase_11 AC 2 ms
4,352 KB
testcase_12 AC 2 ms
4,352 KB
testcase_13 AC 265 ms
4,368 KB
testcase_14 AC 267 ms
4,352 KB
testcase_15 AC 211 ms
4,352 KB
testcase_16 AC 212 ms
4,352 KB
testcase_17 AC 269 ms
4,356 KB
testcase_18 AC 271 ms
4,376 KB
testcase_19 AC 271 ms
4,448 KB
testcase_20 AC 274 ms
4,352 KB
testcase_21 AC 264 ms
4,376 KB
testcase_22 AC 261 ms
4,384 KB
testcase_23 AC 260 ms
4,500 KB
testcase_24 AC 277 ms
4,468 KB
testcase_25 AC 274 ms
4,348 KB
testcase_26 AC 272 ms
4,436 KB
testcase_27 AC 145 ms
4,396 KB
testcase_28 AC 163 ms
4,372 KB
testcase_29 AC 144 ms
4,356 KB
testcase_30 AC 12 ms
4,488 KB
testcase_31 AC 17 ms
4,408 KB
testcase_32 AC 18 ms
4,484 KB
testcase_33 AC 156 ms
4,372 KB
testcase_34 AC 260 ms
4,372 KB
testcase_35 AC 250 ms
4,352 KB
testcase_36 AC 412 ms
4,592 KB
testcase_37 AC 373 ms
4,352 KB
testcase_38 AC 174 ms
4,372 KB
testcase_39 AC 136 ms
4,352 KB
testcase_40 AC 4 ms
4,480 KB
testcase_41 AC 2 ms
4,348 KB
testcase_42 AC 1 ms
4,356 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
#define ALL(x) begin(x),end(x)
#define rep(i,n) for(int i=0;i<(n);i++)
#define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl;
#define mod 1000000007
using ll=long long;
const int INF=1000000000;
const ll LINF=1001002003004005006ll;
int dx[]={1,0,-1,0},dy[]={0,1,0,-1};
// ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;}
template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;}

#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")

struct IOSetup{
    IOSetup(){
        cin.tie(0);
        ios::sync_with_stdio(0);
        cout<<fixed<<setprecision(12);
    }
} iosetup;
 
template<typename T>
ostream &operator<<(ostream &os,const vector<T>&v){
    for(int i=0;i<(int)v.size();i++) os<<v[i]<<(i+1==(int)v.size()?"":" ");
    return os;
}
template<typename T>
istream &operator>>(istream &is,vector<T>&v){
    for(T &x:v)is>>x;
    return is;
}

template<typename flow_t, typename cost_t>
struct PrimalDual{
    struct edge{
        int to;
        flow_t cap;
        cost_t cost;
        int rev;//この辺の逆辺がg[from]の何番目にあるか
        bool isrev;
    };
 
    vector<vector<edge>> graph;
    vector<cost_t> potential,min_cost;
    vector<int> prevv,preve;//点，辺
    const cost_t TINF;
 
    PrimalDual(int V):graph(V),TINF(numeric_limits<cost_t>::max()){}
 
    void add_edge(int from,int to,flow_t cap,cost_t cost){
        graph[from].push_back((edge){to,cap,cost,(int)graph[to].size(),false});
        graph[to].push_back((edge){from,0,-cost,(int)graph[from].size()-1,true});
    }
 
    cost_t min_cost_flow(int s,int t,flow_t f,bool &ok){
        int V=(int)graph.size();
        cost_t ret=0;
        using Pi=pair<cost_t,int>;
        priority_queue<Pi,vector<Pi>,greater<Pi>> que;
        potential.assign(V,0);
        preve.assign(V,-1);
        prevv.assign(V,-1);
 
        while(f>0){
            min_cost.assign(V,TINF);
            que.emplace(0,s);
            min_cost[s]=0;
            //dijkstraパート
            while(!que.empty()){
                Pi p=que.top();que.pop();
                if(min_cost[p.second]<p.first) continue;
                for(int i=0;i<(int)graph[p.second].size();i++){
                    edge &e=graph[p.second][i];
                    cost_t nextCost=min_cost[p.second]+e.cost+potential[p.second]-potential[e.to];
                    if(e.cap>0 and min_cost[e.to]>nextCost){
                        min_cost[e.to]=nextCost;
                        prevv[e.to]=p.second,preve[e.to]=i;
                        que.emplace(min_cost[e.to],e.to);
                    }
                }
            }
            if(min_cost[t]==TINF){
                ok=false;
                return ret;
            }
            // dijkstraの結果に応じてpotentialを調節
            for(int v=0;v<V;v++)potential[v]+=min_cost[v];
            flow_t addflow=f;
            for(int v=t;v!=s;v=prevv[v]){
                addflow=min(addflow,graph[prevv[v]][preve[v]].cap);
            }
            f-=addflow;
            ret+=addflow*potential[t];
            for(int v=t;v!=s;v=prevv[v]){
                edge &e=graph[prevv[v]][preve[v]];
                e.cap-=addflow;
                graph[v][e.rev].cap+=addflow;
            }
        }
        ok=true;
        return ret;
    }
 
    void output(){
        for(int i=0;i<graph.size();i++){
            for(auto &e:graph[i]){
                if(e.isrev)continue;
                auto &rev_e=graph[e.to][e.rev];
                cout<<i<<"->"<<e.to<<" (flow: "<<rev_e.cap<<" / "<<rev_e.cap+e.cap<<")"<<endl;
            }
        }
    }
};

signed main(){
    int n;cin>>n;
    string s;cin>>s;
    vector<ll> v(n);
    cin>>v;

    PrimalDual<ll,ll> flow(3*n+2);
    int src=3*n,sink=3*n+1;

    rep(i,n-1){
        flow.add_edge(i,i+1,flow.TINF,0);
        flow.add_edge(i+n,i+1+n,flow.TINF,0);
        flow.add_edge(i+2*n,i+1+2*n,flow.TINF,0);
    }

    rep(i,n){
        if(s[i]=='y') flow.add_edge(src,i,1,INF-v[i]);
        if(s[i]=='u') flow.add_edge(i,i+n,1,INF-v[i]);
        if(s[i]=='k') flow.add_edge(i+n,i+2*n,1,INF-v[i]);
        if(s[i]=='i') flow.add_edge(i+2*n,sink,1,INF-v[i]);
    }
    flow.add_edge(src,sink,flow.TINF,INF*4ll);
    bool ok;
    cout<<ll(INF)*4*n-flow.min_cost_flow(src,sink,n,ok)<<endl;
    return 0;
}
0