1007D. 航行(高消之把可以递推的丢到外面处理)

news/2024/10/7 17:08:03 标签: dp, 高斯消元

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

前面64分显然,搜个状态,然后转移就行,这不是重点。

考虑我们现在高消的状态数太大,而我们实际上需要的只有 ( x , 0 ) (x,0) (x,0) 的状态,那它们之间能不能提前把关系推出来呢?

假如当前速度是0,我们想知道到下一次速度是0的系数。而在这个过程中,肯定是只往一个方向走的,所以转移不会成环,直接做即可。

#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 998244353
#define N 510
int pw(int a, int b) {
	int ans = 1; 
	while(b) {
		if(b & 1) ans *= a; 
		a *= a; b >>= 1; 
		ans %= mo; a %= mo; 
	}
	return ans; 
}
int pw(int a) { return pw(a, mo - 2); }
inline void Mod(int &a) { if(a >= mo || a <= -mo) a %= mo; if(a < 0) a += mo; }
inline void Add(int &a, int b) { a += b; Mod(a); }
inline void Mul(int &a, int b) { Mod(b); a *= b; Mod(a); } 
int n, m, i, j, k, T, tot;
map<pair<int, int>, int>vis; 
vector<pair<int, int> >G; 
int a[N][N], p[N], fm, Ans[N], ans[N]; 

void print() {
	int i, j; 
	for(i = 1; i <= tot; ++i, debug("\n"))
		for(j = 1; j <= tot + 1; ++j) debug("%lld ", a[i][j]); 
	debug("--------------------------\n"); 
}

void Gauss(int n) {
//	print(); 
	for(i = 1; i <= n; ++i) {
		for(j = i, k = 0; j <= n; ++j) if(a[j][i]) k = j; 
		if(!k) continue; 
		for(j = 1; j <= n + 1; ++j) swap(a[k][j], a[i][j]); 
		for(j = 1, k = pw(a[i][i]); j <= n + 1; ++j) Mul(a[i][j], k); 
//		print(); 
		for(j = 1; j <= n; ++j) if(j != i && a[j][i]) {
			int x = -a[j][i]; 
			for(k = 1; k <= n + 1; ++k) Add(a[j][k], x * a[i][k]); 
		}
	}
	for(i = 1; i <= n; ++i) {
		if(!Ans[i]) continue; 
		if(!a[i][i] || !a[i][n + 1]) { ans[Ans[i]] = -1;  continue; }
		for(j = 1; j <= n; ++j) if(j != i && a[i][j]) break; 
		if(j <= n) { ans[Ans[i]] = -1; continue; }
		ans[Ans[i]] = a[i][n + 1]; 
	}
}

struct node {
	int p, c; 
	void add(node F, int P) {
		Add(p, F.p * P); Add(c, (F.c + F.p) * P % mo); 
	}
}f[N][101];

signed main()
{
	  freopen("sail.in", "r", stdin);
	  freopen("sail.out", "w", stdout);
//	srand(time(NULL));
//	T = read();
//	while(T--) {
//
//	}
	n = read(); fm = pw(100); 
	for(i = 1; i <= n; ++i) {
		k = read(); p[i] = k * fm % mo; Ans[i] = i; 
	}
	for(i = 1; i <= n; ++i) {
		memset(f, 0, sizeof(f)); 
		f[i + 1][1] = {1 - p[i], 1 - p[i]}; f[i - 1][1] = {p[i], p[i]}; 
		for(j = i + 1; j <= n; ++j)
			for(k = 1; k <= 100; ++k) {
				f[min(n + 1, j + k - 1)][k - 1].add(f[j][k], p[j]); 
				f[min(n + 1, j + k + 1)][k + 1].add(f[j][k], 1 - p[j]); 
			}
		for(j = i - 1; j >= 1; --j) 
			for(k = 1; k <= 100; ++k)  {
				f[max(0ll, j - (k - 1))][k - 1].add(f[j][k], 1 - p[j]); 
				f[max(0ll, j - (k + 1))][k + 1].add(f[j][k], p[j]); 
			}
		for(j = 1; j <= n; ++j) 
			if(i != j) a[i][j] = -f[j][0].p, Add(a[i][n + 1], f[j][0].c), debug("%lld\n", f[j][0].c);
		for(j = 0; j <= 100; ++j) Add(a[i][n + 1], f[0][j].c + f[n + 1][j].c), debug("%lld %lld\n", f[0][j].c, f[n + 1][j].c); 
		a[i][i] = 1; 
	}
	tot = n; 
//	debug("== %lld %lld\n", pw(2), pw(4)); 
//	print(); 
	Gauss(tot); 
	for(i = 1; i <= n; ++i) printf("%lld ", ans[i]); 
	return 0;
}


