【NOIP2016提高A组8.12】奇袭
题目
由于各种原因,桐人现在被困在Under World(以下简称UW)中,而UW马上要迎来最终的压力测试——魔界入侵。
唯一一个神一般存在的Administrator被消灭了,靠原本的整合骑士的力量是远远不够的。所以爱丽丝动员了UW全体人民,与整合骑士一起抗击魔族。
在UW的驻地可以隐约看见魔族军队的大本营。整合骑士们打算在魔族入侵前发动一次奇袭,袭击魔族大本营!
为了降低风险,爱丽丝找到了你,一名优秀斥候,希望你能在奇袭前对魔族大本营进行侦查,并计算出袭击的难度。
经过侦查,你绘制出了魔族大本营的地图,然后发现,魔族大本营是一个N×N的网格图,一共有N支军队驻扎在一些网格中(不会有两只军队驻扎在一起)。
在大本营中,每有一个k×k(1≤k≤N)的子网格图包含恰好k支军队,我们袭击的难度就会增加1点。
现在请你根据绘制出的地图,告诉爱丽丝这次的袭击行动难度有多大。
分析
想到,可以把题目简化为在一段数中,选择一段区间,使的该区间的数连续,求方案数。
接着可以进一步转换为在一段数中,选择相邻k个数,使得这k个数中的最大值减去最小值等于k-1,求方案数。
然后考虑如何解决这个问题。
我们用分治的思想。
对于一段区间,左边界为l,右边界为r
那么这一段区间的答案\(ans(l,r)=ans(l,mid)+ans(mid+1,r)+这k个数穿过mid的方案数\)
首先知道一个合法的区间\([i,j]\),\(j-i=区间[i,j]中的最大值-区间[i,j]中的最小值\)
这里分两种情况:
最大最小值都在同一侧
现在假设都在左侧的情况,即在区间\([l,mid]\)中。右侧的情况自己思考。
先定义:
mal[x]:表示区间[x,mid]中的最大值(x在区间[l,mid]中)
mil[x]:表示区间[x,mid]中的最小值(x在区间[l,mid]中)
mar[x]:表示区间[mid+1,x]中的最大值(x在区间[mid+1,r]中)
mir[x]:表示区间[mid+1,x]中的最小值(x在区间[mid+1,r]中)
我们枚举一个i,i从mid向l移动。设j为区间的右边界
因为最大最小值都在左侧,
根据合法区间的定义,
可以轻松求出j,\(j=i+mal[i]-mil[i]\),
但是j不一定是合法的,
1、j必须在区间\([mid+1,r]\)之间
2、mar[j]必须小于mal[i],mir[j]必须大于mil[i],否则最大最小值就不都在左侧了。
右侧的情况类似。
最大最小值在异侧
现在假设最大值在右侧,即在区间\([mid+1,r]\)中;现在假设最小值在右侧,即在区间\([l,mid]\)中。
我们同样枚举一个i,i从mid向l移动。设j为区间的右边界
我们再定义两个指针z和z1从mid+1向r移动。
因为最大值在右侧,所以mal[i]应该小于mar[z],那么当mal[i]>mar[z]时,将指针z向右移;因为mal和mar都是单调的,对于当前的i,因为区间mar[z-1]一定小于mal[i],都是不合法的。
又因为最小值在左侧,所以mil[i]应该小于mir[z1],那么当mir[z1]>mil[i]时,将指针z1向右移;因为mil和mir也都是单调的,对于当前的i,因为区间mir[z1-1]一定大于mil[i],都是合法的。
于是区间[z,z1]中的数都有可能是合法的j。
再根据合法区间的定义,移项得到\(mil[i]-i=mar[j]-j\)
定义一个桶t[],
那么当z移动时,经过的点都是不合法的,就将t[mar[z]-z]减去一;
那么当z1移动时,经过的点都是合法的,就将t[mar[z1]-z1]加上一。
最后,将ans加上t[mil[i]-i]。
桶记住要清零。
另一种情况自己考虑。
#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const int maxlongint=2147483647;
const int mo=1000000007;
const int N=50005;
using namespace std;
int a[N*2],n,mal[N*2],mar[N*2],mil[N*2],mir[N*2],t[N*5],j,z,z1;
int solve(int l,int r)
{
int mid=(l+r)/2,ans=0;
if(l!=r) ans=solve(l,mid)+solve(mid+1,r);
else return 1;
for(int i=l;i<=r;i++) mar[i]=mal[i]=0,mir[i]=mil[i]=maxlongint;
for(int i=mid;i>=l;i--)
{
mil[i]=min(mil[i+1],a[i]);
mal[i]=max(mal[i+1],a[i]);
}
for(int i=mid+1;i<=r;i++)
{
mir[i]=min(mir[i-1],a[i]);
mar[i]=max(mar[i-1],a[i]);
}
//极值都在左边
for(int i=mid;i>=l;i--)
{
j=i+mal[i]-mil[i];
if((j<=mid) || (j>r)) continue;
if(mal[i]>mar[j] && mil[i]<mir[j]) ans++;
}
//极值都在右边
for(int i=mid+1;i<=r;i++)
{
j=i-mar[i]+mir[i];
if((j>=mid+1) || (j<l)) continue;
if(mar[i]>mal[j] && mir[i]<mil[j]) ans++;
}
//min在左,max在右
z=z1=mid+1;
for(int i=mid;i>=l;i--)
{
while(z<=r && mar[z]<mal[i])
{
t[mar[z]-z+N]--;
z++;
}
while(z1<=r && mir[z1]>mil[i])
{
t[mar[z1]-z1+N]++;
z1++;
}
if(t[mil[i]-i+N]>=0)
ans+=t[mil[i]-i+N];
}
for(int i=mid+1;i<=r;i++) t[mar[i]-i+N]=0;
//max在左,min在右
z=z1=mid+1;
for(int i=mid;i>=l;i--)
{
while(z1<=r && mir[z1]>mil[i])
{
t[mir[z1]+z1+N]--;
z1++;
}
while(z<=r && mal[i]>mar[z])
{
t[mir[z]+z+N]++;
z++;
}
if(t[mal[i]+i+N]>=0)
ans+=t[mal[i]+i+N];
}
for(int i=mid+1;i<=r;i++) t[mir[i]+i+N]=0;
return ans;
}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
int x,y;
scanf("%d%d",&x,&y);
a[x]=y;
}
printf("%d",solve(1,n));
}
【NOIP2016提高A组8.12】奇袭的更多相关文章
- [JZOJ4685] 【NOIP2016提高A组8.12】礼物
Description 夏川的生日就要到了.作为夏川形式上的男朋友,季堂打算给夏川买一些生日礼物.商店里一共有种礼物.夏川每得到一种礼物,就会获得相应喜悦值Wi(每种礼物的喜悦值不能重复获得).每次, ...
- 【NOIP2016提高A组8.12】礼物
题目 夏川的生日就要到了.作为夏川形式上的男朋友,季堂打算给夏川买一些生日礼物. 商店里一共有种礼物.夏川每得到一种礼物,就会获得相应喜悦值Wi(每种礼物的喜悦值不能重复获得). 每次,店员会按照一定 ...
- 【NOIP2016提高A组8.12】总结
惨败!!!! 第一题是一道神奇的期望问题. 第二题,发现"如果两个部门可以直接或间接地相互传递消息(即能按照上述方法将信息由X传递到Y,同时能由Y传递到X),我们就可以忽略它们之间的花费&q ...
- 【NOIP2016提高A组8.12】通讯
题目 "这一切都是命运石之门的选择." 试图研制时间机器的机关SERN截获了中二科学家伦太郎发往过去的一条短信,并由此得知了伦太郎制作出了电话微波炉(仮). 为了掌握时间机器的技术 ...
- JZOJ 4732. 【NOIP2016提高A组模拟8.23】函数
4732. [NOIP2016提高A组模拟8.23]函数 (Standard IO) Time Limits: 1500 ms Memory Limits: 262144 KB Detailed ...
- JZOJ 【NOIP2016提高A组集训第16场11.15】兔子
JZOJ [NOIP2016提高A组集训第16场11.15]兔子 题目 Description 在一片草原上有N个兔子窝,每个窝里住着一只兔子,有M条路径连接这些窝.更特殊地是,至多只有一个兔子窝有3 ...
- JZOJ 【NOIP2016提高A组集训第16场11.15】SJR的直线
JZOJ [NOIP2016提高A组集训第16场11.15]SJR的直线 题目 Description Input Output Sample Input 6 0 1 0 -5 3 0 -5 -2 2 ...
- 【NOIP2016提高A组集训第14场11.12】随机游走
题目 YJC最近在学习图的有关知识.今天,他遇到了这么一个概念:随机游走.随机游走指每次从相邻的点中随机选一个走过去,重复这样的过程若干次.YJC很聪明,他很快就学会了怎么跑随机游走.为了检验自己是不 ...
- 【NOIP2016提高A组集训第14场11.12】随机游走——期望+树形DP
好久没有写过题解了--现在感觉以前的题解弱爆了,还有这么多访问量-- 没有考虑别人的感受,没有放描述.代码,题解也写得歪歪扭扭. 并且我要强烈谴责某些写题解的代码不打注释的人,像天书那样,不是写给普通 ...
随机推荐
- Spring+hibernate 配置实例
转自:http://www.cnblogs.com/hongten/archive/2012/03/10/java_spring_hibernate.html 项目结构: http://www.cnb ...
- 七:flask-一些小细节
1.在局域网中,让其他电脑访问我的网站:host参数 如果设置为0.0.0.0,则在局域网中,输入当前项目所在的ip+端口就可以访问这个项目如果host设置为固定的ip,如host=‘'192.168 ...
- 【奇技淫巧】使用 SSH 转发 Sock5 流量
标题:使用 SSH 转发 Sock5 流量 日期:2018-06-27 介绍:使用 ssh 来做个 sock5 的代理,穿透到内网中做后渗透 0x01. 基本信息 在 ubuntu(10.211.55 ...
- 性能测试工具之Apache ab
一.apache ab简介 ab全称ApacheBench,是著名的Web服务器软件apache附带的一款非常简单的压力测试工具,它可以同时模拟多个并发请求,测试Web服务器最大承受压力.Apache ...
- Selenium学习之==>ActionChainsApi接口详解
ActionChains UI自动化测试过程中,经常遇到那种,需要鼠标悬浮后,要操作的才会元素出现的这种场景,那么我们就要模拟鼠标悬浮到某一个位置,做一系列的连贯操作,Selenium给我们提供了Ac ...
- PHP开发一个简单的成绩录入系统
预览界面 源码: <!DOCTYPE html> <html> <head> <meta charset="UTF-8"> < ...
- Html5 Canvas斗地主游戏
过完年来公司,没什么事,主管说研究下html5 游戏,然后主管就给了一个斗地主的demo,随后我就开始看代码, 现在我看了html5以及canvas相关知识和斗地主的demo后,自己用demo上的素材 ...
- 使用cython库对python代码进行动态编译达到加速效果及python第三方包的制作安装
1.测试代码:新建 fib.pyx # coding:utf-8 import matplotlib.pyplot as plt import numpy as np from sklearn.cl ...
- hive数据去重
Hive是基于Hadoop的一个数据仓库工具,可以将结构化的数据文件映射为一张数据库表,并提供类SQL查询功能 hive的元数据存储:通常是存储在关系数据库如 mysql(推荐) , derby(内嵌 ...
- 使用extract-text-webpack-plugin插件后报错
如果你使用的webpack是4+版本,那么尝试运行npm install extract-text-webpack-plugin@next ,即可解决问题 然而最好的解决办法是在webpack4+的版 ...