#include "shortcut.h"
#include <bits/stdc++.h>
#define endl '\n'
#define here cerr<<"=========================================\n"
#define dbg(x) cerr<<#x<<": "<<x<<endl;
#define ll long long
#define pb push_back
#define popb pop_back
#define all(a_) a_.begin(),a_.end()
#define pll pair<ll,ll>
#define sc second
#define fi first
#define llinf 1000000000000000LL
#define ceri(a,l,r) {cerr<<#a<<": ";for(ll i = l;i<=r;i++) cerr<<a[i]<< " ";cerr<<endl;}
using namespace std;
#define maxn 505
#define lg 9
ll n,c,tsz;
vector<pll> g[maxn];
vector<ll> cyc,ps,curr;
ll cyclen;
bool naso = 0;
bool inc[maxn];
bool vis[maxn];
ll len[maxn][maxn];
ll id[maxn];
ll in[maxn],out[maxn];
ll st[maxn][lg];
ll dept[maxn];
ll mx[maxn][2];
ll ti = 0;
ll cur;
void dfs2(ll u,ll par,ll poc){
id[u] = poc;
st[u][0] = par;
in[u] = ++ti;
if(g[u].size()==1) mx[u][0] = 0;
for(pll p : g[u]){
ll s = p.fi;
ll w = len[u][s];
if(inc[s]) continue;
if(s==par) continue;
dept[s] = dept[u] + w;
dfs2(s,u,poc);
vector<ll> v;
v.pb(mx[s][0] + w);
v.pb(mx[s][0] + w);
v.pb(mx[u][0]);
v.pb(mx[u][1]);
cur = max(cur,mx[u][0] + mx[s][0] + w);
sort(all(v));
mx[u][0] = v[3];
mx[u][1] = v[2];
}
out[u] = ti-1;
}
ll dis[maxn];
void dfs3(ll u,ll par){
for(pll p : g[u]){
ll s = p.fi;
ll w = len[u][s];
if(s==par) continue;
dis[s] = dis[u] + w;
dfs3(s,u);
}
}
ll naj(ll x){
for(ll i = 1;i<=tsz;i++) dis[i] = llinf;
dis[x] = 0;
dfs3(x,x);
ll ans = x;
for(ll i = 1;i<=tsz;i++) if(dis[i]>dis[ans]) ans = i;
return ans;
}
ll find_shortcut(int N, vector<int> L, vector<int> D, int C)
{
n = N;
tsz = n;
c = C;
for(ll i = 0;i<maxn;i++) for(ll j = 0;j<maxn;j++) len[i][j] = llinf;
for(ll i = 1;i<n;i++){
g[i].pb({i+1,L[i-1]});
g[i+1].pb({i,L[i-1]});
len[i][i+1] = len[i+1][i] = L[i-1];
}
for(ll i = 1;i<=n;i++) if(D[i-1]){
g[i].pb({++tsz,D[i-1]});
g[tsz].pb({i,D[i-1]});
len[i][tsz] = len[tsz][i] = D[i-1];
}
ll ans = llinf;
for(ll i = 1;i<=n;i++){
for(ll j = i+2;j<=n;j++){
g[i].pb({j,c});
g[j].pb({i,c});
len[i][j] = len[j][i] = c;
ti = 0;
for(ll k = 1;k<=tsz;k++) vis[k] = inc[k] = id[k] = in[k] = out[k] = dept[k] = 0;
for(ll k = 1;k<=tsz;k++) mx[k][0] = mx[k][1] = -llinf;
cyclen = 0;
cyc.clear();
ps.clear();
ps.pb(0);
for(ll k = i;k<=j;k++) cyc.pb(k),inc[k] = 1;
for(ll k = 1;k<cyc.size();k++) ps.pb(ps[k-1] + len[cyc[k]][cyc[k-1]]);
cyclen = ps.back() + len[cyc[0]][cyc.back()];
cur = 0;
for(ll k = 0;k<cyc.size();k++){
ll y = cyc[k];
st[y][0] = y;
dfs2(y,y,k);
if(g[y].size()==2) mx[y][0] = 0;
}
for(ll e = 0;e<cyc.size();e++){
for(ll f = e+1;f<cyc.size();f++){
ll idx = e,idy = f;
if(idx<idy) swap(idx,idy);
ll r = ps[idx] - ps[idy];
cur = max(cur,mx[cyc[e]][0] + mx[cyc[f]][0] + min(r,cyclen-r));
}
}
ans = min(ans,cur);
len[i][j] = len[j][i] = llinf;
g[i].popb();
g[j].popb();
}
}
for(ll i = 1;i<n;i++){
len[i][i+1] = len[i+1][i] = min((ll)L[i-1],c);
ll x = naj(1);
ll y = naj(x);
ans = min(ans,dis[y]);
len[i][i+1] = len[i+1][i] = L[i-1];
}
return ans;
}
/**
4 10
10 20 20
0 40 0 30
9 30
10 10 10 10 10 10 10 10
20 0 30 0 0 40 0 40 0
3 3
1 1
1 1 1
5 1
1 1 1 1
0 0 0 0 0
**/
