洛谷—— P3576 [POI2014]MRO-Ant colony

https://www.luogu.org/problem/show?pid=3576

题目描述

The ants are scavenging an abandoned ant hill in search of food.

The ant hill has nn chambers and n-1n−1 corridors connecting them.

We know that each chamber can be reached via a unique path from every other chamber.

In other words, the chambers and the corridors form a tree.

There is an entrance to the ant hill in every chamber with only one corridor leading into (or out of) it.

At each entry, there are gg groups of m_1,m_2,\cdots,m_gm1,m2,⋯,mg ants respectively.

These groups will enter the ant hill one after another, each successive group entering once there are no ants inside.

Inside the hill, the ants explore it in the following way:

Upon entering a chamber with dd outgoing corridors yet unexplored by the group,the group divides into dd groups of equal size. Each newly created group follows one of the d corridors.If d=0d=0, then the group exits the ant hill.
If the ants cannot divide into equal groups, then the stronger ants eat the weaker until a perfect division is possible.Note that such a division is always possible since eventually the number of ants drops down to zero.Nothing can stop the ants from allowing divisibility - in particular, an ant can eat itself, and the last one remaining will do so if the group is smaller than dd.

The following figure depicts mm ants upon entering a chamber with three outgoing unexplored corridors, dividing themselves into three (equal) groups of \left \lfloor m/3 \right \rfloor⌊m/3⌋ ants each.

A hungry anteater dug into one of the corridors and can now eat all the ants passing through it.

However, just like the ants, the anteater is very picky when it comes to numbers.

It will devour a passing group if and only if it consists of exactly kk ants.

We want to know how many ants the anteater will eat.

给一棵树，对于每个叶子节点，都有g群蚂蚁要从外面进来，每群蚂蚁在行进过程中只要碰到岔路，就将平均地分成岔路口数-1那么多份，然后平均地走向剩下的那些岔路口，余下的蚂蚁自动消失，树上有一个关键边，假如有一群蚂蚁通过了这条边且数量恰好为k，这k只蚂蚁就被吃掉，问一共有多少只蚂蚁被吃掉

输入输出格式

输入格式：

The first line of the standard input contains three integers nn, gg, kk(2\le n,g\le 1\ 000\ 0002≤n,g≤1 000 000, 1\le k\le 10^91≤k≤109), separated by single spaces.

These specify the number of chambers, the number of ant groups and the number of ants the anteater devours at once. The chambers are numbered from 1 to nn.

The second line contains gg integers m_1,m_2,\cdots,m_gm1,m2,⋯,mg (1\le m_i\le 10^91≤mi≤109), separated by single spaces, where m_imi gives the number of ants in the ii-th group at every entrance to the ant hill. The n-1n−1 lines that follow describe the corridors within the ant hill;the ii-th such line contains two integers a_iai,b_ibi (1\le a_i,b_i\le n1≤ai,bi≤n), separated by a single space, that indicate that the chambers no. a_iai and b_ibi are linked by a corridor. The anteater has dug into the corridor that appears first on input.

输出格式：

Your program should print to the standard output a single line containing a single integer: the number of ants eaten by the anteater.

输入输出样例

输入样例#1：

输出样例#1：

说明

从关键路的端点DFS，统计出到达每个点的最大和最小的蚂蚁数，（只有当最小值不比最大的蚁群数时，才继续搜下一层）

二分统计出每个点能得到的最大蚁群数、

 #include <algorithm>

 #include <cstdio>

 #define LL long long

 const int N();

 inline void read(LL &x)

 {

     x=; register char ch=getchar();

     for(; ch>''||ch<''; ) ch=getchar();

     for(; ch>=''&&ch<=''; ch=getchar()) x=x*+ch-'';

 }

 LL n,g,k,s1,s2,gi[N];

 int head[N],sumedge;

 struct Edge {

     int v,next;

     Edge(int v=,int next=):v(v),next(next){}

 }edge[N<<];

 inline void ins(int u,int v)

 {

     edge[++sumedge]=Edge(v,head[u]);

     head[u]=sumedge;

     edge[++sumedge]=Edge(u,head[v]);

     head[v]=sumedge;

 }

 #define min(a,b) (a<b?a:b)

 #define max(a,b) (a>b?a:b)

 int du[N],dad[N],minn[N],maxx[N];

 void DFS(int u)

 {

     for(int v,i=head[u]; i; i=edge[i].next)

     {

         v=edge[i].v;

         if(dad[u]==v) continue;

         dad[v]=u; du[u]++;

     }

     for(int v,i=head[u]; i; i=edge[i].next)

     {

         v=edge[i].v;

         if(dad[u]==v) continue;

         minn[v]=minn[u]*du[u];

         maxx[v]=(maxx[u]+)*du[u]-;

         maxx[v]=min(maxx[v],gi[g]);

         if(minn[v]<=gi[g]) DFS(v);

     }

 }

 LL l,r,mid,ans;

 LL check(LL x)

 {

     LL ret=;

     for(l=,r=g; l<=r; )

     {

         mid=l+r>>;

         if(gi[mid]<x)

         {

             ret=mid;

             l=mid+;

         }

         else r=mid-;

     }

     return ret;

 }

 int Presist()

 {

     read(n),read(g),read(k);

     for(int i=; i<=g; ++i) read(gi[i]);

     read(s1);read(s2);

     for(LL u,v,i=; i<n; ++i)

         read(u),read(v),ins(u,v);

     std::sort(gi+,gi+g+);

     maxx[s1]=maxx[s2]=minn[s1]=minn[s2]=k;

     DFS(s1); DFS(s2);

     for(int i=; i<=n; ++i)

         if(!du[i]) ans+=check(maxx[i]+)-check(minn[i]);

     printf("%lld\n",ans*k);

     return ;

 }

 int Aptal=Presist();

 int main(){;}