解析

二分的單調性較為明顯，一路推導下去的性質都不算太難想，正解的思路還是不難想到的
但從頭到尾都實現很考驗思維的嚴密性和代碼能力
然后我就雙重被考驗掛了qwq
第一交的時候一個地方把倍增的dis數組寫成pl，判斷封口也有問題…
但竟然有70
真實考試的時候絕對不能出打錯數組這種錯！提交前還是要謹慎一些！
另外那個判封口…只能說路漫漫其修遠兮了…

我對于最后根下方兒子（暫且叫它次根好了）與軍隊貪心匹配的地方寫的和題解不太一樣

題解的實現是：爬到根下方的軍隊不考慮，然后把所有軍隊提出來一起貪心

我的做法是先給有軍隊的次根一個剩余最少的軍隊，再維護一個堆看是否需要用外面的軍隊把里面的軍隊替換出來

這兩個實現我覺得差不多吧，談不上孰優孰劣

然而我就寫掛了qwq

代碼

#include<bits/stdc++.h> using namespace std; #define ll long long const int N=1e5+100; ll read() {ll x=0,f=1;char c=getchar();while(!isdigit(c)) {if(c=='-')f=-1;c=getchar();}while(isdigit(c)) {x=x*10+(c^48);c=getchar();}return x*f; }int n,m; int x[N]; struct node{int to,nxt,w; }p[N<<1]; int fi[N],ecnt; inline void addline(int x,int y,int w){p[++ecnt]=(node){y,fi[x],w};fi[x]=ecnt;return; } ll dis[N][20]; int pl[N][20]; int que[N],num,jd[N]; void dfs(int x,int f){for(int k=1;k<=16;k++){pl[x][k]=pl[pl[x][k-1]][k-1];dis[x][k]=dis[x][k-1]+dis[pl[x][k-1]][k-1];}jd[x]=1;for(int i=fi[x];~i;i=p[i].nxt){int to=p[i].to;if(to==f) continue;jd[x]=0;if(x==1) que[++num]=to;pl[to][0]=x;dis[to][0]=p[i].w;dfs(to,x);}//printf("x=%d jd=%d\n",x,jd[x]);return; }bool vis[N],tag[N]; vector<ll>v[N]; int solve(int x,int f){//printf("solve:x=%d vis=%d jd=%d\n",x,vis[x],jd[x]);if(vis[x]) return tag[x]=1;if(jd[x]) return tag[x]=0;tag[x]=1;for(int i=fi[x];~i;i=p[i].nxt){int to=p[i].to;if(to==f) continue;tag[x]&=solve(to,x);}//printf("x=%d tag=%d\n",x,tag[x]);return tag[x]; } ll w[N],cnt,add[N],o; priority_queue<ll>q; struct aaa{ll ned,val;bool operator < (const aaa oth)const{return ned<oth.ned;} }a[N];; int numa; bool cmp(ll a,ll b){return a>b;} void init(){memset(vis,0,sizeof(vis));memset(tag,0,sizeof(tag));for(int i=1;i<=num;i++){v[que[i]].clear();v[que[i]].shrink_to_fit();}cnt=o=numa=0;while(!q.empty()) q.pop(); } bool check(ll tot){init();for(int o=1;o<=m;o++){int pos=x[o];ll lft=tot;for(int k=16;k>=0;k--){if(dis[pos][k]>lft||pl[pos][k]==0||pl[pos][k]==1) continue;lft-=dis[pos][k];pos=pl[pos][k];} //printf("pos=%d lft=%lld\n",pos,lft);if(pl[pos][0]==1){v[pos].push_back(lft);}else vis[pos]=1;}if(solve(1,0)) return true;for(int i=1;i<=num;i++){int now=que[i];if(v[now].size()==0&&!tag[now]){w[++cnt]=dis[now][0];continue;}sort(v[now].begin(),v[now].end());if(!tag[now]){if(v[now][0]-dis[now][0]>0){//printf("now=%d v=%lld dis=%lld\n",now,v[now][0],dis[now][0]);a[++numa]=(aaa){dis[now][0],v[now][0]-dis[now][0]};}for(int j=1,tp=v[now].size();j<tp;j++){if(v[now][j]>dis[now][0]) add[++o]=v[now][j]-dis[now][0];}}else{for(int j=0,tp=v[now].size();j<tp;j++){if(v[now][j]>dis[now][0]) add[++o]=v[now][j]-dis[now][0];}}}if(o<cnt) return false;sort(a+1,a+1+numa);sort(add+1,add+1+o);sort(w+1,w+1+cnt,cmp);//printf("w: ");for(int i=1;i<=cnt;i++) printf("%lld ",w[i]);putchar('\n');//printf("add: ");for(int i=1;i<=o;i++) printf("%lld ",add[i]);putchar('\n');//printf("aaa: ");for(int i=1;i<=numa;i++) printf("(%lld %lld) ",a[i].ned,a[i].val);putchar('\n');int pos=0;for(int k=1;k<=o;k++){while(pos<numa&&add[k]>=a[pos+1].ned){pos++;q.push(a[pos].val);}if(q.empty()) continue;ll now=q.top();q.pop();if(add[k]<now) swap(add[k],now);q.push(now);}sort(add+1,add+1+o,cmp);//printf("add: ");for(int i=1;i<=o;i++) printf("%lld ",add[i]);putchar('\n');for(int i=1;i<=cnt;i++){if(w[i]>add[i]) return false;}return true; }int main() { #ifndef ONLINE_JUDGE//freopen("a.in","r",stdin);//freopen("a.out","w",stdout); #endifmemset(fi,-1,sizeof(fi));ecnt=-1;n=read();for(int i=1;i<n;i++){int x=read(),y=read(),w=read();addline(x,y,w);addline(y,x,w);}m=read();for(int i=1;i<=m;i++) x[i]=read();dfs(1,0);//printf("check=%d\n",check(8));//return 0;ll st=0,ed=1e14;while(st<ed){ll mid=(st+ed)>>1;if(check(mid)) ed=mid;else st=mid+1;}printf("%lld\n",st<1e14?st:-1);return 0; } /* 4 1 2 10 1 3 1 4 1 5 3 3 4 4 */

總結

以上是生活随笔為你收集整理的NOIP2012:疫情控制（二分、贪心、树上倍增）的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。