当前位置: 首页 > article >正文

BZOJ3688. 折线统计(dp+ds)

https://hydro.ac/d/bzoj/p/3688

一个很显然的dp, f ( x , k , 0 / 1 ) f(x,k,0/1) f(x,k,0/1) 表示现在末尾的数为 x x x,已经有 k k k 段线段,之前一直在上/下的方案数,转移显然。

然后前面一维我们遍历 i i i 时只会修改一个 x x x,同时查询其他前后缀的和,那直接树状数组即可。

#include<bits/stdc++.h>
using namespace std;
#ifdef LOCAL
 #define debug(...) fprintf(stdout, ##__VA_ARGS__)
#else
 #define debug(...) 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 100007
#define N 50010
struct node {
	int x, y; 
}a[N]; 
int b[N], f0, f1, ans; 
int n, m, i, j, k, T, x;

struct Binay_tree {
	int cnt[N]; 
	void add(int x, int k) {
		while(x < N) (cnt[x] += k) %= mo, x += x & -x; 
	}
	int cha(int x) {
		int ans = 0; 
		while(x) (ans += cnt[x]) %= mo, x -= x & -x; 
		return ans; 
	}
	int qry(int o, int x) {
		if(!o) return cha(x); 
		return cha(n) - cha(x); 
	}
}Bin[11][2]; 

signed main()
{
	#ifdef LOCAL
	  freopen("in.txt", "r", stdin);
	  freopen("out.txt", "w", stdout);
	#endif
//	srand(time(NULL));
//	T=read();
//	while(T--) {
//
//	}
	n = read(); m = read(); 
	for(i = 1; i <= n; ++i) a[i].x = read(), a[i].y = b[i] = read(); 
	sort(a + 1, a + n + 1, [] (node x, node y) { return x.x < y.x; }); 
	sort(b + 1, b + n + 1); 
	for(i = 1; i <= n; ++i) a[i].y = lower_bound(b + 1, b + n + 1, a[i].y) - b; 
	for(i = 1; i <= n; ++i) {
		x = a[i].y; 
//		printf("# %d\n", x); 
		for(k = m; k >= 1; --k) {
			f0 = f1 = 0; 
			f0 = Bin[k][0].qry(0, x); 
			f1 = Bin[k][1].qry(1, x); 
				f0 += Bin[k - 1][1].qry(0, x); 
				f1 += Bin[k - 1][0].qry(1, x); 
			Bin[k][0].add(x, f0); Bin[k][1].add(x, f1); 
		}
		Bin[0][0].add(x, 1); Bin[0][1].add(x, 1); 
	}
	ans = Bin[m][0].cha(n) + Bin[m][1].cha(n); 
	printf("%lld", (ans % mo + mo) % mo); 
	return 0;
}




http://www.kler.cn/a/299124.html

相关文章:

  • 进网许可认证、交换路由设备检测项目更新25年1月起
  • 电子应用设计方案-63:智能升降茶几系统方案设计
  • 【Datawhale AI 冬令营】如何动手微调出自己的大模型
  • 电商数据采集电商,行业数据分析,平台数据获取|稳定的API接口数据
  • springboot449教学资源共享平台(论文+源码)_kaic
  • 揭秘区块链隐私黑科技:零知识证明如何改变未来
  • 全国计算机二级考试C语言篇3——选择题
  • FANUC 数控 A06B-6058-H227 伺服放大器
  • AlphaNovel的身份验证失败了..........
  • 第 5 篇 Chart 文件结构详解
  • Flutter动画—雷达扫描效果
  • HarmonyOS 是如何实现一次开发多端部署 -- HarmonyOS自学1
  • Java集合——Queue详解
  • ActiveMQ 的网络连接及消息回流机制
  • Redis:处理缓存穿透的两种方法
  • web项目如何部署到服务器上并运行呢?——使用Tomcat插件
  • 数学建模笔记—— 蒙特卡罗法
  • 剖析Cookie的工作原理及其安全风险
  • git查看、修改、删除属性,配置提交commit用户名和邮箱
  • shiny 的认识和理解(必要)
  • 【JavaEE】TCP协议 (TCP-传输层协议 万字详解)
  • PHP轻松创建高效收集问卷调查小程序系统源码
  • 百度Apollo打通与ROS的通信,扩展自动驾驶系统生态
  • 8.10Laplacian算子
  • YOLO-v8:yolov8配置文件(yolov8.yaml)说明
  • 自用小问答