雅礼集训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. freemark2pdf

    freemarker+ITextRenderer 生成html转pdf 博客分类: ITextRenderer ITextRenderer  网上已经有比较多的例子 写这个 但是很多都是简单的 dem ...

  2. log4j:WARN Please initialize the log4j system properly.

    在tomcat启动的时候,出现这个警告: log4j:WARN No appenders could be found for logger (org.apache.commons.digester. ...

  3. [概念理解] MVC模式和C++的实现

    [转]学习可以是一件很快乐的事,特别是当你发现以前所学的点点滴滴慢慢地能够串起来或者变成了一个环,这种感觉真好.这篇文章就这么来的. 从MVC架构开始说起吧.这两天系统了解了一下MVC架构的内容,主要 ...

  4. 【selenium+Python WebDriver API】之复选框顺序正选和顺序反选

    from selenium import webdriver from selenium.webdriver.common.by import By import os,time driver = w ...

  5. go web的简单服务器

    1)简单web服务器: package main import ( "fmt" "net/http" ) func sayHelloName(w http.Re ...

  6. TP 自动验证规则

    #自动验证 protected $_validate=array( #参数最后代表1 表示必须验证,0表示当这个字段存在的时候验证 array('username','require','账号不能为空 ...

  7. Mysql 索引增加与删除

    [1]索引 索引,通俗理解,即目录. 之前说过,计算机是对现实世界的模拟.目录应用在数据库领域,即所谓的索引. 目录的作用显而易见,所以建立索引可以大大提高检索的速度. 但是,会降低更新表的速度,如对 ...

  8. Python—发邮件总结

    来自: http://my.oschina.net/jhao104/blog/613774 1.登录SMTP服务器 首先使用网上的方法(这里使用163邮箱,smtp.163.com是smtp服务器地址 ...

  9. 【BZOJ4861】[Beijing2017]魔法咒语 矩阵乘法+AC自动机+DP

    [BZOJ4861][Beijing2017]魔法咒语 题意:别看BZ的题面了,去看LOJ的题面吧~ 题解:显然,数据范围明显的分成了两部分:一个是L很小,每个基本词汇长度未知:一个是L很大,每个基本 ...

  10. EasyPlayerPro windows播放器本地配置文件配置方法介绍

    需求背景 应EasyPlayerPro某客户需求,在EasyPlayerPro启动时,自动播放指定的url源, 不需要每次都去手动填写, 且实现自动播放,不需要手动的单击播放按钮: 为响应该需求,特增 ...