集训Day2
雅礼集训2017Day2
T1
给你一个水箱,水箱里有n-1个挡板,水遵循物理定律
给你m个条件,表示第i个格子上面y+1高度的地方有或没有水
现在给你无限的水从任意地方往下倒,问最多满足多少条件
n,m 1e5
SOL:
考虑答案的表示方法,肯定是类似于dp["区间1到n"]这种的
“区间1到n”这个东西我们可以状压,就过了20%
另外10%的数据只有“有水”的条件 我们cout<<m
另外30%的n平方做法肯定是要用到dp(我没有想出来这个东西怎么搞到n平方)
。。。
这里是全知全能的XiongGod提供的一种n平方做法
对于每一个条件,可以将其视为l到r区间的最高/最低高度为x;
把所有条件记录下来后就成了区间选择问题的版题,分成上下界两种讨论dp就好了
考虑暴力碾标算
一个区间能能不能“有水”取决于中间最长的那个木板
于是可以考虑递归建一棵“线段树”
这样一个区间的问题就可以递归到左右儿子解决
且因为这是棵线段树,最多nlogn个点
然后树形DP
$d[i]$表示i节点水没有漫出去的情况
$f[i]$表示i节点水漫出去的情况
式子由于...这是要给学妹(实际只有学弟吧?)看的
所以你们自己推吧
T2
一个100*100的棋盘,有地方有障碍
Alice和Bob两个人玩游戏,Alice放一个棋子,Bob先手二人轮流移动棋子
要求:不能移到障碍上且走过的地方不能走
不能动就输了
求Alice的必胜点
SOL:
原题啊...
将棋盘黑白染色建出二分图 有障碍直接跳过就行了
考虑对于一个完美匹配,Bob(先手)按匹配边走就可以了
对于不完美的最大匹配,Alice可以把棋放在“不一定是最大匹配”的地方,让Bob走非匹配边,自己走匹配边
现在就是要考虑如何求这个“不一定是最大匹配”的点
跑完dinic,从S开始跑未满流的边,跑到左边的且非S的点就是“不一定是最大匹配”的
T开始跑未满流的边,跑到右边且非T的点就是“不一定是最大匹配的”
为什么正确可以通过交错路定理来证明一下
时间复杂度是O(您用的二分图匹配算法 + 棋盘大小)
100 + 100 + 0 = 200?
T3前30随便写可是没时间写了...
话说是不是上午时间不够啊
如果再给我半个小时就是100 + 100 + 30 = 230了(再次大众分
#include<bits/stdc++.h>
#define LL long long
inline int read() {
int x = ,f = ;
char ch = getchar();
for(; !isdigit(ch); ch = getchar())if(ch == '-')f = -f;
for(; isdigit(ch); x = * x + ch - '',ch = getchar());
return x * f;
}
const int maxn = 2e5 + ,inf = ;
int n,m,N;
int fa[maxn],bot[maxn],top[maxn],dx[maxn];
int ST[][maxn];
int son[maxn][];
int f[maxn],d[maxn]; struct block {
int h,pos;
} hs[maxn];
bool cmp(const block &a,const block &b) {
if(a.h == b.h)return a.pos < b.pos;
return a.h < b.h;
}
struct info {
int h,type;
bool operator < (const info &a)const {
if(h == a.h)return type < a.type;
return h < a.h;
}
}; std::vector<info> vec[maxn]; inline int find(int x) {
return x == fa[x] ? x : fa[x] = find(fa[x]);
} int main() {
int T = read();
while(T--) {
n = read(),m = read();
N = n;
memset(dx,,sizeof(dx));
memset(fa,,sizeof(fa));
memset(f,,sizeof(f));
memset(d,,sizeof(d));
int S,T,sum,sz,tmp;
for(int i=; i<n; i++) {
hs[i].h = read();
hs[i].pos = i;
}
std::sort(hs + ,hs + n,cmp);
bot[] = inf;
for(int i=; i<=n; i++)top[i] = i,fa[i] = i;
for(int i=; i<n; i++) {
int fx = find(hs[i].pos),fy = find(hs[i].pos + );
bot[++N] = hs[i].h;
ST[][N] = ;
ST[][top[fx]] = N;
ST[][top[fy]] = N;
son[N][] = top[fx],son[N][] = top[fy];
fa[fy] = fx;
top[fx] = N;
}
for(int i=; i<=N; i++)vec[i].clear();
for(int i=; i<=; i++)
for(int j=; j<=N; j++)
ST[i][j]=ST[i-][ST[i-][j]];
int x,y,k;
for(int i=; i<=m; i++) {
x = read(),y = read(),k = read();
for(int j=; j>=; j--)
if(bot[ST[j][x]] <= y)x = ST[j][x];
vec[x].push_back((info) {y,k});
dx[x] += !k;
}
for(int i=; i<=N; i++)std::sort(vec[i].begin(),vec[i].end());
for(int i=; i<=N; i++) {
if(!vec[i].empty()) {
sz = vec[i].size();
S = ;
d[i] = sum = dx[i] + ( i > n ? f[son[i][]] + f[son[i][]] : );
while(S < sz) {
T = S;
tmp = (vec[i][T].type ? : -);
while(T + < sz && vec[i][T+].h == vec[i][T].h)
++T,tmp += (vec[i][T].type ? : -);
sum += tmp;
d[i] = std::max(d[i],sum);
S = T + ;
}
f[i] = sum;
}
if(i > n) {
d[i] = std::max(d[i],dx[i] + d[son[i][]] + d[son[i][]]);
f[i] = std::max(f[i],f[son[i][]] + f[son[i][]]);
}
}
printf("%d\n",d[N]);
}
}
T1(压行太严重就用了Astyle)
#include<bits/stdc++.h>
using namespace std;
const int _MAX = ,MAX_POINT = ,MAXE = ,inf = ;
const int dx[] = {,,-,,};
const int dy[] = {,,,,-};
int ans[MAX_POINT];
int m,n;
int S = ,T = MAX_POINT - ;
char s[_MAX][_MAX];
int pos[_MAX][_MAX],col[MAX_POINT],vis[MAX_POINT];
struct DINIC
{
int first[MAX_POINT],cnt;
int next[MAXE],to[MAXE],caps[MAXE];
int deep[MAX_POINT];
DINIC(){cnt = ;}
void add(int x,int y,int f)
{
next[++cnt] = first[x];
to[cnt] = y;
caps[cnt] = f;
first[x] = cnt;
}
void insert(int x,int y,int f)
{
add(x,y,f);
add(y,x,);
}
bool BFS()
{
queue<int> q;
memset(deep,,sizeof(deep));
deep[S] = ;
q.push(S);
while(!q.empty())
{
int x = q.front(); q.pop();
for(int i = first[x]; i; i = next[i])
if(caps[i] && !deep[to[i]])
{
deep[to[i]] = deep[x] + ;
q.push(to[i]);
if(to[i] == T) return true;
}
}
return false;
}
int dfs(int x,int f)
{
if(x == T) return f;
int temp = f;
for(int i = first[x]; i; i = next[i])
if(caps[i] && deep[to[i]] == deep[x] + && temp)
{
int w = dfs(to[i],min(caps[i],temp));
if(!w) deep[to[i]] = ;
caps[i] -= w;
caps[i^] += w;
temp -= w;
}
return f - temp;
}
}G; int ctt;
void DFS(int x,int f)
{
vis[x] = ;
if(col[x] == f && x != S && x != T)ans[++ctt] = x;
for(int i = G.first[x]; i; i = G.next[i])
if(G.caps[i] == f && !vis[G.to[i]])DFS(G.to[i],f);
}
void solveGraph()
{
for(int i = ; i <= m; i++)
for(int j = ; j <= n; j++)
{
if(s[i][j] == '#') continue;
if(!((i + j)&))G.insert(S,pos[i][j],),col[pos[i][j]] = ;
else
{
G.insert(pos[i][j],T,);
continue;
}
for(int k = ; k <= ; k++)
{
int fx = i + dx[k],fy = j + dy[k];
if(fx < || fy < || fx > m || fy > n) continue;
if(s[fx][fy] == '.')G.insert(pos[i][j],pos[fx][fy],);
}
}
while(G.BFS())G.dfs(S,inf);
DFS(S,);memset(vis,,sizeof(vis));DFS(T,);
}
int main()
{
scanf("%d%d",&m,&n);
for(int i = ; i <= m; i++)
scanf("%s",s[i] + );
int cnt = ;
for(int i = ; i <= m; i++)
for(int j = ; j <= n; j++)
pos[i][j] = ++cnt;
solveGraph();
printf("%d\n",ctt);
sort(ans + ,ans + ctt + );
for(int i = ; i <= ctt; i++)printf("%d %d\n",(ans[i] - ) / n + ,(ans[i] - ) % n + );
return ;
}
T2
集训Day2的更多相关文章
- WC集训DAY2笔记 组合计数 part.1
目录 WC集训DAY2笔记 组合计数 part.1 基础知识 组合恒等式 错排数 卡特兰数 斯特林数 伯努利数 贝尔数 调和级数 后记 补完了几天前写的东西 WC集训DAY2笔记 组合计数 part. ...
- 2019暑期金华集训 Day2 线性代数
自闭集训 Day2 线性代数 高斯消元 做实数时,需要找绝对值最大的作为主元,以获取更高精度. 在欧几里得环(简单例子是模合数)意义下也是对的.比如模合数意义下可以使用辗转相除法消元. 欧几里得环:对 ...
- 暑假集训Day2 互不侵犯(状压dp)
这又是个状压dp (大型自闭现场) 题目大意: 在N*N的棋盘里面放K个国王,使他们互不攻击,共有多少种摆放方案.国王能攻击到它上下左右,以及左上左下右上右下八个方向上附近的各一个格子,共8个格子. ...
- 2022寒假集训day2
day1:学习seach和回溯,初步了解. day2:深度优化搜索 T1 洛谷P157:https://www.luogu.com.cn/problem/P1157 题目描述 排列与组合是常用的数学方 ...
- 考前停课集训 Day2 非
因为太长了 所以一天一天分开发 Day2 昨天晚上没开黑车 没脱衣服就睡了 可能是我难受了…… 新的一天. 早上好. 我没去晨跑,早上先和团长集合了,没看见rkbudlo来 于是就先吃饭了 去机房的时 ...
- 长沙雅礼中学集训-------------------day2
怎么说呢,今天的题特别的神奇,因为emmmmmm,T1看不懂(一直到现在还没有理解明白期望概率什么的),T2题面看不懂+扩展欧几里得求逆元怎么求我忘了,T3哇,终于看懂一题了,然而写了个50分的程序但 ...
- HZNU-ACM寒假集训Day2小结 二分答案
Day2 ---二分 这里直接给出模板 两种对应不同的情况 可以借助数轴理解 int bsearch_1(int l, int r) { while (l < r) { ; if (check( ...
- [日常训练]常州集训day2
T1 Description 给定$N$个点,问这$N$个点能构成的三角形个数. Input 第一行一个整数$N$,代表点数. 接下来$N$行,每行两个非负整数$X,Y$,表示一个点的坐标. Outp ...
- FJ省队集训DAY2 T2
思路:我们可以考虑三角剖分,这样问题就变成考虑三角形的选取概率和三角形内有多少个点了. 先用树状数组预处理出三角剖分的三角形中有多少个点,然后用线段树维护,先用原点极角排序,然后枚举i,再以i极角排序 ...
随机推荐
- Lua学习三----------Lua数据类型
© 版权声明:本文为博主原创文章,转载请注明出处 Lua数据类型 - Lua是动态类型语言,不需要为变量定义类型,只需要为变量赋值 - Lua有8中基本数据类型:nil.boolean.number. ...
- HDU 2473 Junk-Mail Filter 删点并查集
题目来源:pid=2473">HDU 2473 Junk-Mail Filter 题意:2中操作 M x, y 将x,y 合并到一个集合 S x 将x从所在的集合去掉 自己成为一个集合 ...
- Fragment小结
Fragment是Android3.0之后增加的新特性,通常人们叫它碎片.可是,我认为把它理解成一个View模块比較好,尽管它不是继承自View.假设阅读过源代码就知道它是内置View对象从而实现Vi ...
- js错误: Unexpected number in JSON at position 2792 value里面有双引号怎么解决
源头 出现这个报错提示,大家从错误就可以看的出来,这就是json的错误,一般来说都是json格式出现了错误,本人遇到比较多的情况就是json字符串里面出现了一些会影响json格式的符号,这次出现这个 ...
- HDFS源码分析数据块复制之PendingReplicationBlocks
PendingReplicationBlocks实现了所有正在复制的数据块的记账工作.它实现以下三个主要功能: 1.记录此时正在复制的块: 2.一种对复制请求进行跟踪的粗粒度计时器: 3.一个定期识别 ...
- 源码维护基本命令diff_patch_quilt
源码维护基本命令 一. diff--生成补丁 diff [命令行选项] 源文件 新文件 -r 递归处理相应目录 -N 包含新文件到patch -u 输出统一格式(unified format),这种格 ...
- LINUX线程初探
LINUX程序设计最重要的当然是进程与线程.本文主要以uart程序结合键盘输入控制uart的传输. 硬件平台:树莓派B+ 软件平台:raspberry 须要工具:USB转TTL(PL2303)+ ...
- vim与sublime
vim与sublime 对程序员来说,写代码是再熟悉不过的事情了,windows系统自带有记事本软件,能写写小规模的代码,可是代码量大了,它的局限性就暴露得很明显了:没有语法高亮,没有自动提示,不支持 ...
- sublime3 支持 jsx 语法
添加几个插件即可在js中快速写html babel 可以识别React,并高亮显示ES6 command+shift+p -> install package -> babel 使用 在打 ...
- nginx教程1:location 匹配规则
worker_process # 表示工作进程的数量,一般设置为cpu的核数 worker_connections # 表示每个工作进程的最大连接数 server{} # 块定义了虚拟主机 liste ...