题目描述

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的更多相关文章

  1. 6423. 【NOIP2019模拟11.11】画

    题目描述 Description Input Output Sample Input 3 2 3 3 6 5 1 2 1 3 Sample Output 15 Data Constraint 题解 迫 ...

  2. 6407. 【NOIP2019模拟11.05】小 D 与随机

    题目描述 Description Input 第一行两个个整数 n,k. 之后 n -1 行,第 i 行两个整数 ui, vi, 表示一条树边. 保证输入的数据构成一棵树. Output 一行一个数表 ...

  3. 6402. 【NOIP2019模拟11.01】Cover(启发式合并)

    题目描述 Description 小 A 现在想用

  4. 6411. 【NOIP2019模拟11.06】上网

    题目描述 Description Input Output 若无解,则输出"Impossible". 否则第一行输出"Possible",第二行输出 n 个正整 ...

  5. 6409. 【NOIP2019模拟11.06】困难的图论(Tarjan求点双)

    题目描述 Description 给定由 n 个点 m 条边组成的无向连通图,保证没有重边和自环. 你需要找出所有边,满足这些边恰好存在于一个简单环中.一个环被称为简单环,当且仅当它包含的所有点都只在 ...

  6. 【NOIP2019模拟11.01】Game(贪心+线段树)

    Description: ​ 小 A 和小 B 在玩一个游戏,他们两个人每人有

  7. How to Write and Publish a Scientific Paper: 7th Edition(科技论文写作与发表教程)(11.04更新)

    How to Write and Publish a Scientific Paper: 7th Edition(科技论文写作与发表教程)(11.04更新) 重要通知: 最近开题报告已差不多告一段落, ...

  8. 安装qt5.3.2后,qtcreator在ubuntu 11.04无法启动的问题

    在官方网站下载.run文件安装后,qtcreator启动失败,然后找到命令行启动,失败原因如下: shr@shr-Sieyuan:~/Qt5.3.2/Tools/QtCreator/bin$ ./qt ...

  9. ubuntu 11.04 源 更新不了,全显示ign、404

    原文地址:http://blog.csdn.net/enjio/article/details/11603373   ubuntu 11.04 源 更新不了 分类: 开发相关2013-09-12 14 ...

随机推荐

  1. 系统分析与设计HW1

    软件工程的定义 1993年,电气电子工程师学会(IEEE)给出了一个定义:"将系统化的.规范的.可度量的方法用于软件的开发.运行和维护的过程,即将工程化应用于软件开发中". 阅读经 ...

  2. C# 实现IDisposable

    #region 实现IDisposable public void Dispose() { Dispose(true); GC.SuppressFinalize(this);//防止Finalize调 ...

  3. 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 ...

  4. dbvisualizer安装

    1. 下载DbVisualizer安装包. 2.解压. 无论是哪个版本的dbvisualizer破解版, 都可以找到安装程序(例dbvis_windows-x64_10_0_10.exe), dbvi ...

  5. C# 编写的webservice 怎样返回XML数据

    [WebMethod] public string GetXml() { string sConStr = ConfigurationManager.ConnectionStrings["c ...

  6. [转帖]16nm国内最先进 兆芯展示x86 KX-6000八核处理器

    16nm国内最先进 兆芯展示x86 KX-6000八核处理器 https://www.cnbeta.com/articles/tech/847125.htm 在近日的2019北京国际互联网科技博览会暨 ...

  7. 使用CyclicBarrier+线程池,按总页数分批次开多线程执行逻辑

    通过CyclicBarrier+线程池的方式,同步的方式分页分批次并发高效处理逻辑,将总页数分成多个批次并发执行每页逻辑,每个批次处理DO_MAX_SIZE个页,每个批次等待DO_MAX_SIZE个页 ...

  8. P1056排坐椅

    这是2008普及组真题,是一个提高—的模拟. 仔细读完题便有了思路:累放在i行能隔开wi个,比较排序wi,输出即可.所以在这里遇到了结构体排序的问题与手写cmp的问题.对于两个语法知识掌握得都不好,所 ...

  9. bfs(同一最短路径)

    http://oj.jxust.edu.cn/contest/Problem?id=1702&pid=7 题意:现在有两个相同大小的地图,左上角为起点,右下角问终点.问是否存在同一条最短路径. ...

  10. 洛谷 P1541 乌龟棋 & [NOIP2010提高组](dp)

    传送门 解题思路 一道裸的dp. 用dp[i][j][k][kk]表示用i个1步,j个2步,k个3步,kk个4步所获得的最大价值,然后状态转移方程就要分情况讨论了(详见代码) 然后就是一开始统计一下几 ...