RUPC day3: AA グラフ (AA Graph)

問題概要

https://onlinejudge.u-aizu.ac.jp/#/challenges/sources/VPC/HUPC/2869

アスキーアートでグラフが与えられるので最短経路を求める問題。
個人的にすごく好き。

注意するとこ

　　○｜○○○
　　○｜○G○
　　○｜○○○
この場合はうまく弾かないと行けない。

#include<bits/stdc++.h>
using ll = long long;
#define int ll
#define INF 1e9
#define EPS 0.0000000001
#define rep(i,n) for(int i=0;i<n;i++)
#define all(in) in.begin(), in.end()
int dx[]={0,0,2,-2};
int dy[]={2,-2,0,0};
int _dx[]={0,0,1,-1};
int _dy[]={1,-1,0,0};
int W,H;
bool inrange(int x,int y){return (0<=x&&x<W)&&(0<=y&&y<H);}
using namespace std;
typedef pair<int, int> pii;
#define MP make_pair
signed main(){
    char s,t;
    cin >> H >> W >> s >> t;
    vector<string>grid(H);
    vector<vector<pii>>edge(100);
    map<int,pii>memo;
    rep(i,H) cin >> grid[i];
    rep(i,H)rep(j,grid[i].size()){
        if('A' <= grid[i][j] && grid[i][j] <= 'Z'  ){
            memo[grid[i][j] - 'A'] = MP(j,i);
        }
    }
    for(auto itr : memo){
        int vertex = itr.first;
        pii pos = itr.second;
        rep(i,4){
            pii temp = MP(pos.first+dx[i],pos.second+dy[i]);
            if(!inrange(temp.first,temp.second))continue;
            if(grid[temp.second][temp.first] != '-' && grid[temp.second][temp.first] != '|') continue;
            if(grid[temp.second][temp.first] == '|' && ( i > 1)) continue;
            if(grid[temp.second][temp.first] == '-' && ( i <= 1))continue;
            pii now = temp;
            while(true){
                if('A' <= grid[now.second][now.first] && grid[now.second][now.first] <= 'Z'  )
                    break;
                now.first += _dx[i];
                now.second += _dy[i];
                assert(inrange(now.first, now.second));
            }
            int next =  grid[now.second][now.first] - 'A';
            edge[vertex].push_back(pii(next,1));
        }
    }
    // dijkstra
    vector<int>d(t,INF);
    d[s-'A'] = 0;
    priority_queue<pii,vector<pii>,greater<pii>>que;
    que.push(pii(0,s-'A'));
    while(!que.empty()){
        pii now = que.top();
        que.pop();
        if(now.first > d[now.second])continue;
        rep(i,edge[now.second].size()){
            int next  = edge[now.second][i].first;
            if(d[next] < 1 + now.first) continue;
            d[next] = 1 + now.first;
            que.push(pii(d[next],next));
        }
    }
    cout << d[t-'A'] << endl;
}

CS Academy Array Removal

csacademy.com

問題概要

要素の数列が与えられ、要素のsubarrayが与えられる。
回連続する数列の和の最大値を出力する。ただし、subarrayの番目の数字は使えなくなる。

解く前に考えてたこと

• DP :　不可能そう。使えない数字が出るたびに引くと間に合わない
• StarrySkyTree:区間がどれほど重複してるか数える必要が出てきて、無理となった。

解法

union-findとsubarrayを逆からみていくと良い。
番目の数字を使えるとすると、、番目の要素とuniteしていく。各集合のparentをkeyとする和の配列sumを持っていると良い。
答えとなりうるのは、番目の要素を足す前の最大値か、足した後の和となるので、そのあたりはうまいことやる。

