jzoj6404. 【NOIP2019模拟11.04】B
题目描述
Description
Input
从文件b.in中读入数据. 第丬行三个正整数 n, m, K. 接下来 n 行每行 m 个正整数, 表示矩阵A.
Output
输出到文件b.out中. 不行, 两个数分别表示机大值和和.
Sample Input
3 5 2
1 5 3 3 3
4 1 3 3 4
4 2 4 4 3
Sample Output
4 20
Data Constraint
题解
从左往右扫,维护一个宽为K的区域
对于一个位置(i,j),求出bz[i][j]表示(i,j+1)~(i,j+K)之中是否有a[i][j]
那么在求以每个点为左上角时,区域内的点的纵坐标不会影响到结果
所以维护每种权值出现的行,0-->1就直接加,1-->0就是在删掉一个bz[i][j]=0的值时
只需要在删/加的时候求出一种值上的一个位置的前/后继
可以线段树,也可以用bitset的_Find_next()
然而NOIP应该不能用
所以显然手写bitset(
code
#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#define fo(a,b,c) for (a=b; a<=c; a++)
#define fd(a,b,c) for (a=b; a>=c; a--)
#define change(T,t) b[T][(t)/64]^=p[(t)%64]
#define pd(T,t) (b[T][(t)/64]&p[(t)%64])
#define min(a,b) (a<b?a:b)
#define max(a,b) (a>b?a:b)
#define low(x) ((x)&-(x))
#define Len 100000
using namespace std;
unsigned long long p[64];
int a[3001][3001];
int f[3002][3001];
bool bz[3001][3001];
int num[100001];
unsigned long long b[200001][47];
char st[72000001];
char *Ch=st;
int N,n,m,K,i,j,k,l;
long long ans1,ans2;
int getint()
{
int x=0;
while (*Ch<'0' || *Ch>'9') *++Ch;
while (*Ch>='0' && *Ch<='9') x=x*10+(*Ch-'0'),*++Ch;
return x;
}
int find(int T,int t)
{
int i,j,s=t/64;t%=64;
unsigned long long S=b[T][s];
if (t<63)
S>>=t+1;
else
S=0;
if (S)
return floor(log2(low(S))+0.1)+64*s+t+1;
fo(i,s+1,N)
if (b[T][i])
return floor(log2(low(b[T][i]))+0.1)+64*i;
return -1;
}
void add(int I,int i,int j,int s)
{
int k,l,L,R;
change(a[i][j],i);
change(a[i][j]+Len,n-i+1);
k=find(a[i][j]+Len,n-i+1);if (k!=-1) k=n-k+1;
l=find(a[i][j],i);
L=max(i-K+1,1);
R=i;
if (k!=-1) L=max(L,k+1);
if (l!=-1) R=min(R,l-K);
if (L<=R)
f[L][I]+=s,f[R+1][I]-=s;
}
int main()
{
// freopen("b53.in","r",stdin);
freopen("b.in","r",stdin);
freopen("b.out","w",stdout);
fread(st,1,72000001,stdin);
p[0]=1;
fo(i,1,63)
p[i]=p[i-1]<<1;
n=getint(),m=getint(),K=getint();N=n/64;
fo(i,1,n)
{
fo(j,1,m)
a[i][j]=getint();
}
memset(num,127,sizeof(num));
fo(i,1,n)
{
fd(j,m,1)
{
bz[i][j]=(num[a[i][j]]-j)<=K;
num[a[i][j]]=j;
}
fo(j,1,m)
num[a[i][j]]=2133333333;
}
fo(j,1,K)
{
fo(i,1,n)
if (!pd(a[i][j],i))
add(1,i,j,1);
}
fo(j,2,m-K+1)
{
fo(i,1,n)
f[i][j]=f[i][j-1];
fo(i,1,n)
{
if (!pd(a[i][j+K-1],i))
add(j,i,j+K-1,1);
if (!bz[i][j-1])
add(j,i,j-1,-1);
}
}
fo(i,1,n-K+1)
{
fo(j,1,m-K+1)
{
f[i][j]+=f[i-1][j];
ans1=max(ans1,f[i][j]);
ans2+=f[i][j];
}
}
printf("%lld %lld\n",ans1,ans2);
fclose(stdin);
fclose(stdout);
return 0;
}
jzoj6404. 【NOIP2019模拟11.04】B的更多相关文章
- 6423. 【NOIP2019模拟11.11】画
题目描述 Description Input Output Sample Input 3 2 3 3 6 5 1 2 1 3 Sample Output 15 Data Constraint 题解 迫 ...
- 6407. 【NOIP2019模拟11.05】小 D 与随机
题目描述 Description Input 第一行两个个整数 n,k. 之后 n -1 行,第 i 行两个整数 ui, vi, 表示一条树边. 保证输入的数据构成一棵树. Output 一行一个数表 ...
- 6402. 【NOIP2019模拟11.01】Cover(启发式合并)
题目描述 Description 小 A 现在想用
- 6411. 【NOIP2019模拟11.06】上网
题目描述 Description Input Output 若无解,则输出"Impossible". 否则第一行输出"Possible",第二行输出 n 个正整 ...
- 6409. 【NOIP2019模拟11.06】困难的图论(Tarjan求点双)
题目描述 Description 给定由 n 个点 m 条边组成的无向连通图,保证没有重边和自环. 你需要找出所有边,满足这些边恰好存在于一个简单环中.一个环被称为简单环,当且仅当它包含的所有点都只在 ...
- 【NOIP2019模拟11.01】Game(贪心+线段树)
Description: 小 A 和小 B 在玩一个游戏,他们两个人每人有
- How to Write and Publish a Scientific Paper: 7th Edition(科技论文写作与发表教程)(11.04更新)
How to Write and Publish a Scientific Paper: 7th Edition(科技论文写作与发表教程)(11.04更新) 重要通知: 最近开题报告已差不多告一段落, ...
- 安装qt5.3.2后,qtcreator在ubuntu 11.04无法启动的问题
在官方网站下载.run文件安装后,qtcreator启动失败,然后找到命令行启动,失败原因如下: shr@shr-Sieyuan:~/Qt5.3.2/Tools/QtCreator/bin$ ./qt ...
- ubuntu 11.04 源 更新不了,全显示ign、404
原文地址:http://blog.csdn.net/enjio/article/details/11603373 ubuntu 11.04 源 更新不了 分类: 开发相关2013-09-12 14 ...
随机推荐
- 系统分析与设计HW1
软件工程的定义 1993年,电气电子工程师学会(IEEE)给出了一个定义:"将系统化的.规范的.可度量的方法用于软件的开发.运行和维护的过程,即将工程化应用于软件开发中". 阅读经 ...
- C# 实现IDisposable
#region 实现IDisposable public void Dispose() { Dispose(true); GC.SuppressFinalize(this);//防止Finalize调 ...
- hadoop报错WARN util.NativeCodeLoader: Unable to load native-hadoop library for your platform... using builtin-java classes where applicable
19/06/14 10:44:58 WARN common.Util: Path /opt/hadoopdata/hdfs/name should be specified as a URI in c ...
- dbvisualizer安装
1. 下载DbVisualizer安装包. 2.解压. 无论是哪个版本的dbvisualizer破解版, 都可以找到安装程序(例dbvis_windows-x64_10_0_10.exe), dbvi ...
- C# 编写的webservice 怎样返回XML数据
[WebMethod] public string GetXml() { string sConStr = ConfigurationManager.ConnectionStrings["c ...
- [转帖]16nm国内最先进 兆芯展示x86 KX-6000八核处理器
16nm国内最先进 兆芯展示x86 KX-6000八核处理器 https://www.cnbeta.com/articles/tech/847125.htm 在近日的2019北京国际互联网科技博览会暨 ...
- 使用CyclicBarrier+线程池,按总页数分批次开多线程执行逻辑
通过CyclicBarrier+线程池的方式,同步的方式分页分批次并发高效处理逻辑,将总页数分成多个批次并发执行每页逻辑,每个批次处理DO_MAX_SIZE个页,每个批次等待DO_MAX_SIZE个页 ...
- P1056排坐椅
这是2008普及组真题,是一个提高—的模拟. 仔细读完题便有了思路:累放在i行能隔开wi个,比较排序wi,输出即可.所以在这里遇到了结构体排序的问题与手写cmp的问题.对于两个语法知识掌握得都不好,所 ...
- bfs(同一最短路径)
http://oj.jxust.edu.cn/contest/Problem?id=1702&pid=7 题意:现在有两个相同大小的地图,左上角为起点,右下角问终点.问是否存在同一条最短路径. ...
- 洛谷 P1541 乌龟棋 & [NOIP2010提高组](dp)
传送门 解题思路 一道裸的dp. 用dp[i][j][k][kk]表示用i个1步,j个2步,k个3步,kk个4步所获得的最大价值,然后状态转移方程就要分情况讨论了(详见代码) 然后就是一开始统计一下几 ...