\[だから妄想感傷代償連盟
\]
\[愛を懐いて理想を叫んだ
\]
\[行き場のない愚者のメロディー
\]
\[再挑戦•転生•テレポーテーション
\]
\[何回だって 重ねて逝くんだ
\]
\[終わりなき愛の隨に さあ
\]
\[愛や厭...
\]
提要：打比賽的時候這東西在我腦子裡放了四個小時
A.黎明與螢火
容易想到無解的判斷：每次刪除節點都會至少導致減少 \(4\) 的總度數，且減少的度數總會是 \(4\) 的倍數，因此總度數不為 \(4\) 的倍數的即為無解，搜一遍即可.
隨後輸出方案數. 考慮到葉節點只能被其父節點刪除導致刪除，因此可以先從上到下刪除葉節點的父節點，再從上到下刪除葉節點即可. 實際實現的時候比較簡單，兩遍搜即可.
複雜度 \(O(N)\)
#include<bits/stdc++.h>
using namespace std;
int n;
vector<int>e[200001];
vector<int>son[200001];
int degree[200001];
bool vis[200001];
int fa[200001];
void dfs(int now,int last){
fa[now]=last;
son[last].push_back(now);
for(int i:e[now]){
if(i!=last) dfs(i,now);
}
}
void print1(int now){
if(vis[now]){
for(int i:son[now]) print1(i);
return;
}
for(int i:son[now]){
print1(i);
}
if(degree[now]%2==0){
vis[now]=true;
printf("%d\n",now);
degree[fa[now]]--;
for(int i:son[now]){
degree[i]--;
}
}
}
void print2(int now){
if(vis[now]){
for(int i:son[now]) print2(i);
return;
}
if(degree[now]%2==0){
vis[now]=true;
printf("%d\n",now);
degree[fa[now]]--;
for(int i:son[now]){
degree[i]--;
}
}
for(int i:son[now]){
print2(i);
}
}
int main(){
// freopen("test.in","r",stdin);
// freopen("test.out","w",stdout);
scanf("%d",&n);int root=0;
for(int i=1;i<=n;++i){
int x;scanf("%d",&x);
if(x){
e[i].push_back(x);
e[x].push_back(i);
}
}
for(int i=1;i<=n;++i){
if(e[i].size()%2==0){
root=i;
break;
}
}
dfs(root,0);
int res=0;
for(int i=1;i<=n;++i){
res+=son[i].size();
degree[i]=son[i].size()+1;
}
degree[root]--;
if((res+n-1)%4!=0){
cout<<"NO"<<'\n';
return 0;
}
cout<<"YES"<<'\n';
print1(root);
print2(root);
}
B.Darling Dance
這個 Darling Dance 曲繪怎麼一股兔子洞味，但是完全不是同一個 P 主
好題
考慮先找出哪些邊會在一些節點到 \(1\) 的最短路上，可以發現這些邊會構成一顆以 \(1\) 為根的樹
證明：在最短路上走環一定不優，因此不存在環，故為樹
又因為該樹有全部的 \(n\) 個節點，可以推至樹上共 \(n-1\) 條邊.
因此當 \(k\ge n-1\) 時，直接全選即可.
否則，考慮到不選邊緣的邊可以使答案更優，直接從根節點搜下來即可
#include<bits/stdc++.h>
using namespace std;
#define int long long
int n,m,k;
struct edge{
int to,w,id;
};
vector<edge>e[300001];
int dis[300001];
bool vis[300001];
int pre[300001];
vector<int>e2[300001];
struct node{
int id,dis;
bool operator <(const node &A)const{
return dis>A.dis;
}
};
priority_queue<node>p;
bool istree[300001];
void dij(int s){
memset(vis,0,sizeof vis);
memset(dis,0x3f,sizeof dis);
dis[s]=0;
p.push({s,dis[s]});
while(!p.empty()){
node u=p.top();
p.pop();
if(vis[u.id]) continue;
istree[pre[u.id]]=true;
vis[u.id]=true;
for(edge i:e[u.id]){
if(!vis[i.to]){
if(dis[i.to]>dis[u.id]+i.w){
dis[i.to]=dis[u.id]+i.w;
p.push({i.to,dis[i.to]});
pre[i.to]=i.id;
}
}
}
}
}
vector<int>ans;
void dfs(int now,int last){
if((int)ans.size()==k) return;
// cout<<"dfs "<<now<<" "<<last<<endl;
for(edge i:e[now]){
if((int)ans.size()==k) return;
if(i.to!=last and istree[i.id]){
ans.push_back(i.id);
dfs(i.to,now);
}
}
}
signed main(){
// freopen("T2.in","r",stdin);
// freopen("out.out","w",stdout);
scanf("%lld %lld %lld",&n,&m,&k);
for(int i=1;i<=m;++i){
int x,y,z;
scanf("%lld %lld %lld",&x,&y,&z);
e[x].push_back({y,z,i});
e[y].push_back({x,z,i});
}
dij(1);
// for(int i=1;i<=m;++i){
// cout<<i<<" "<<istree[i]<<endl;
// }
dfs(1,0);
cout<<ans.size()<<endl;
for(int i:ans){
cout<<i<<" ";
}
}
C.Non-breath oblige
可以把操作離線下來搞掃描線，但是我不會，所以學了另一種方法
首先可以想到操作能直接用暴力線段樹維護
1 操作就相當於兩次單點查詢於修改
2 操作是區間賦值
3 操作是單點查詢
因此可以直接建線段樹，但是複雜度太高，考慮最佳化。
注意到可以直接放棄暴力修改，直接記下操作一改變後的 \(pos\)，然後每次查詢直接跳到對應區間查詢，這是一種可行的思路，提議轉化成詢問在 \(l\) 到 \(r\) 之間 \(type_i=3\) 且 \(pos_i\geq l\) 的所有 \(i\) 的 \(val_i\) 之和。
接下來可以上一顆樹狀陣列，維護操作變化，把詢問按 \(l\) 降序排序，然後按 \(pos\) 降序依次把 \(val\) 加入樹狀陣列，拿雙指標搞，詢問的時候直接查 \(query(r)-query(l-1)\) 就行了。
#include<bits/stdc++.h>
using namespace std;
template<typename T>
inline void read(T& x){
x=0;bool sym=0;char c=getchar();
while(!isdigit(c)){sym^=(c=='-');c=getchar();}
while(isdigit(c)){x=x*10+c-48;c=getchar();}
if(sym)x=-x;
}
template<size_t N>
inline void read(char (&str)[N]){
size_t n=0;char c=getchar();
while(n<N-1&&!isspace(c)){str[n]=c;c=getchar();++n;}
str[n]=0;
}
template<typename T,size_t N>
inline void read(T (&a)[N],int range=N){
for(int i=1;i<=range;++i){read(a[i]);}
}
template<typename T,typename... Args>
inline void read(T& x,Args&... args){
read(x);read(args...);
}
template<typename T,typename T2>
inline void readarray(T& x,T2& args){
read(x);read(args,x);
}
template<typename func,typename... Args>
inline void readact(int x,function<func>fu,Args&... args){
for(int i=1;i<=x;++i){
read(args...);
fu(args...);
}
}
#define inread(x) int (x);read(x)
inline void write(int A){if(A<0){putchar('-');A=-A;}if(A>9){write(A/10);}putchar(A%10+'0');}
inline void write(long long A){if(A<0){putchar('-');A=-A;}if(A>9){write(A/10);}putchar(A%10+'0');}
inline void write(char A){putchar(A);}
int n,m,Q;
struct tree{
int l,r;
int w;
}t[4000001];
#define tol (id*2)
#define tor (id*2+1)
#define mid(l,r) mid=((l)+(r))/2
inline void pushdown(int id){
if(t[id].w){
t[tol].w=t[tor].w=t[id].w;
t[id].w=0;
}
}
void build(int id,int l,int r){
t[id].l=l;t[id].r=r;
if(l==r){
t[id].w=0;
return;
}
int mid(l,r);
build(tol,l,mid);
build(tor,mid+1,r);
}
void change(int id,int l,int r,int val){
// cout<<id<<" "<<l<<" "<<r<<" "<<t[id].l<<" "<<t[id].r<<" "<<val<<endl;
if(l<=t[id].l and t[id].r<=r){
t[id].w=val;
return;
}
pushdown(id);
if(l<=t[tol].r) change(tol,l,r,val);
if(t[tor].l<=r) change(tor,l,r,val);
}
int ask(int id,int pos){
if(t[id].l==t[id].r or t[id].w) return t[id].w;
int mid(t[id].l,t[id].r);
if(pos<=mid) return ask(tol,pos);
else return ask(tor,pos);
}
long long sum[1000001];
auto lowbit=[](int x){return x&(-x);};
inline void add(int x,int val){
if(!x) return;
while(x<=m){
sum[x]+=val;
x+=lowbit(x);
}
}
inline long long ask(int x){
if(x<1) return 0;
long long ans=0;
while(x){
ans+=sum[x];
x-=lowbit(x);
}
return ans;
}
struct area{
int l,r;
};
vector<area>q[1000001];
long long ans[1000001];
struct operation{
int op,x,y,val;
}o[1000001];
int main(){
read(n,m,Q);
build(1,1,n);
bool has2=false;
for(int i=1;i<=m;++i){
read(o[i].op);
if(o[i].op==1){
read(o[i].x,o[i].y);
int hx=ask(1,o[i].x),
hy=ask(1,o[i].y);
change(1,o[i].x,o[i].x,hy);
change(1,o[i].y,o[i].y,hx);
}
else if(o[i].op==2){
has2=true;
read(o[i].x,o[i].y,o[i].val);
change(1,o[i].x,o[i].y,i);
}
else{
read(o[i].x);
o[i].y=ask(1,o[i].x);
}
}
if(!has2){
for(int i=1;i<=Q;++i){
write(0);
putchar('\n');
}
return 0;
}
for(int i=1;i<=Q;++i){
inread(l);inread(r);
q[r].push_back({l,i});
}
// cout<<"?"<<endl;
for(int i=1;i<=m;++i){
// cout<<i<<" "<<m<<endl;
if(o[i].op==3){
add(o[i].y,o[o[i].y].val);
// cout<<i<<" "<<m<<"?"<<endl;
}
for(area j:q[i]){
ans[j.r]=ask(i)-ask(j.l-1);
}
}
// cout<<"?"<<endl;
for(int i=1;i<=Q;++i){
write(ans[i]);
putchar('\n');
}
}