#include <bits/stdc++.h>
using ll = long long;
#define int ll
#define INF 1e9
#define EPS 0.0000000001
#define rep(i,n) for(int i=0;i<n;i++)
#define all(in) in.begin(), in.end()
template<class T, class S> void cmin(T &a, const S &b) { if (a > b)a = b; }
template<class T, class S> void cmax(T &a, const S &b) { if (a < b)a = b; }
using namespace std;
#define MAX 9999999
class unionfind {
    vector<int> par, rank, size_ ,sum;
public:
    unionfind(int n) :par(n), rank(n), size_(n, 1),sum(n,0) {
        iota(all(par), 0);
    }
    int find(int x) {
        if (par[x] == x)return x;
        return par[x] = find(par[x]);
    }
    void unite(int x, int y) {
        
        x = find(x), y = find(y);
        if (x == y)return;
        if (rank[x] < rank[y])swap(x, y);
        par[y] = x;
        size_[x] += size_[y];
        if (rank[x] == rank[y])rank[x]++;
    }
    bool same(int x, int y) {
        return find(x) == find(y);
    }
    int size(int x) {
        return size_[find(x)];
    }
};
signed main(){
    int n, maxi = 0;
    cin >> n;
    vector<int> v(n,0),sub(n),ans,sum(2*n,0);
    vector<bool> used(n,false);
    unionfind uni(2*n);
    rep(i,n) cin >> v[i];
    rep(i,n) cin >> sub[i];
    rep(i,n) {
        uni.unite(i,i+n);
        sum[uni.find(i)] = v[i];
    }
    reverse(all(sub));
    rep(i,sub.size()){
        int index = sub[i]-1;
        maxi = max(maxi,v[index]);
        if(index+1 < n){
            if(used[index+1]){
                int val = sum[uni.find(index)]+sum[uni.find(index+1)];
                uni.unite(index, index+1);
                sum[uni.find(index)] =val;
            }
        }
        if(index-1 >=0){
            if(used[index-1]){
                int val = sum[uni.find(index)]+sum[uni.find(index-1)];
                uni.unite(index,index-1);
                sum[uni.find(index)] = val;
                
            }
        }
        cmax(maxi, sum[uni.find(index)]);
        
        ans.push_back(maxi);
        used[index] = true;
    }
    reverse(all(ans));
    rep(i,ans.size()) cout << ans[i] << endl;
}

RUPC2018 参加記

RUPC2018に参加しました。

Day 0

今回は作問チームに関わってもおらず、何も準備してません。ごめんなさい。

Day 1

家で寝てました。

Day2

・立命館大学　プロジェクト団体 RiPP○が9時に来いって言われたので早起きした。

・久方ぶりにサークルに顔を出すと団体長がｱになってた。

・ @tubo28さん、@ferin_tech15さん,@kjuner8さんとチームを組みました。

・A問題を解いて、その後はE問題をひたすら考えてた。到達可能回数を二分探+シミュレーションしようとしたものの、最善手が選べないことに気づきDijkstraをferinさんと実装するもバグが取れず。悲しかったし申し訳なかったです。

・懇親会にも参加しました。

下の写真はixmel先輩です。

Day3

・集合時間が8時30分は無理です。つらすぎる。

・大悪魔さん(@ukuuku09)といらいざさん（@Eliza_0x)とチームを組みました。

・B問題を通してDのDPをひたすら考えてたんですが、これまた場合分けとかがややこしい事になり通らずに悲しい結果に。

・大悪魔さんが最小カットを✝やるだけ✝といって爆速で通していてプロ。

・実装力を鍛えたいねと思ったコンテストでした

最後に

普段Twitter上で見る方々と実際にあえて楽しかったです。

今回は運営側にもかかわらずほんとに何もしていないので来年はそういったところでも貢献したいね。

あと強くなりたいと思いました。

D- People on a Line

D - People on a Line

コンテスト中に解けなくて悲しい思いをした問題。

コンテスト中に考えてたこと

・全点対最短路を求めることができたら良さそう。
・有効辺のみなので、dfs(コンテスト中はDijkstraをした)時にusedとdを持つと解ける。
（つまりdfsごとに到達と距離の配列を持つ必要がある）
・グラフは非連結の場合がある。↑の方針だとグラフの数だけ配列が必要になって解けそうにない。

正解のときに考えたこと

