提交记录 9579


用户 题目 状态 得分 用时 内存 语言 代码长度
zrz0520 noip18f. 【NOIP2018】保卫王国 Wrong Answer 0 2 s 588896 KB C++ 4.67 KB
提交时间 评测时间
2019-06-10 11:09:01 2020-08-01 01:40:39
#pragma GCC optimize(3)
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+5;
const ll INF=1e11;
int Head[N],Next[N<<1],vet[N<<1];
int top[N],fa[N],ch[N][2],son[N],sz[N],lsz[N],stk[N],Fa[N],pre[N];
ll dp[N][2],g[N][2],p[N];
int n,m,u,v,edgenum,a,x,b,y,root,cnt;
char s[10];
inline ll Min(const ll& a, const ll& b){ return a<b?a:b; }
struct Matrix{ ll a[2][2]; };
inline Matrix Get(const ll& x, const ll& y){
    Matrix c;
    c.a[0][0]=INF; c.a[0][1]=x;
    c.a[1][0]=c.a[1][1]=y;
    return c;
}
Matrix tree[N],Node[N];
Matrix operator * (const Matrix& a, const Matrix& b){
    Matrix c;
    c.a[0][0]=Min(a.a[0][0]+b.a[0][0],a.a[0][1]+b.a[1][0]);
    c.a[0][1]=Min(a.a[0][0]+b.a[0][1],a.a[0][1]+b.a[1][1]);
    c.a[1][0]=Min(a.a[1][0]+b.a[0][0],a.a[1][1]+b.a[1][0]);
    c.a[1][1]=Min(a.a[1][0]+b.a[0][1],a.a[1][1]+b.a[1][1]);
    return c;
}
inline void addedge(const int& u, const int& v){
    vet[++edgenum]=v;
    Next[edgenum]=Head[u];
    Head[u]=edgenum;
}
// inline void pushup(const int& u){ tree[u]=tree[ch[u][0]]*Node[u]*tree[ch[u][1]]; }
inline void update(int u){
    Node[u]=Get(g[u][0],g[u][1]);
    for (; u; u=fa[u]){
        int f=fa[u];
        if (f && ch[f][0]!=u && ch[f][1]!=u){
            g[f][0]-=tree[u].a[1][1]; g[f][1]-=Min(tree[u].a[0][1],tree[u].a[1][1]);
            // pushup(u);
            tree[u]=tree[ch[u][0]]*Node[u]*tree[ch[u][1]];
            g[f][0]+=tree[u].a[1][1]; g[f][1]+=Min(tree[u].a[0][1],tree[u].a[1][1]);
            Node[f]=Get(g[f][0],g[f][1]);
        } else tree[u]=tree[ch[u][0]]*Node[u]*tree[ch[u][1]];//pushup(u);
    }
}
int build(int l, int r){
    if (l>r) return 0;
    int hav=(pre[r]+pre[l-1]+1)/2;
    int k=lower_bound(pre+l,pre+r+1,hav)-pre,u=stk[k];
    fa[ch[u][0]=build(l,k-1)]=u;
    fa[ch[u][1]=build(k+1,r)]=u;
    tree[u]=tree[ch[u][0]]*Node[u]*tree[ch[u][1]];//pushup(u);
    return u;
}
void dfs1(const int& u, const int& f){
    sz[u]=1; dp[u][1]=p[u]; fa[u]=Fa[u]=f;
    for (int e=Head[u]; e; e=Next[e]){
        int v=vet[e];
        if (v==f) continue;
        dfs1(v,u); sz[u]+=sz[v];
        dp[u][0]+=dp[v][1];
        dp[u][1]+=Min(dp[v][0],dp[v][1]);
        if (sz[v]>sz[son[u]]) son[u]=v;
    }
    lsz[u]=sz[u]-sz[son[u]];
}
void dfs2(const int& u, const int& f, const int& st){
    top[u]=st; g[u][1]=p[u];
    if (son[u]) dfs2(son[u],u,st);
    for (int e=Head[u]; e; e=Next[e]){
        int v=vet[e];
        if (v==son[u] || v==f) continue;
        dfs2(v,u,v);
        g[u][0]+=dp[v][1];
        g[u][1]+=Min(dp[v][0],dp[v][1]);
    }
    Node[u]=Get(g[u][0],g[u][1]);
    if (top[u]==u){
        cnt=0;
        for (int i=u; i; i=son[i]) stk[++cnt]=i,pre[cnt]=pre[cnt-1]+lsz[i];
        int las=fa[u];
        fa[root=build(1,cnt)]=las;
    }
}
/*
inline int read(){
    int num=0,fu=1; char ch=getchar();
    while (ch<'0' || ch>'9'){ if (ch=='-') fu=0; ch=getchar(); }
    while (ch>='0' && ch<='9'){ num=(num<<3)+(num<<1)+ch-'0'; ch=getchar(); }
    return fu?num:-num;
}
*/
inline char read() {
	static const int IN_LEN = 1000000;
	static char buf[IN_LEN], *s, *t;
	return (s==t?t=(s=buf)+fread(buf,1,IN_LEN,stdin),(s==t?-1:*s++):*s++);
}
template<class T>
inline void read(T &x) {
	static bool iosig;
	static char c;
	for (iosig=false, c=read(); !isdigit(c); c=read()) {
		if (c == '-') iosig=true;
		if (c == -1) return;
	}
	for (x=0; isdigit(c); c=read()) x=((x+(x<<2))<<1)+(c^'0');
	if (iosig) x=-x;
}
const int OUT_LEN = 10000000;
char obuf[OUT_LEN], *ooh=obuf;
inline void print(char c) {
	if (ooh==obuf+OUT_LEN) fwrite(obuf, 1, OUT_LEN, stdout), ooh=obuf;
	*ooh++=c;
}
template<class T>
inline void print(T x) {
	static int buf[30], cnt;
	if (x==0) print('0');
	else {
		if (x<0) print('-'), x=-x;
		for (cnt=0; x; x/=10) buf[++cnt]=x%10+48;
		while(cnt) print((char)buf[cnt--]);
	}
}
inline void flush() { fwrite(obuf, 1, ooh - obuf, stdout); }

