BZOJ 3872 Ant colony
Description
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 \(g\) groups of \(m_{1},m_{2},...,m_{g}\) 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 \(d\) outgoing corridors yet unexplored by the group, the group divides into \(d\) groups of equal size. Each newly created group follows one of the \(d\) corridors. If \(d=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 \(d\).
The following figure depicts \(m\) ants upon entering a chamber with three outgoing unexplored corridors, dividing themselves into three (equal) groups of (\(\lfloor \frac{m}{3} \rfloor\)) 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 \(k\) ants. We want to know how many ants the anteater will eat.
给定一棵有\(n\)个节点的树。在每个叶子节点,有\(g\)群蚂蚁要从外面进来,其中第\(i\)群有\(m[i]\)只蚂蚁。这些蚂蚁会相继进入树中,而且要保证每一时刻每个节点最多只有一群蚂蚁。这些蚂蚁会按以下方式前进:
·在即将离开某个度数为\(d+1\)的点时,该群蚂蚁有\(d\)个方向还没有走过,这群蚂蚁就会分裂成\(d\)群,每群数量都相等。如果\(d=0\),那么蚂蚁会离开这棵树。
·如果蚂蚁不能等分,那么蚂蚁之间会互相吞噬,直到可以等分为止,即一群蚂蚁有\(m\)只,要分成\(d\)组,每组将会有\(\lfloor \frac{m}{d} \rfloor\)只,如下图。
一只饥饿的食蚁兽埋伏在一条边上,如果有一群蚂蚁通过这条边,并且数量恰为\(k\)只,它就会吞掉这群蚂蚁。请计算一共有多少只蚂蚁会被吞掉。
Input
The first line of the standard input contains three integers \(n, g, k\) (\(2 \le n,g \le 1000000,1 \le k \le 10^9\)), 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 \(n\).
The second line contains g integers \(m[1],m[2],...,m[g]\)(\(1 \le m[i] \le 10^{9}\)), separated by single spaces, where \(m[i]\) gives the number of ants in the \(i\)-th group at every entrance to the ant hill. The \(n-1\) lines that follow describe the corridors within the ant hill; the \(i\)-th such line contains two integers \(a[i],b[i] (1 \le a[i],b[i] \le n)\), separated by a single space, that indicate that the chambers no.\(a[i]\) and \(b[i]\) are linked by a corridor. The anteater has dug into the corridor that appears first on input.
第一行包含三个整数\(n,g,k\),表示点数、蚂蚁群数以及\(k\)。
第二行包含\(g\)个整数\(m[1],m[2],...,m[g]\),表示每群蚂蚁中蚂蚁的数量。
接下来\(n-1\)行每行两个整数,表示一条边,食蚁兽埋伏在输入的第一条边上。
Output
Your program should print to the standard output a single line containing a single integer: the number of ants eaten by the anteater.
一个整数,即食蚁兽能吃掉的蚂蚁的数量。
Sample Input
7 5 3
3 4 1 9 11
1 2
1 4
4 3
4 5
4 6
6 7
Sample Output
21
一道比较水的题目吧(我都想得出,还不水吗???)
从食蚁兽潜伏的那两条线段开始bfs,初步推出食蚁兽能吃掉\(k\)只蚂蚁,需要有几只蚂蚁到达该点(该范围一定是一段连续的区间)。递推到叶子之后统计答案即可(lower_bound区间左端点\(l\),upper_bound区间右端点\(r\),\(r-l\)即为对于答案的贡献)。
#include<queue>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
using namespace std;
typedef long long ll;
#define maxn (1000010)
int side[maxn],toit[maxn*2],next[maxn*2],m[maxn],g,k,n;
int S,T,cnt,d[maxn];
bool in[maxn]; ll arr[maxn][2],ans,inf;
inline int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
inline void add(int a,int b)
{
next[++cnt] = side[a]; side[a] = cnt;
toit[cnt] = b; ++d[a];
}
inline void ins(int a,int b) { add(a,b); add(b,a); }
inline void bfs()
{
queue <int> team;
arr[S][0] = (ll)k*(ll)max((d[S]-1),1); arr[S][1] = min((ll)(k+1)*(ll)max((d[S]-1),1)-1LL,inf);
arr[T][0] = (ll)k*(ll)max((d[T]-1),1); arr[T][1] = min((ll)(k+1)*(ll)max((d[T]-1),1)-1LL,inf);
in[S] = in[T] = true;
if (arr[S][0] < inf) team.push(S); if (arr[T][0] < inf) team.push(T);
while (!team.empty())
{
int now = team.front(); team.pop();
if (d[now] == 1)
{
int l = lower_bound(m+1,m+g+1,arr[now][0])-m,r = upper_bound(m+1,m+g+1,arr[now][1])-m;
ans += (ll)(r-l)*(ll)k;
}
for (int i = side[now];i;i = next[i])
{
int v = toit[i];
if (in[toit[i]]) continue;
in[toit[i]] = true;
arr[v][0] = arr[now][0]*(ll)(max(d[v]-1,1));
arr[v][1] = min((arr[now][1]+1LL)*(ll)(max(d[v]-1,1))-1LL,inf);
if (arr[v][0] < inf) team.push(v);
}
}
}
int main()
{
freopen("3872.in","r",stdin);
freopen("3872.out","w",stdout);
n = read(); g = read(); k = read();
for (int i = 1;i <= g;++i) m[i] = read();
sort(m+1,m+g+1); inf = m[g]+1;
S = read(); T = read();
ins(S,T);
for (int i = 2;i < n;++i) ins(read(),read());
bfs(); printf("%lld",ans);
fclose(stdin); fclose(stdout);
return 0;
}
BZOJ 3872 Ant colony的更多相关文章
- bzoj 3872: [Poi2014]Ant colony -- 树形dp+二分
3872: [Poi2014]Ant colony Time Limit: 30 Sec Memory Limit: 128 MB Description There is an entranc ...
- [BZOJ3872][Poi2014]Ant colony
[BZOJ3872][Poi2014]Ant colony 试题描述 There is an entrance to the ant hill in every chamber with only o ...
- Codeforces 474 F. Ant colony
线段树求某一段的GCD..... F. Ant colony time limit per test 1 second memory limit per test 256 megabytes inpu ...
- Codeforces Round #271 (Div. 2) F. Ant colony 线段树
F. Ant colony time limit per test 1 second memory limit per test 256 megabytes input standard input ...
- 【BZOJ3872】Ant colony(二分,动态规划)
[BZOJ3872]Ant colony(二分,动态规划) 题面 又是权限题... Description There is an entrance to the ant hill in every ...
- 【BZOJ3872】[Poi2014]Ant colony 树形DP+二分
[BZOJ3872][Poi2014]Ant colony Description 给定一棵有n个节点的树.在每个叶子节点,有g群蚂蚁要从外面进来,其中第i群有m[i]只蚂蚁.这些蚂蚁会相继进入树中, ...
- CodeForces 474F Ant colony ST+二分
Ant colony 题解: 因为一个数是合法数,那么询问区间内的其他数都要是这个数的倍数,也就是这个区间内的gcd刚好是这个数. 对于这个区间的gcd来说,不能通过前后缀来算. 所以通过ST表来询问 ...
- Codeforces G. Ant colony
题目描述: F. Ant colonytime limit per test1 secondmemory limit per test256 megabytesinputstandard inputo ...
- bzoj 3872 [Poi2014]Ant colony——二分答案
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3872 可以倒推出每个叶子节点可以接受的值域.然后每个叶子二分有多少个区间符合即可. 注意一开 ...
随机推荐
- Some code changes cannot be hot swapped into a running virtual machine,
翻译一下:不能热交换到运行虚拟机,一些代码变化不能热交换到运行虚拟机,如更改名称或介绍的方法错误运行代码.解决方法:增加.删除类文件或者在一个类中增加.删除方法时,是不能够热部署到服务上的.这时候需要 ...
- Log4J2基本配置
[1]. Log4J2入门: <1>. 导入Jar包: log4j-api-2.0-beta9.jar log4j-core-2.0-beta9.jar <2>. 编写代码: ...
- Android(java)学习笔记202:Handler消息机制的原理和实现
联合学习 Android 异步消息处理机制 让你深入理解 Looper.Handler.Message三者关系 1. 首先我们通过一个实例案例来引出一个异常: (1)布局文件activity_m ...
- linux系统时间同步更新
linux中的时钟分为系统时钟和硬件时钟.硬件时钟通过bioss进行设置,系统时钟通过linux kernel进行设置.当linux启动时,bioss会读取系统时钟的时钟设置.在linux中涉及时间设 ...
- js获取图片高度
js获取图片高度时经常会获取的图片高度为0,原因是图片未加载完毕.第一次加载时,显示0(火狐等部分浏览器显示24).待加载完毕后,再刷新,显示图片高度258. var oImg = document. ...
- c#正则表达式采集数据
protected void Page_Load(object sender, EventArgs e){ StringBuilder MyStringBuilder = new StringBuil ...
- 与数据库打交道的Adapter----SimpleCursorAdapter
http://www.cnblogs.com/wenjiang/p/3196486.html 程序员是这个世界上最神奇的职业,因为几乎所有其他职业的人都能转到该行来,只要他智力正常,有接受过正规的编程 ...
- bootstrap 下的 validation插件
http://reactiveraven.github.io/jqBootstrapValidation/
- Java封装的与当前时间比较,得到多少年,多少月,多少天前,多少小时前,多小分钟前
public class CalendarCal { /** * 与当前时间比较,得到多少年,多少月,多少天前,多少小时前,多小分钟前 * * @param calendar * ...
- UITableViewCell 添加长按手势
UILongPressGestureRecognizer *lpgr = [[UILongPressGestureRecognizer alloc] initWithTarget:self actio ...