・1つの座標を決めることで、全ての座標が決まる。
・逆辺を貼ることで、一度のdfsで到達したか、どうかを判定する配列を持つ必要がない。
・これトポロジカル順にソートすると逆辺はらなくてよくない？
・訪問した頂点の全ての辺について探索を行うので、再訪問の必要はない（よって間に合う。

トポロジカルソートでもうまいことやると解けるらしい。

#include<bits/stdc++.h>
using ll = long long;
#define int ll
#define INF 100000000000000LL
#define rep(i,n) for(int i=0;i<n;i++)
using namespace std;
using pii = pair<int,int>;
vector<bool>used;
vector<vector<pii>>edge;
vector<int>d;
bool judge=true;

void dfs(int v,int val){
    used[v]=true;
    rep(i,edge[v].size()){
        int next=edge[v][i].first;
        int cost=edge[v][i].second;
        if(d[next]==INF){
            d[next]=val+cost;
            if(!used[next])dfs(next,d[next]);
        }else if(d[next]!=cost+val){
            judge=false;
        }
    }
}
signed main(){
    int n,m; cin>>n>>m;
    
    edge = vector<vector<pii>>(n);
    used = vector<bool>(n,false);
    d    = vector<int>(n,INF);
    
    rep(i,m){
        int a,b,c;cin>>a>>b>>c;
        --a,--b;
        edge[a].push_back(make_pair(b,c));
        edge[b].push_back(make_pair(a,-c));
    }
    rep(i,n)if(!used[i])dfs(i,0);
    cout << (judge ? "Yes" : "No") << endl;
}

2019/6月追記
コードがひどすぎたので直した。

RiPProを休部する話

RiPProを休部することにしました。

※この記事は完全にポエムです。

理由をメンヘラにならないように簡潔に書く。

 

①留学の準備のために1日3,4時間ほど勉強する必要がでてきた。

留学のために2017年の前半を溶かしたといっても過言ではないので、留学はしっかり結果を残したい。そのためにも、準備をきちんとすることにした。

②競プロが辛い。

*1個人的に競プロの世界は弱肉強食だと思っているので、レートが上がらない、強くなれてる気がしないという状況に、そろそろ辛くなって来た。

精進が足りないと言われれば、それまでなんだけど、*2考えが変わるまで自分のペースでゆったり競プロをやることにした。

来年の4月までいません。

迷惑をかけるかもしれませんが、オネガイシマス。

*1:異論あると思う

*2:競プロは楽しいものだと思ってます

D-塗り絵

　
D: 塗り絵 - AtCoder Beginner Contest 036 | AtCoder

個の島があり、個の橋がある。
島を白または黒に塗っていく。ただし、両端の島を黒で塗ることはできない。
色の塗り方が何通りあるか求める問題です。

グラフ上の頂点を、黒色で塗るか、白色で塗るかをメモ化再帰で求めていきます。

頂点を で塗ったときの場合の数

とします。（1が白、0が黒）

(1)頂点を白色に塗るとき(はの子の数)

の子は白、黒の2つの色に塗ることができるので

(2)頂点を黒色に塗るとき

の子は白色に塗ることができるので

ただし、が葉となる場合はとします。

#include<bits/stdc++.h>
#define vi vector<int>
#define vvi vector<vector<int> >
#define vl vector<ll>
#define vvl vector<vector<ll>>
#define vb vector<bool>
#define vc vector<char>
#define vs vector<string>
using ll = long long;
using ld =long double;
#define int ll
#define INF 1e9
#define EPS 0.0000000001
#define rep(i,n) for(int i=0;i<n;i++)
#define loop(i,s,n) for(int i=s;i<n;i++)
#define all(in) in.begin(), in.end()
template<class T, class S> void cmin(T &a, const S &b) { if (a > b)a = b; }
template<class T, class S> void cmax(T &a, const S &b) { if (a < b)a = b; }
#define MAX 9999999
using namespace std;
typedef pair<int, int> pii;
typedef pair<int,pii> piii;
#define mp make_pair
int MOD = INF+7;
int memo[100005][2];
vector<vector<int>>edge;

