8.22 T3 escape from whk 3(2次幂相关)

http://cplusoj.com/d/senior/p/SS240709C

考虑一个 [ l , r ] [l,r] [l,r] 区间,我们有什么策略?

性质1:我们从大到小,能选就选,肯定最优

如果按照这样子,我们可以发现选出来的数肯定长成这个样子

在这里插入图片描述

因此我们可以:

在这里插入图片描述

这样子对于一个 [ l , r ] [l,r] [l,r] 的复杂度是 O ( log ⁡ ) O(\log) O(log) 的,我们解决掉了 n u m = 0 num=0 num=0

考虑统计全局,假设我们现在求 [ 1 , r ] [1,r] [1,r],肯定长成这样:

在这里插入图片描述

而当 l l l 变化时,还是选 l l l 右边的蓝色段:

在这里插入图片描述

所以我们在固定 r r r 的时候,可以直接考虑蓝色短的贡献。

对于一个蓝色格子,它会在所有 l l l 在它前面时有贡献。也就是说 i i i 的贡献恰好为 i i i,这启示我们:

在这里插入图片描述

所以对于一个整个蓝色短,我们直接等差数列求和必然就是其贡献

这个可以递归,也可以直接dp实现

#include<bits/stdc++.h>
using namespace std;
#ifdef LOCAL
 #define debug(...) fprintf(stdout, ##__VA_ARGS__)
 #define debag(...) fprintf(stderr, ##__VA_ARGS__)
#else
 #define debug(...) void(0)
 #define debag(...) void(0)
#endif
#define int long long
inline int 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<<1)+
(x<<3)+(ch^48);ch=getchar();}return x*f;}
#define Z(x) (x)*(x)
#define pb push_back
#define fi first
#define se second
//#define M
//#define mo
#define N 300010
int n, m, i, j, k, T;
int In[N], f[N], Q, o, l, r, ans; 

int solve(int l, int r) {
	int mid = In[r]; 
	if(l > r) return 0; 
	if(mid < l) return r - l + 1; 
	int ans = r - mid + 1, k = mid - ans; 
	return ans + solve(l, k); 
}

int calc(int l, int r) {
	debug("[%lld %lld]\n", l, r); 
	return (l + r) * (r - l + 1) / 2; 
}

//int solve2(int r) {
//	int mid = In[r]; 
//	debug("r is %llld\n", r); 
//	if(1 > r) return 0; 
//	if(mid <= 1) return calc(1, r); 
//	int ans = calc(mid, r), k = mid - (r - mid + 1); 
//	return ans + solve2(k); 
//}

signed main()
{
freopen("kuhu.in", "r", stdin);
	  freopen("kuhu.out", "w", stdout);
//	srand(time(NULL));
//	T=read();
//	while(T--) {
//
//	}
	n = read(); Q = read(); o = read(); 
	for(i = 1, k = 1; i <= n; ++i) {
		if((k << 1) <= i) k <<= 1; 
		In[i] = k; 
	}
	for(i = 1; i <= Q; ++i) {
		l = read(); r = read(); 
		printf("%lld\n", solve(l, r)); 
	}
	if(!o) return printf("0"), 0; 
	f[1] = 1; f[2] = 3; ans = 4; 
	for(i = 3; i <= n; ++i) {
		j = In[i]; k = max(0ll, j - (i - j + 1)); 
		f[i] = f[k] + calc(j, i); ans += f[i]; 
	}
	printf("%lld", ans); 
	return 0;
}



  • 8
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值