int main(){
    freopen("defense.in","r",stdin);
    freopen("defense.out","w",stdout);
    read(n); read(m); scanf("%s",s);
    ll *P=p;
    for (register int i=1; i<=n; i++) read(*(P+i));
    for (register int i=1; i<n; i++){
        read(u); read(v);
        addedge(u,v); addedge(v,u);
    }
    tree[0].a[0][0]=tree[0].a[1][1]=0;
    tree[0].a[0][1]=tree[0].a[1][0]=INF;
    dfs1(1,0); dfs2(1,0,1);
    for (; m; m--){
        read(a); read(x); read(b); read(y);
        if (!x && !y && (Fa[a]==b || Fa[b]==a)){ puts("-1"); continue; }
        x^=1; y^=1;
        ll last1=g[a][x],last2=g[b][y];
        g[a][x]=INF; update(a);
        g[b][y]=INF; update(b);
        print(Min(tree[root].a[0][1],tree[root].a[1][1]));
        g[a][x]=last1; update(a);
        g[b][y]=last2; update(b);
    }
    flush();
    return 0;
}
/*
5 3 C3
2 4 1 3 9
1 5
5 2
5 3
3 4
1 0 3 0
2 1 3 1
1 0 5 0
*/

CompilationN/AN/ACompile OKScore: N/A

Testcase #146.44 us156 KBWrong AnswerScore: 0

Testcase #252.16 us156 KBWrong AnswerScore: 0

Testcase #350.67 us116 KBWrong AnswerScore: 0

Testcase #444.41 us116 KBWrong AnswerScore: 0

Testcase #554.03 us172 KBWrong AnswerScore: 0

Testcase #651.3 us172 KBWrong AnswerScore: 0

Testcase #761.04 us148 KBWrong AnswerScore: 0

Testcase #8196.03 us296 KBWrong AnswerScore: 0

Testcase #9194.68 us296 KBWrong AnswerScore: 0

Testcase #10426.21 us628 KBWrong AnswerScore: 0

Testcase #11425.97 us628 KBWrong AnswerScore: 0

Testcase #1267.126 ms575 MB + 96 KBRuntime ErrorScore: 0

Testcase #1367.693 ms575 MB + 96 KBRuntime ErrorScore: 0

Testcase #1455.62 ms26 MB + 792 KBWrong AnswerScore: 0

Testcase #1567.798 ms575 MB + 96 KBRuntime ErrorScore: 0

Testcase #1667.774 ms575 MB + 96 KBRuntime ErrorScore: 0

Testcase #1767.842 ms575 MB + 96 KBRuntime ErrorScore: 0

Testcase #1826.559 ms17 MB + 808 KBWrong AnswerScore: 0

Testcase #192 s8 MB + 4 KBTime Limit ExceededScore: 0

Testcase #202 s12 MB + 152 KBTime Limit ExceededScore: 0

Testcase #2167.55 ms25 MB + 440 KBWrong AnswerScore: 0

Testcase #2220.979 ms23 MB + 496 KBWrong AnswerScore: 0

Testcase #2320.995 ms23 MB + 484 KBWrong AnswerScore: 0

Testcase #2418.395 ms21 MB + 676 KBWrong AnswerScore: 0

Testcase #252 s555 MB + 1008 KBTime Limit ExceededScore: 0


Judge Duck Online | 评测鸭在线
Server Time: 2024-04-24 04:50:06 | Loaded in 2 ms | Server Status
个人娱乐项目,仅供学习交流使用