[JSOI2008]球形空间产生器
Code
经过简单的转换后高斯消元即可
好久没打过高斯消元了,最近的还是 \text{two years ago......}
#include <cstdio> #include <iostream> #define RE register using namespace std; int n; double a[15][15], b[15][15]; int main() { scanf("%d", &n); for(RE int i = 1; i <= n + 1; i++) { for(RE int j = 1; j <= n; j++) scanf("%lf", &b[i][j]); if (i > 1) for(RE int j = 1; j <= n; j++) a[i - 1][n + 1] += -1.0 * (b[i - 1][j] + b[i][j]) * (b[i - 1][j] - b[i][j]), a[i - 1][j] = -2.0 * (b[i - 1][j] - b[i][j]); } for(RE int i = 1; i <= n; i++) { for(RE int j = i; j <= n; j++) if (a[i][j] != 0) { for(RE int k = 1; k <= n + 1; k++) swap(a[i][k], a[j][k]); break; } for(RE int j = i + 1; j <= n + 1; j++) a[i][j] /= a[i][i]; a[i][i] = 1.0; for(RE int j = 1; j <= n; j++) if (i ^ j) { double p = a[j][i]; for(RE int k = 1; k <= n + 1; k++) a[j][k] -= a[i][k] * p; } } for(RE int i = 1; i <= n; i++) printf("%.3lf ", a[i][n + 1]); }
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 深入理解 Mybatis 分库分表执行原理
· 如何打造一个高并发系统?
· .NET Core GC压缩(compact_phase)底层原理浅谈
· 现代计算机视觉入门之:什么是图片特征编码
· .NET 9 new features-C#13新的锁类型和语义
· Sdcb Chats 技术博客:数据库 ID 选型的曲折之路 - 从 Guid 到自增 ID,再到
· 语音处理 开源项目 EchoSharp
· 《HelloGitHub》第 106 期
· Spring AI + Ollama 实现 deepseek-r1 的API服务和调用
· 使用 Dify + LLM 构建精确任务处理应用