http://www.niftyadmin.cn/n/5693099.html

相关文章

[C语言]第十一节 函数递归一基础知识到高级技巧的全景探索

目录 11.1. 递归是什么&#xff1f; 11.1.1 递归的思想&#xff1a; 11.2 递归的限制条件 举例1&#xff1a;求n的阶乘 画图推演 举例2&#xff1a;顺序打印⼀个整数的每⼀位 画图推演 11.3. 递归与迭代 举例3&#xff1a;求第n个斐波那契数 11.1. 递归是什么&#xff…

Android阶段学习思维导图

前言 记录下自己做的一个对Android原生应用层的思维导图&#xff0c;方便个人记忆扩展&#xff1b;这里只露出二级标题。 后语 虽然有些内容只是初步了解&#xff0c;但还是记录了下来&#xff1b;算是对过去一段学习的告别。

程序化交易的接口从哪里获取?怎么选择?

炒股自动化&#xff1a;申请官方API接口&#xff0c;散户也可以 python炒股自动化&#xff08;0&#xff09;&#xff0c;申请券商API接口 python炒股自动化&#xff08;1&#xff09;&#xff0c;量化交易接口区别 Python炒股自动化&#xff08;2&#xff09;&#xff1a;获取…

Java性能调优:实战技巧与最佳实践

引言 Java作为企业级应用开发的首选语言之一&#xff0c;其性能直接影响到系统的响应速度和用户体验。性能调优是一项复杂的工作&#xff0c;涉及多个层面的知识和技术。本文将通过具体的示例&#xff0c;探讨一些常见的性能调优技巧及最佳实践。 1. 了解你的应用程序 示例&…

程序化交易与非程序化交易者盈利能力孰优孰劣

炒股自动化&#xff1a;申请官方API接口&#xff0c;散户也可以 python炒股自动化&#xff08;0&#xff09;&#xff0c;申请券商API接口 python炒股自动化&#xff08;1&#xff09;&#xff0c;量化交易接口区别 Python炒股自动化&#xff08;2&#xff09;&#xff1a;获取…

苹果电脑磁盘满了怎么清理内存?必看清理秘籍

对于很多Mac用户来说&#xff0c;随着时间的推移&#xff0c;电脑逐渐变慢并出现磁盘空间不足的提示是一件非常头疼的事情。正确理解内存和存储的区别&#xff0c;并采用有效的清理方法&#xff0c;对于保持Mac性能至关重要。本文将深入探讨如何有效清理Mac上的磁盘空间&#x…

AcWing 662:点的坐标 ← 结构体 or 三目运算符

【题目来源】https://www.acwing.com/problem/content/664/【题目描述】 给定两个保留一位小数的浮点数 X,Y&#xff0c;用来表示一个点的横纵坐标。 请你判断该点在坐标系中的位置。 【输入格式】 共一行&#xff0c;包含两个浮点数 X,Y&#xff0c;表示点的横纵坐标。【输出格…

k8s 中存储之 NFS 卷

目录 1 NFS 卷的介绍 2 NFS 卷的实践操作 2.1 部署一台 NFS 共享主机 2.2 在所有k8s节点中安装nfs-utils 2.3 部署nfs卷 2.3.1 生成 pod 清单文件 2.3.2 修改 pod 清单文件增加 实现 NFS卷 挂载的 参数 2.3.3 声明签单文件并查看是否创建成功 2.3.4 在 NFS 服务器 创建默认发布…