結果

問題 No.1288 yuki collection
ユーザー mugen_1337mugen_1337
提出日時 2020-11-28 20:38:53
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 4,450 bytes
コンパイル時間 2,708 ms
コンパイル使用メモリ 218,928 KB
実行使用メモリ 4,592 KB
最終ジャッジ日時 2023-10-11 01:57:36
合計ジャッジ時間 15,418 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,352 KB
testcase_01 AC 2 ms
4,352 KB
testcase_02 AC 2 ms
4,348 KB
testcase_03 AC 2 ms
4,356 KB
testcase_04 AC 1 ms
4,356 KB
testcase_05 AC 2 ms
4,352 KB
testcase_06 AC 1 ms
4,348 KB
testcase_07 AC 2 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,348 KB
testcase_13 AC 289 ms
4,376 KB
testcase_14 AC 291 ms
4,352 KB
testcase_15 AC 232 ms
4,364 KB
testcase_16 AC 231 ms
4,396 KB
testcase_17 AC 293 ms
4,584 KB
testcase_18 AC 298 ms
4,580 KB
testcase_19 AC 298 ms
4,380 KB
testcase_20 AC 298 ms
4,376 KB
testcase_21 AC 291 ms
4,352 KB
testcase_22 AC 287 ms
4,372 KB
testcase_23 AC 285 ms
4,356 KB
testcase_24 AC 303 ms
4,456 KB
testcase_25 AC 297 ms
4,348 KB
testcase_26 AC 298 ms
4,372 KB
testcase_27 AC 163 ms
4,396 KB
testcase_28 AC 181 ms
4,592 KB
testcase_29 AC 160 ms
4,452 KB
testcase_30 AC 14 ms
4,480 KB
testcase_31 AC 19 ms
4,404 KB
testcase_32 AC 20 ms
4,352 KB
testcase_33 TLE -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
testcase_41 -- -
testcase_42 -- -
権限があれば一括ダウンロードができます

ソースコード

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,-v[i]);
        if(s[i]=='u') flow.add_edge(i,i+n,1,-v[i]);
        if(s[i]=='k') flow.add_edge(i+n,i+2*n,1,-v[i]);
        if(s[i]=='i') flow.add_edge(i+2*n,sink,1,-v[i]);
    }
    flow.add_edge(src,sink,flow.TINF,0);
    bool ok;
    cout<<-flow.min_cost_flow(src,sink,n,ok)<<endl;
    return 0;
}
0