int dfs(int v,int from,int clr){
    if(memo[v][clr]!=INF)return memo[v][clr];
    memo[v][clr]=0;
    bool leaf =  true;
    int value=1;
    rep(i,edge[v].size()){
        int nv=edge[v][i];
        if(nv == from)continue;
        leaf = false;
        if(clr){
            value*=(dfs(nv,v,0)+dfs(nv,v,1));
            value%= MOD;
        }else{
            value*=dfs(nv,v,1);
            value%=(MOD);
        }
    }
    if(leaf)memo[v][clr]=1;
    else memo[v][clr]=value;
    return memo[v][clr];
}
signed main(){
    int n;
    cin>>n;
    edge.resize(n);
    rep(i,n-1){
        int a,b; cin>>a>>b;
        --a,--b;
        edge[a].push_back(b);
        edge[b].push_back(a);
    }
    rep(i,n)rep(j,2)memo[i][j]=INF;
    cout<<(dfs(0,0,1)+dfs(0,0,0))%(MOD)<<endl;
    
}

D - Mixing Experiment

D: Mixing Experiment - AtCoder Beginner Contest 054 | AtCoder

薬品をにする最小費用を求める問題。
制約も小さいので、でナップサックをすると解ける。
更新をを薬品の最小の比だと思い込んでて時間を溶かした。
具体的には、漸化式を

 int div=(int)gcd((ll)(j+val[i+1].first),(ll)(k+val[i+1].second));
                pii ni=mp((j+val[i+1].first)/div,(k+val[i+1].second)/div);
                cmin(dp[i+1][ni.first][ni.second],dp[i][j][k]+c[i+1]);

こんな形で更新してた。
頭を冷やしてみると、比の足し算はやばい事がわかるし、比はとの間にしか成り立たない関係ということもわかる。
反省。

#include<bits/stdc++.h>
#define vi vector<int>
#define vvi vector<vector<int> >
#define vl vector<ll>
#define vvl vector<vector<ll>>
#define vb vector<bool>
#define vc vector<char>
#define vs vector<string>
using ll = long long;
using ld =long double;
#define int ll
#define INF 1e9
#define EPS 0.0000000001
#define rep(i,n) for(int i=0;i<n;i++)
#define loop(i,s,n) for(int i=s;i<n;i++)
#define all(in) in.begin(), in.end()
template<class T, class S> void cmin(T &a, const S &b) { if (a > b)a = b; }
template<class T, class S> void cmax(T &a, const S &b) { if (a < b)a = b; }
#define MAX 9999999
using namespace std;
typedef pair<int, int> pii;
typedef pair<int,pii> piii;
#define mp make_pair
long long gcd(long long a, long long b) {return b ? gcd(b, a % b) : a;}
long long lcm(long long m, long long n ){
    if((0 == m )||( 0 == n ))return 0;
    return ((m / gcd(m, n)) * n);
}
int dp[41][401][401];
signed main(){
    int n; cin>>n;
    pii comp; cin>>comp.first>>comp.second;
    rep(i,n)rep(j,401)rep(k,401)dp[i][j][k]=INF;
    dp[0][0][0]=0;
    vector<pii>val;
    vector<int>c(n);
    rep(i,n){
        pii temp; cin>>temp.first>>temp.second;
        cin>>c[i];
        val.push_back(mp(temp.first,temp.second));
        dp[i][val[i].first][val[i].second]=c[i];
    }
    rep(i,n-1){
        rep(j,401){
            rep(k,401){
                if(dp[i][j][k]==INF)continue;
                cmin(dp[i+1][j][k],dp[i][j][k]);
                pii ni=mp((j+val[i+1].first),(k+val[i+1].second));
                cmin(dp[i+1][ni.first][ni.second],dp[i][j][k]+c[i+1]);
            }
        }
    }
    int ans=INF;
    rep(i,401)rep(j,401){
        int div=gcd(i,j);
        if(i&&j){
        if(i/div==comp.first&&j/div==comp.second)cmin(ans,dp[n-1][i][j]);
        }
    }
    if(ans==INF)ans=-1;
    cout<<ans<<endl;
}