[BZOJ 3681]Ariettad
终于是打完了 CH 上的数据结构专场了……
不过看样子还有一套 5555
传送门: http://ch.ezoj.tk/contest/CH%20Round%20%2351%20-%20Shinrein%E7%A5%AD%20%231
直接连边会 TLE。但是,连出去的点都是一个连续点集……
——于是我们机智的用了 “可合并线段树”(好,) 以空间换取时间
在网络流的图中塞一棵可持久化线段树进去,这样区间连边时就只会连 O(logn) 的边了
智商++,而且这道题让我重新认识了线段树的合并,真是受益匪浅~
不过——
我恨卡内存!!!
#include <cstdio>
#include <cstring>
#include <vector>
const int inf=0x7FFFFFFF;
const int sizeOfSegmemt=;
const int sizeOfEdge=;
const int sizeOfPoint=;
const int sizeOfNote=; int n, m;
int P[sizeOfNote], H[sizeOfNote];
std::vector<int> son[sizeOfPoint];
inline int min(int, int);
inline int getint();
inline void putint(int); struct seg {int p; seg * l, * r;};
seg memory_seg[sizeOfSegmemt], * port_seg=memory_seg;
seg * t[sizeOfPoint];
inline seg * newseg();
void insert(seg *& , int, int, int, int);
seg * merge(seg * , seg * );
void query(seg * , int, int, int, int, int);
void dfs(int); int S, T;
int V, E;
struct edge {int point, flow; edge * next, * pair;};
edge memory_edge[sizeOfEdge], * port_edge=memory_edge;
edge * e[sizeOfPoint];
int h[sizeOfPoint], gap[sizeOfPoint];
inline edge * newedge(int, int, edge * );
inline void link(int, int, int);
int head, tail, queue[sizeOfPoint];
inline bool bfs();
edge * r[sizeOfPoint], * p[sizeOfPoint];
int aug[sizeOfPoint];
inline int isap(); int main()
{
int D, L, R, M; n=getint(), m=getint();
S=; T=n+; V=n+;
for (int i=;i<=n;i++)
{
P[i]=getint();
son[P[i]].push_back(i);
}
for (int i=;i<=n;i++)
H[i]=getint(); dfs(); for (int i=;i<=m;i++)
{
L=getint(), R=getint(), D=getint(), M=getint();
link(S, ++V, M);
query(t[D], , n, L, R, V);
} putint(isap()); return ;
}
inline int min(int x, int y)
{
return x<y?x:y;
}
inline int getint()
{
register int num=;
register char ch;
do ch=getchar(); while (ch<'' || ch>'');
do num=num*+ch-'', ch=getchar(); while (ch>='' && ch<='');
return num;
}
inline void putint(int num)
{
char stack[];
register int top=;
if (num==) stack[top=]='';
for ( ;num;num/=) stack[++top]=num%+'';
for ( ;top;top--) putchar(stack[top]);
putchar('\n');
}
inline edge * newedge(int point, int flow, edge * next)
{
edge * ret=port_edge++;
ret->point=point; ret->flow=flow; ret->next=next;
return ret;
}
inline void link(int u, int v, int f)
{
e[u]=newedge(v, f, e[u]); e[v]=newedge(u, , e[v]);
e[u]->pair=e[v]; e[v]->pair=e[u];
}
inline bool bfs()
{
memset(h, 0xFF, sizeof(h)); h[T]=;
head=tail=;
for (queue[tail++]=T;head<tail;head++)
{
int u=queue[head];
gap[h[u]]++;
for (edge * i=e[u];i;i=i->next) if (h[i->point]==-)
{
h[i->point]=h[u]+;
queue[tail++]=i->point;
}
}
return h[S]>-;
}
inline int isap()
{
int flow=;
int hmin=; if (!bfs()) return ; memcpy(r, e, sizeof(e));
memset(p, , sizeof(p));
aug[S]=inf;
for (int u=S;h[S]<V; )
{
if (u==T)
{
flow+=aug[T];
for (edge * i=p[T];i;i=p[i->point])
i->pair->flow-=aug[T], i->flow+=aug[T], aug[i->point]-=aug[T];
for (edge * i=p[T];i;i=p[i->point]) if (aug[i->point])
{
u=i->point;
break;
}
} edge *& i=r[u];
for ( ;i && (!i->flow || h[i->point]+!=h[u]);i=i->next);
if (i)
{
p[i->point]=i->pair; aug[i->point]=min(aug[u], i->flow);
u=i->point;
}
else
{
if (!--gap[h[u]]) break;
hmin=V;
for (edge * j=e[u];j;j=j->next) if (j->flow && h[j->point]+<hmin)
{
r[u]=j;
hmin=h[j->point]+;
}
++gap[h[u]=hmin];
u=u==S?S:p[u]->point;
}
} return flow;
} inline seg * newseg()
{
seg * ret=port_seg++;
ret->p=++V; ret->l=ret->r=NULL;
return ret;
}
void insert(seg *& t, int l, int r, int p, int v)
{
if (!t) t=newseg();
if (l==r) link(t->p, v, inf);
else
{
int m=(l+r)>>;
if (p<=m) insert(t->l, l, m, p, v), link(t->p, t->l->p, inf);
else insert(t->r, m+, r, p, v), link(t->p, t->r->p, inf);
}
}
seg * merge(seg * x, seg * y)
{
if (!x) return y;
if (!y) return x;
seg * t=newseg();
bool leaf=true;
if ((t->l=merge(x->l, y->l))) link(t->p, t->l->p, inf), leaf=false;
if ((t->r=merge(x->r, y->r))) link(t->p, t->r->p, inf), leaf=false;
if (leaf) link(t->p, x->p, inf), link(t->p, y->p, inf);
return t;
}
void query(seg * t, int l, int r, int ql, int qr, int p)
{
if (!t) return ;
if (l==ql && r==qr) link(p, t->p, inf);
else
{
int m=(l+r)>>;
if (qr<=m) query(t->l, l, m, ql, qr, p);
else if (ql>=m+) query(t->r, m+, r, ql, qr, p);
else query(t->l, l, m, ql, m, p), query(t->r, m+, r, m+, qr, p);
}
}
void dfs(int u)
{
link(u, T, );
insert(t[u], , n, H[u], u);
for (unsigned int i=;i<son[u].size();i++)
{
dfs(son[u][i]);
t[u]=merge(t[u], t[son[u][i]]);
}
}
我恨卡内存++
[BZOJ 3681]Ariettad的更多相关文章
- bzoj 3681 Arietta
一棵有根树,每个点有一个音高,有 $m$ 中弹奏方法,每种方法可以弹奏 $d$ 子树中音高在 $[l,r]$ 间的音符,每种方法最多弹 $t$ 次 求最多能弹出多少个音符 $n \leq 10000$ ...
- BZOJ 3681 线段树合并+网络流
思路: 暴力建图有n*m条边 考虑怎么优化 (那就只能加个线段树了呗) 然后我就不会写了..... 抄了一波题解 //By SiriusRen #include <bits/stdc++.h&g ...
- 【BZOJ 3681】Arietta
传送门 题目描述 Arietta 的命运与她的妹妹不同,在她的妹妹已经走进学院的时候,她仍然留在山村中. 但是她从未停止过和恋人 Velding 的书信往来.一天,她准备去探访他. 对着窗外的阳光,临 ...
- BZOJ 2127: happiness [最小割]
2127: happiness Time Limit: 51 Sec Memory Limit: 259 MBSubmit: 1815 Solved: 878[Submit][Status][Di ...
- BZOJ 3275: Number
3275: Number Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 874 Solved: 371[Submit][Status][Discus ...
- BZOJ 2879: [Noi2012]美食节
2879: [Noi2012]美食节 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 1834 Solved: 969[Submit][Status] ...
- bzoj 4610 Ceiling Functi
bzoj 4610 Ceiling Functi Description bzoj上的描述有问题 给出\(n\)个长度为\(k\)的数列,将每个数列构成一个二叉搜索树,问有多少颗形态不同的树. Inp ...
- BZOJ 题目整理
bzoj 500题纪念 总结一发题目吧,挑几道题整理一下,(方便拖板子) 1039:每条线段与前一条线段之间的长度的比例和夹角不会因平移.旋转.放缩而改变,所以将每条轨迹改为比例和夹角的序列,复制一份 ...
- 【sdoi2013】森林 BZOJ 3123
Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整数N,M,T,分别表示节点数.初始边数.操作数.第三行包含N个非负整数 ...
随机推荐
- sac 文档使用
目前我遇到的问题是我想要得到BHE,BHN 方向的数据,但发现IRIS下载的数据都是BH1,BH2 方向的,很困惑,请教大神后发现,原来IRIS之所以提供BH1,BH2方向是因为很多时候台站的水平方向 ...
- 利用calc计算宽度
width:calc(100% - 40px)可用 + - * / 进行计算(ie9+) 注:计算符号前后必须跟上空格.
- Remove openjdk in Ubuntu/Configure jdk and running adb in 64-bit Ubuntu
sudo apt-get autoremove openjdk-7-jre sudo apt-get purge openjdk* java -version No openjdk available ...
- The Pragmatic Programmer Quick Reference Guide
1.关心你的技艺 Care About Your Craft 如果不在乎能否漂亮地开发出软件,你又为何要耗费生命去开发软件呢? 2.思考!你的工作 Think! About Your Work 关掉自 ...
- Permutation
(M) Permutations (M) Permutations II (M) Permutation Sequence (M) Palindrome Permutation II
- 微信公众账号 Senparc.Weixin.MP SDK 开发教程 索引
Senparc.Weixin.MP SDK从一开始就坚持开源的状态,这个过程中得到了许多朋友的认可和支持. 目前SDK已经达到比较稳定的版本,这个过程中我觉得有必要整理一些思路和经验,和大家一起分享. ...
- vue吃进去的object已经变了样,不在是原来的!
直接上代码: class data { public list:number[]; constructor() { this.list=[0,0,0]; } } class methods exten ...
- How can I protect derived classes from breaking when I change the internal parts of the base class?
How can I protect derived classes from breaking when I change the internal parts of the base class? ...
- 基于HTML5的网络拓扑图
电信网管系统中,设备状态信息的实时展示非常重要,通常会挂载一堆图标来展示状态或告警信息,图标的信息量有限,有时需要更详细的面板,甚至以图表的形式展现,本文将结合最近客户提到的需求,使用 Qunee1. ...
- WPF与WinForm开发有什么区别?
转自http://hi.baidu.com/leoliu83/blog/item/1d1a4a66dcb41134aa184cfd.html WPF开发于WinForm之后,从技术发展的角度,WPF比 ...