QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#687776 | #9492. 树上简单求和 | xlwang | 0 | 0ms | 0kb | C++14 | 8.1kb | 2024-10-29 21:10:16 | 2024-10-29 21:10:17 |
answer
#include<bits/stdc++.h>
#define ll long long
#define ull unsigned long long
#define fr(i,j,k) for(register int i=j;i<=k;++i)
#define rf(i,j,k) for(register int i=j;i>=k;--i)
#define foredge(i,j) for(register int i=head[j];i;i=e[i].nxt)
#define pb push_back
#define Times printf("Time:%.3lf\n",clock()/CLOCKS_PER_SEC)
#define pii pair<int,int>
#define mk make_pair
using namespace std;
inline int read(){
int x=0;
bool f=0;
char c=getchar();
while(!isdigit(c)) f|=(c=='-'),c=getchar();
while(isdigit(c)) x=(x<<3)+(x<<1)+(c^48),c=getchar();
return f?-x:x;
}
inline void write(int x){
if(x<0){putchar('-');x=-x;}
if(x>9)write(x/10);
putchar(x%10+'0');
}
inline void writeln(int x){write(x); puts("");}
inline void writepl(int x){write(x); putchar(' ');}
mt19937 rnd(chrono::steady_clock::now().time_since_epoch().count());
inline int randfind(int l,int r){return rnd()%(r-l+1)+l;}
//inline void init(){
// int t=read();
// while(t--) work();
//}
const int Maxn=2e5+10,B=500;
int n,q;
ull a[Maxn];
int toid[Maxn];
struct sgt{
ull s[Maxn<<2],tag[Maxn<<2];
inline int ls(int x){return x<<1;}
inline int rs(int x){return x<<1|1;}
inline void pushup(int x){s[x]=s[ls(x)]+s[rs(x)];}
inline void change(int x,int l,int r,ull k){s[x]+=(r-l+1)*k;tag[x]+=k;}
inline void pushdown(int x,int l,int r){
if(tag[x]){
int mid;mid=(l+r)>>1;
change(ls(x),l,mid,tag[x]);
change(rs(x),mid+1,r,tag[x]);
tag[x]=0;
}
}
inline void update(int ql,int qr,int l,int r,int x,ull k){
// cout<<ql<<' '<<qr<<' '<<l<<' '<<r<<' '<<x<<' '<<k<<endl;
if(ql<=l && r<=qr) {change(x,l,r,k);return;}
pushdown(x,l,r);int mid;mid=(l+r)>>1;
if(ql<=mid) update(ql,qr,l,mid,ls(x),k);
if(mid<qr) update(ql,qr,mid+1,r,rs(x),k);
pushup(x);
}
inline ull query(int ql,int qr,int l,int r,int x){
if(ql<=l && r<=qr) return s[x];
int mid;mid=(l+r)>>1;pushdown(x,l,r);ull ans=0;
if(ql<=mid) ans+=query(ql,qr,l,mid,ls(x));
if(mid<qr) ans+=query(ql,qr,mid+1,r,rs(x));
return ans;
}
inline void clear(int x,int l,int r){
s[x]=tag[x]=0;
if(l==r) return;
int mid;mid=(l+r)>>1;
clear(ls(x),l,mid);clear(rs(x),mid+1,r);
}
}SGT;
struct BLOCK{
ull S[Maxn],s[Maxn];
int id[Maxn];
int L[Maxn],R[Maxn];
inline void into(){
fr(i,1,n){
id[i]=(i-1)/B+1;
if(!L[id[i]]) L[id[i]]=i;
R[id[i]]=i;
}
// fr(i,1,n){
// s[i]=a[toid[i]];
// if(i!=L[id[i]]) s[i]+=s[i-1];
// }
// fr(i,1,n) if(id[i]!=id[n]) S[id[i]+1]+=a[toid[i]];
// fr(i,1,id[n]) S[i]+=S[i-1];
}
inline void update(int x,ull k){
// cout<<x<<' '<<k<<endl;
fr(i,x,R[id[x]]) s[i]+=k;
fr(i,id[x]+1,id[n]) S[i]+=k;
}
inline ull getans(int x){return s[x]+S[id[x]];}
inline ull query(int l,int r){if(l>r) return 0;return getans(r)-getans(l-1);}
}BL;
struct Tree{
vector<int> vc[Maxn];
int f[20][Maxn<<1],dep[Maxn];
int dfn[Maxn],idx;
int num,id[Maxn];
int lg[Maxn<<1];
int father[Maxn];
int siz[Maxn];
inline void dfs(int x,int fa){
dep[x]=dep[fa]+1;f[0][++idx]=x;dfn[x]=idx;father[x]=fa;id[x]=++num;siz[x]=1;
for(auto y:vc[x]) if(y!=fa){
dfs(y,x);
f[0][++idx]=x;siz[x]+=siz[y];
}
}
inline int getmin(int x,int y){if(dep[x]<dep[y]) return x;return y;}
inline int getlca(int x,int y){
x=dfn[x],y=dfn[y];if(x>y) swap(x,y);
int ln=lg[y-x+1];
// cout<<ln<<' '<<x<<' '<<y<<' '<<f[ln][x]<<' '<<f[ln][y-(1<<ln)+1]<<endl;
return getmin(f[ln][x],f[ln][y-(1<<ln)+1]);
}
inline void Read(){
fr(i,1,n-1){
int x,y;
cin>>x>>y;
vc[x].pb(y);vc[y].pb(x);
}
dfs(1,0);
lg[0]=-1;fr(i,1,idx) lg[i]=lg[i/2]+1;
fr(j,1,lg[idx]) fr(i,1,idx){
if(i+(1<<j)-1>idx) break;
f[j][i]=getmin(f[j-1][i],f[j-1][i+(1<<(j-1))]);
// cout<<j<<' '<<i<<' '<<f[j][i]<<endl;
}
}
}T1,T2;
int vis[Maxn];
vector<int> ID;
struct node{int x,y;ull z;}que[Maxn];
ull ans[Maxn];
inline bool cmp(int a,int b){return T2.id[a]<T2.id[b];}
inline void init(){
cin>>n>>q;
fr(i,1,n) cin>>a[i];
T1.Read();T2.Read();
// fr(i,1,n) vis[i]=1;
fr(i,1,n) SGT.update(T2.id[i],T2.id[i]+T2.siz[i]-1,1,n,1,a[i]);
fr(i,1,n) ID.pb(i);
shuffle(ID.begin(),ID.end(),rnd);
fr(i,0,min(n-1,B)) vis[ID[i]]=1;vis[1]=1;
ID.clear();
// fr(i,1,n) vis[i]=1;
fr(i,1,n) if(vis[i]) ID.pb(i);
sort(ID.begin(),ID.end(),cmp);
int ln=ID.size();
fr(i,0,ln-2) ID.pb(T2.getlca(ID[i],ID[i+1]));
sort(ID.begin(),ID.end(),cmp);
ID.erase(unique(ID.begin(),ID.end()),ID.end());
for(auto x:ID) vis[x]=1;
fr(i,1,q) cin>>que[i].x>>que[i].y>>que[i].z;
BL.into();
}
vector<pii> vc[Maxn];
inline void update(int x,int y,ull z){
int llca=T1.getlca(x,y);
// cout<<llca<<endl;
BL.update(T1.id[x],z);BL.update(T1.id[y],z);
BL.update(T1.id[llca],-2*z);
// cout<<"**"<<endl;
SGT.update(T2.id[llca],T2.id[llca]+T2.siz[llca]-1,1,n,1,z);
a[llca]+=z;
}
inline ull query(int x){return BL.query(T1.id[x],T1.id[x]+T1.siz[x]-1);}
inline ull getval(int x,int y){
int llca=T2.getlca(x,y);
ull ans=0;
ans=SGT.query(T2.id[x],T2.id[x],1,n,1);
ans+=SGT.query(T2.id[y],T2.id[y],1,n,1);
ans-=2*SGT.query(T2.id[llca],T2.id[llca],1,n,1);
ans+=a[llca];
return ans;
}
int id1[Maxn];
inline void query(int x,int y,int id){
// cout<<"query:"<<x<<' '<<y<<' '<<id<<endl;
ans[id]=getval(x,y);
// cout<<"ans:"<<ans[id]<<endl;
while(1){
// cout<<x<<' '<<y<<endl;
if(vis[x] && vis[y]) break;
if(x==y){
ans[id]+=query(x),x=T2.father[x];
return;
}
if(vis[y]) ans[id]+=query(x),x=T2.father[x];
else if(vis[x]) ans[id]+=query(y),y=T2.father[y];
else {
if(T2.dep[x]<T2.dep[y]) swap(x,y);
// cout<<x<<endl;
ans[id]+=query(x),x=T2.father[x];
// cout<<"ans:"<<ans[id]<<endl;
}
}
int llca=T2.getlca(x,y);
// cout<<"llca:"<<x<<' '<<y<<' '<<llca<<endl;
vc[x].pb(mk(id,1));vc[y].pb(mk(id,1));vc[llca].pb(mk(id,-2));
id1[id]=llca;
}
inline void Update(int x,int y,ull z){
// cerr<<x<<' '<<y<<' '<<z<<endl;
SGT.update(T1.id[x],T1.id[x],1,n,1,z);
SGT.update(T1.id[y],T1.id[y],1,n,1,z);
int llca=T1.getlca(x,y);
SGT.update(T1.id[llca],T1.id[llca],1,n,1,-2*z);
// cerr<<"end\n";
}
inline ull getpoint(int x){return SGT.query(T1.id[x],T1.id[x]+T1.siz[x]-1,1,n,1)+a[x];}
int fl[Maxn],s[Maxn];
inline void dfs(int x,int fa){
s[x]=s[fa]+fl[x];
for(auto y:T1.vc[x]) if(y!=fa) dfs(y,x);
}
ull Sum=0;
inline void Add(int x,int y,ull z){
int sum=0,llca=T1.getlca(x,y);
sum=s[x]+s[y]-2*s[llca];
Sum+=z*sum;
}
inline void into(int X){
// cout<<"into:"<<X<<endl;
int pre=X;
fr(i,1,n) fl[i]=0;Sum=0;
while(X) fl[X]=1,X=T2.father[X];
dfs(1,0);
int tim=0;
for(auto p:vc[pre]){
fr(i,tim+1,p.first) Add(que[i].x,que[i].y,que[i].z);
ans[p.first]+=Sum*p.second;
tim=p.first;
}
}
inline void work(){
// cout<<"dfn:\n";
// fr(i,1,n) cout<<i<<' '<<T1.id[i]<<endl;
fr(i,1,q) update(que[i].x,que[i].y,que[i].z),query(que[i].x,que[i].y,i);
// cerr<<"**"<<endl;
SGT.clear(1,1,n);fr(i,1,n) a[i]=0;
// cerr<<"**"<<endl;
// cout<<ans[1]<<endl;
fr(i,1,q) Update(que[i].x,que[i].y,que[i].z),ans[i]+=getpoint(id1[i]);
// cout<<ans[1]<<endl;
for(auto x:ID) into(x);
fr(i,1,q) cout<<ans[i]<<'\n';
}
signed main(){
// freopen("input.in","r",stdin);
// freopen("output.out","w",stdout);
ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
init();work();
// printf("\nTIME:%.3lf",(double)clock()/CLOCKS_PER_SEC);
return 0;
}
詳細信息
Subtask #1:
score: 0
Runtime Error
Test #1:
score: 0
Runtime Error
input:
3000 3000 7236742292501328495 17973811477309806363 16075782662531676171 17971236571771878676 11392080645527132110 3685563455925680459 9773593720088356683 8313828403245053795 7736401634567449043 1634817828009987181 6951124933529719486 12775126714635387213 15460977209223753216 397573676785925632 31372...
output:
result:
Subtask #2:
score: 0
Skipped
Dependency #1:
0%
Subtask #3:
score: 0
Skipped
Dependency #2:
0%
Subtask #4:
score: 0
Runtime Error
Test #21:
score: 0
Runtime Error
input:
200000 200000 622783158027686223 2242697872372232537 8481648430436878777 10092474834140799044 15403999682625301609 12614289513474949582 9180944589267018841 7823784919308285798 8257785171198951273 5134508521895120821 8041682272181381093 3835432206618893170 2653803171409877650 5589823419153460372 1007...
output:
result:
Subtask #5:
score: 0
Runtime Error
Test #27:
score: 0
Runtime Error
input:
200000 200000 1958469220619413759 14991498002015735322 6054491201406941902 18206143187746582567 15082377615826460430 2936248617457291604 10073577150351675920 16534472678586906457 2207599132486246393 10301540360769075442 1492580560381080472 551692353431379140 13238280352539145808 8462626987240986565 ...
output:
result:
Subtask #6:
score: 0
Runtime Error
Test #34:
score: 0
Runtime Error
input:
200000 200000 6794776813641982926 1561596256197101737 10910039723053043515 7892247858295192798 12233819960547881004 17695389034783066733 9173201689566865598 17626618141377486739 7358781671024283919 6787559733384974662 3884392438269280436 14872846228351316833 9037842441501571648 14299818404271084016 ...
output:
result:
Subtask #7:
score: 0
Skipped
Dependency #1:
0%