雅礼集训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的更多相关文章

  1. WC集训DAY2笔记 组合计数 part.1

    目录 WC集训DAY2笔记 组合计数 part.1 基础知识 组合恒等式 错排数 卡特兰数 斯特林数 伯努利数 贝尔数 调和级数 后记 补完了几天前写的东西 WC集训DAY2笔记 组合计数 part. ...

  2. 2019暑期金华集训 Day2 线性代数

    自闭集训 Day2 线性代数 高斯消元 做实数时,需要找绝对值最大的作为主元,以获取更高精度. 在欧几里得环(简单例子是模合数)意义下也是对的.比如模合数意义下可以使用辗转相除法消元. 欧几里得环:对 ...

  3. 暑假集训Day2 互不侵犯(状压dp)

    这又是个状压dp (大型自闭现场) 题目大意: 在N*N的棋盘里面放K个国王,使他们互不攻击,共有多少种摆放方案.国王能攻击到它上下左右,以及左上左下右上右下八个方向上附近的各一个格子,共8个格子. ...

  4. 2022寒假集训day2

    day1:学习seach和回溯,初步了解. day2:深度优化搜索 T1 洛谷P157:https://www.luogu.com.cn/problem/P1157 题目描述 排列与组合是常用的数学方 ...

  5. 考前停课集训 Day2 非

    因为太长了 所以一天一天分开发 Day2 昨天晚上没开黑车 没脱衣服就睡了 可能是我难受了…… 新的一天. 早上好. 我没去晨跑,早上先和团长集合了,没看见rkbudlo来 于是就先吃饭了 去机房的时 ...

  6. 长沙雅礼中学集训-------------------day2

    怎么说呢,今天的题特别的神奇,因为emmmmmm,T1看不懂(一直到现在还没有理解明白期望概率什么的),T2题面看不懂+扩展欧几里得求逆元怎么求我忘了,T3哇,终于看懂一题了,然而写了个50分的程序但 ...

  7. HZNU-ACM寒假集训Day2小结 二分答案

    Day2 ---二分 这里直接给出模板 两种对应不同的情况 可以借助数轴理解 int bsearch_1(int l, int r) { while (l < r) { ; if (check( ...

  8. [日常训练]常州集训day2

    T1 Description 给定$N$个点,问这$N$个点能构成的三角形个数. Input 第一行一个整数$N$,代表点数. 接下来$N$行,每行两个非负整数$X,Y$,表示一个点的坐标. Outp ...

  9. FJ省队集训DAY2 T2

    思路:我们可以考虑三角剖分,这样问题就变成考虑三角形的选取概率和三角形内有多少个点了. 先用树状数组预处理出三角剖分的三角形中有多少个点,然后用线段树维护,先用原点极角排序,然后枚举i,再以i极角排序 ...

随机推荐

  1. Lua学习三----------Lua数据类型

    © 版权声明:本文为博主原创文章,转载请注明出处 Lua数据类型 - Lua是动态类型语言,不需要为变量定义类型,只需要为变量赋值 - Lua有8中基本数据类型:nil.boolean.number. ...

  2. HDU 2473 Junk-Mail Filter 删点并查集

    题目来源:pid=2473">HDU 2473 Junk-Mail Filter 题意:2中操作 M x, y 将x,y 合并到一个集合 S x 将x从所在的集合去掉 自己成为一个集合 ...

  3. Fragment小结

    Fragment是Android3.0之后增加的新特性,通常人们叫它碎片.可是,我认为把它理解成一个View模块比較好,尽管它不是继承自View.假设阅读过源代码就知道它是内置View对象从而实现Vi ...

  4. js错误: Unexpected number in JSON at position 2792 value里面有双引号怎么解决

    源头  出现这个报错提示,大家从错误就可以看的出来,这就是json的错误,一般来说都是json格式出现了错误,本人遇到比较多的情况就是json字符串里面出现了一些会影响json格式的符号,这次出现这个 ...

  5. HDFS源码分析数据块复制之PendingReplicationBlocks

    PendingReplicationBlocks实现了所有正在复制的数据块的记账工作.它实现以下三个主要功能: 1.记录此时正在复制的块: 2.一种对复制请求进行跟踪的粗粒度计时器: 3.一个定期识别 ...

  6. 源码维护基本命令diff_patch_quilt

    源码维护基本命令 一. diff--生成补丁 diff [命令行选项] 源文件 新文件 -r 递归处理相应目录 -N 包含新文件到patch -u 输出统一格式(unified format),这种格 ...

  7. LINUX线程初探

     LINUX程序设计最重要的当然是进程与线程.本文主要以uart程序结合键盘输入控制uart的传输. 硬件平台:树莓派B+ 软件平台:raspberry 须要工具:USB转TTL(PL2303)+ ...

  8. vim与sublime

    vim与sublime 对程序员来说,写代码是再熟悉不过的事情了,windows系统自带有记事本软件,能写写小规模的代码,可是代码量大了,它的局限性就暴露得很明显了:没有语法高亮,没有自动提示,不支持 ...

  9. sublime3 支持 jsx 语法

    添加几个插件即可在js中快速写html babel 可以识别React,并高亮显示ES6 command+shift+p -> install package -> babel 使用 在打 ...

  10. nginx教程1:location 匹配规则

    worker_process # 表示工作进程的数量,一般设置为cpu的核数 worker_connections # 表示每个工作进程的最大连接数 server{} # 块定义了虚拟主机 liste ...