1006 Xor sum
题目意思
给出一个长度为n的数组,和一个数k,问异或和不小于k的连续子序列最短是谁
题目思路
我们要求的是区间[l,r]的异或和为k,那么我们可以将它分解成区间[1,l-1]和区间[1,r]的异或值。
那么我们就可以枚举右端点r,将1到r-1的前缀异或和都放进字典树。
我们每次放入一个新的异或和时,我们沿着路径涂上这个新异或和的右端点。
我们在查找时,取出k的当前二进制位和[1,r]异或和的当前二进制位。
如果k的当前位为0,那么我们要往字典树上与[1,r]异或和的当前位值相同的方向走。
为1时,往相反的位置走就好了。
k当前位为0时,往相反方向走最终取出来的数与[1,r]异或的值x一定大于k(之前的位置异或后都相同,当k的当前位为0,而x的当前位为1),所以我们也要记录反方向取出来的数的原区间右端点。
ac代码
#include <stdio.h>
#include <iostream>
#include <algorithm>
#include <math.h>
#include <string.h>
#include <vector>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include <utility>
#define pi 3.1415926535898
#define ll long long
#define lson rt<<1
#define rson rt<<1|1
#define eps 1e-6
#define ms(a,b) memset(a,b,sizeof(a))
#define legal(a,b) a&b
#define print1 printf("111\n")
#define pb(x) push_back(x)
#define pair4 pair<pair<int,int>,pair<int,int> >
#define fi first
#define se second
using namespace std;
const int maxn = 500000+10;
const int inf = 1e9+10;
const ll llinf =1e18+10;
const ll mod = 1e9+7;
int a[40*maxn][4];
int s[maxn];
int tag[40*maxn];
int n,k;
int res;
int cnt;
void inserts(int x,int num)
{
int u=0;
for(int j=31;j>=0;j--)
{
int c=(x>>j)&1;
if(!a[u][c])
{
a[u][c]=++cnt;tag[cnt]=-1;
a[cnt][1]=a[cnt][0]=0;
}
u=a[u][c];
tag[u]=max(tag[u],num);
}
}
int finds(int x)
{
int u=0;
for(int j=31;j>=0;j--)
{
int c=(x>>j)&1;
int kk=(k>>j)&1;
if(!(kk&1))
{
if(a[u][c^1])res=max(res,tag[a[u][c^1]]);//记录相反方向的值的右端点,更新答案
u=a[u][c];
}
else
u=a[u][c^1];
if(!u)
{
return -1;
}
}
return tag[u];
}
int main()
{
int _;
scanf("%d",&_);
while(_--)
{
scanf("%d%d",&n,&k);
for(int i=1;i<=n;i++)
{
scanf("%d",&s[i]);
}
s[0]=0;
for(int i=1;i<=n;i++)
s[i]^=s[i-1];
for(int i=0;i<=cnt;i++)
{
a[i][1]=0;
a[i][0]=0;
tag[i]=0;
}
cnt=0;
int ansl=-1,ansr=n+1;
for(int i=0;i<=n;i++)//从0开始储存区间 防止[1,1]区间不计入答案
{
res=-1;
int x=finds(s[i]);
if(x!=-1)
{
res=max(x,res);
}
if(i-res<ansr-ansl&&res!=-1)
{
ansl=res;
ansr=i;
}
inserts(s[i],i);
}
if(ansl!=-1)
printf("%d %d\n",ansl+1,ansr);
else
printf("-1\n");
}
}
}
1008 Maximal submatrix
题目意思
给出一个n*m的矩阵,求每一列的值从上往下不递减的最大子矩阵
题目思路
我们先枚举每一列以每个位置为起点,求出不递减串的长度dp[i][j]。
我们先单看每一行。
我们将dp[i][j]的值看作该位置的长为1的矩形的高度,我们要求的最大子矩阵其实也就是这些矩形构成的最大矩形,这个问题我们用单调栈就能解决的了。
然后我们每一行统计完毕后维护一个最大值就好了。
ac代码
int a[maxn][maxn];
ll dp[maxn][maxn];
int q[maxn],num[maxn];
int main()
{
int _;
scanf("%d",&_);
while(_--)
{
int n,m;
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
dp[i][j]=1;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
scanf("%d",&a[i][j]);
}
}
for(int j=1;j<=m;j++)
{
for(int i=n;i>=1;i--)
{
if(a[i][j]<=a[i+1][j])
{
dp[i][j]=max(dp[i][j],dp[i+1][j]+1);
}
}
}
ll ans=1;
for(int i=1;i<=n;i++)
dp[i][m+1]=0;
for(int i=1;i<=n;i++)
{
int x=n-i+1;
int head=1,tail=0;
for(int j=1;j<=m+1;j++)
{
int temp=0;
while(head<=tail&&q[tail]>dp[i][j])
{
temp+=num[tail];
ans=max(ans,q[tail]*(ll)temp);
tail--;
}
q[++tail]=dp[i][j];
num[tail]=temp+1;
}
}
printf("%lld\n",ans);
}
}
1010 zoto
题目思路
我们用莫队处理x轴,根据坐标排序,把区间里每个横坐标对应的纵坐标加上1。
然后我们单看y轴,我每次询问就相当于问y轴的[l,r]区间里面有多少个fi序列里的点。
我们将区间分块,sum记录每个区间里面大于零的点的个数,num表示每个位置的值,
根据这两个数组我们就能求出答案了
y轴分块k选316左右(sqrt(100000))
ac代码
const int maxn = 2e5+10;
const int inf = 1e9+10;
const ll llinf = 1e18+10;
const ll mod = 1e9+7;
int n,q;
int k=313;
int siz,len;
int a[maxn];
int belong[maxn];
int num[maxn],sum[maxn];
int now;
int ans[maxn];
struct node
{
int l,r,L,R,id;
}e[maxn];
bool cmp(node a,node b)
{
return (belong[a.l]^belong[b.l])?belong[a.l]<belong[b.l]:((belong[a.l]&1)?a.r<b.r:a.r>b.r);
}
void add(int x)
{
num[x]++;
if(num[x]==1)sum[x/k]++;
}
void del(int x)
{
num[x]--;
if(num[x]==0)sum[x/k]--;
}
int calc(int x)
{
int ans=0;
for(int i=0;i<x/k;i++)
ans+=sum[i];
for(int i=(x/k)*k;i<=x;i++)
ans+=(num[i]!=0);
return ans;
}
int main()
{
int _;
scanf("%d",&_);
while(_--)
{
scanf("%d%d",&n,&q);
ms(sum,0);
ms(num,0);
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
}
siz=sqrt(n);
len=ceil((double)n/siz);
for(int i=1;i<=siz;i++)
{
for(int j=1;j<=len;j++)
{
belong[j+(i-1)*len]=i;
}
}
for(int i=1;i<=q;i++)
{
scanf("%d%d%d%d",&e[i].l,&e[i].L,&e[i].r,&e[i].R);
e[i].id=i;
}
sort(e+1,e+1+q,cmp);
int l=1,r=0;
now=0;
for(int i=1;i<=q;i++)
{
int ql=e[i].l;
int qr=e[i].r;
while(l<ql)del(a[l]),l++;
while(l>ql)add(a[--l]);
while(r>qr)del(a[r]),r--;
while(r<qr)add(a[++r]);
ans[e[i].id]=calc(e[i].R)-calc(e[i].L-1);
}
for(int i=1;i<=q;i++)
{
printf("%d\n",ans[i]);
}
}
}