ZJOI2013K大数查询bzoj3110

树套树,写的太弱了15s才过,不断维护就好啦,权值套区间,似乎我的常熟很渣……后来才发现本题数据是没有负数的嗯嗯喵,学长说:T_T省选第一题辣么简单还用讲?其实是几天前在JZOJ上面交过了……然后去刷了下BZOJ……

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define LL long long
using namespace std;
LL read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int a,b,c;
int n,m,sz;
int root[200005];
int ls[20000005],rs[20000005],sum[20000005],lazy[20000005];

void push(int k,int l,int r)
{
	if(!lazy[k]||l==r)return;
	if(!ls[k])ls[k]=++sz;
	if(!rs[k])rs[k]=++sz;
	lazy[ls[k]]+=lazy[k];lazy[rs[k]]+=lazy[k];
	int mid=(l+r)>>1;
	sum[ls[k]]+=(mid-l+1)*lazy[k];
	sum[rs[k]]+=(r-mid)*lazy[k];
	lazy[k]=0;
}
void modify(int &k,int l,int r,int a,int b)
{
	if(!k)k=++sz;
	push(k,l,r);
	if(l==a&&r==b)
	{
		sum[k]+=r-l+1;
		lazy[k]++;
		return;
	}
	int mid=(l+r)>>1;
        if(b<=mid)modify(ls[k],l,mid,a,b);
	else if(a>mid)modify(rs[k],mid+1,r,a,b);
	else
	{
		modify(ls[k],l,mid,a,mid);modify(rs[k],mid+1,r,mid+1,b);
	}
	sum[k]=sum[ls[k]]+sum[rs[k]];
}

int query(int k,int l,int r,int a,int b)
{
        if (!k) return 0;
        push(k,l,r);
        if (l==a&&r==b) return sum[k];
        int mid=(l+r)>>1;
        if (b<=mid) return query(ls[k],l,mid,a,b);
        else if (a>mid) return query(rs[k],mid+1,r,a,b);
        else return query(ls[k],l,mid,a,mid)+query(rs[k],mid+1,r,mid+1,b);
}
/*
void ins()
{
        int k=1,l=1,r=n;
        while (l!=r)
        {
                int mid=(l+r)>>1;
                modify(root[k],1,n,a,b);
                if (c<=mid) r=mid,k=k<<1;
                else l=mid+1,k=k<<1+1;
        }
        modify(root[k],1,n,a,b);
}
*/
void ins()
{
    int k=1,l=1,r=n;
	while(l!=r)
	{
		int mid=(l+r)>>1;
		modify(root[k],1,n,a,b);
		if(c<=mid)r=mid,k=k<<1;
		else l=mid+1,k=k<<1|1;
	}
	modify(root[k],1,n,a,b);
}
/*
int solve()
{
        int k=1,l=1,r=n;
        while (l!=r)
        {
                int mid=(l+r)>>1;
                int t=query(root[k<<1],1,n,a,b);
                if (t>=c) r=mid,k=k<<1;
                else l=mid+1,k=k<<1+1,c=c-t;
        }
        return l;
}*/
int solve()
{
	int l=1,r=n,k=1;
	while(l!=r)
	{
		int mid=(l+r)>>1;
		int t=query(root[k<<1],1,n,a,b);
		if(t>=c)r=mid,k<<=1;
		else l=mid+1,k=k<<1|1,c-=t;
	}
	return l;
}

int main()
{
//        freopen("data.in","r",stdin);
 //       freopen("data.out","w",stdout);
        n=read();m=read();
        while (m--)
        {
                int f=read();a=read();b=read();c=read();
                if (f==1)
                {
                        c=n-c+1;ins();
                }
                else printf("%d\n",n-solve()+1);
        }
        return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
题目描述 有一个 $n$ 个点的棋盘,每个点上有一个数字 $a_i$,你需要从 $(1,1)$ 走到 $(n,n)$,每次只能往右或往下走,每个格子只能经过一次,路径上的数字和为 $S$。定义一个点 $(x,y)$ 的权值为 $a_x+a_y$,求所有满足条件的路径中,所有点的权值和的最小值。 输入格式 第一行一个整数 $n$。 接下来 $n$ 行,每行 $n$ 个整数,表示棋盘上每个点的数字。 输出格式 输出一个整数,表示所有满足条件的路径中,所有点的权值和的最小值。 数据范围 $1\leq n\leq 300$ 输入样例 3 1 2 3 4 5 6 7 8 9 输出样例 25 算法1 (形dp) $O(n^3)$ 我们可以先将所有点的权值求出来,然后将其看作是一个有权值的图,问题就转化为了在这个图中求从 $(1,1)$ 到 $(n,n)$ 的所有路径中,所有点的权值和的最小值。 我们可以使用形dp来解决这个问题,具体来说,我们可以将这个图看作是一棵,每个点的父节点是它的前驱或者后继,然后我们从根节点开始,依次向下遍历,对于每个节点,我们可以考虑它的两个儿子,如果它的两个儿子都被遍历过了,那么我们就可以计算出从它的左儿子到它的右儿子的路径中,所有点的权值和的最小值,然后再将这个值加上当前节点的权值,就可以得到从根节点到当前节点的路径中,所有点的权值和的最小值。 时间复杂度 形dp的时间复杂度是 $O(n^3)$。 C++ 代码 算法2 (动态规划) $O(n^3)$ 我们可以使用动态规划来解决这个问题,具体来说,我们可以定义 $f(i,j,s)$ 表示从 $(1,1)$ 到 $(i,j)$ 的所有路径中,所有点的权值和为 $s$ 的最小值,那么我们就可以得到如下的状态转移方程: $$ f(i,j,s)=\min\{f(i-1,j,s-a_{i,j}),f(i,j-1,s-a_{i,j})\} $$ 其中 $a_{i,j}$ 表示点 $(i,j)$ 的权值。 时间复杂度 动态规划的时间复杂度是 $O(n^3)$。 C++ 代码
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值