历届试题 三体攻击

标题:三体攻击

【题目描述】
三体人将对地球发起攻击。为了抵御攻击,地球人派出了 A × B × C 艘战舰,在太空中排成一个 A 层 B 行 C 列的立方体。其中,第 i 层第 j 行第 k 列的战舰(记为战舰 (i, j, k))的生命值为 d(i, j, k)。

三体人将会对地球发起 m 轮“立方体攻击”,每次攻击会对一个小立方体中的所有战舰都造成相同的伤害。具体地,第 t 轮攻击用 7 个参数 lat, rat, lbt, rbt, lct, rct, ht 描述;
所有满足 i ∈ [lat, rat],j ∈ [lbt, rbt],k ∈ [lct, rct] 的战舰 (i, j, k) 会受到 ht 的伤害。如果一个战舰累计受到的总伤害超过其防御力,那么这个战舰会爆炸。

地球指挥官希望你能告诉他,第一艘爆炸的战舰是在哪一轮攻击后爆炸的。

【输入格式】
从标准输入读入数据。

第一行包括 4 个正整数 A, B, C, m;
第二行包含 A × B × C 个整数,其中第 ((i − 1)×B + (j − 1)) × C + (k − 1)+1 个数为 d(i, j, k);
第 3 到第 m + 2 行中,第 (t − 2) 行包含 7 个正整数 lat, rat, lbt, rbt, lct, rct, ht。

【输出格式】
输出到标准输出。

输出第一个爆炸的战舰是在哪一轮攻击后爆炸的。保证一定存在这样的战舰。

【样例输入】
2 2 2 3
1 1 1 1 1 1 1 1
1 2 1 2 1 1 1
1 1 1 2 1 2 1
1 1 1 1 1 1 2

【样例输出】
2

【样例解释】
在第 2 轮攻击后,战舰 (1,1,1) 总共受到了 2 点伤害,超出其防御力导致爆炸。

【数据约定】
对于 10% 的数据,B = C = 1;
对于 20% 的数据,C = 1;
对于 40% 的数据,A × B × C, m ≤ 10, 000;
对于 70% 的数据,A, B, C ≤ 200;
对于所有数据,A × B × C ≤ 10^6, m ≤ 10^6, 0 ≤ d(i, j, k), ht ≤ 10^9。


资源约定:
峰值内存消耗(含虚拟机) < 256M
CPU消耗  < 2000ms


请严格按要求输出,不要画蛇添足地打印类似:“请您输入...” 的多余内容。

注意:
main函数需要返回0;
只使用ANSI C/ANSI C++ 标准;
不要调用依赖于编译环境或操作系统的特殊函数。
所有依赖的函数必须明确地在源文件中 #include <xxx>
不能通过工程设置而省略常用头文件。

提交程序时,注意选择所期望的语言类型和编译器类型。

思路:别看着题目很庞大,实现就是一步步跟着走而已,不能慌

#include<stdio.h>
#include<iostream>
#include<string>
using namespace std;
#define NUM 1000001
int fx[NUM];//存储生命值 
int key=0;//判断是否找到了结果,可节省时间 
int result=0; 
int A,B,C;//意义同题目 
void sam(int a,int b,int c,int d,int e,int f,int g,int num){
	int i,j,k,sa;
	for(i=a;i<=b;i++){
		for(j=c;j<=d;j++){
			for(k=e;k<=f;k++){
				sa=B*C*(i-1)+C*(j-1)+(k-1);
				if(sa<A*B*C){
					fx[sa]-=g;
					if(fx[sa]<0){
						result=num+1;
						key=1;
						return;
					}	
				}
				
			}
		}
	}
}
int main(){
	int m,i;
	scanf("%d %d %d %d",&A,&B,&C,&m);
	for(i=0;i<A*B*C;i++){
		scanf("%d",&fx[i]);
	} 
	int lat,rat,lbt,rbt,lct,rct,ht;
	for(i=0;i<m;i++){
		
		scanf("%d %d %d %d %d %d %d",&lat,&rat,&lbt,&rbt,&lct,&rct,&ht);
		if(key==1)
			continue;
		sam(lat,rat,lbt,rbt,lct,rct,ht,i);
	}
	printf("%d",result);
	return 0;
} 

 

三体运动是指三个物体之间的相互作用和运动,可以使用Python进行模拟。其中,可以使用scipy库中的odeint函数求解微分方程组,来模拟三体运动。 以下是一个简单的三体运动模拟代码示例: ```python import numpy as np from scipy.integrate import odeint # 定义微分方程组 def three_body_equations(w, t, G, m1, m2, m3): x1, y1, vx1, vy1, x2, y2, vx2, vy2, x3, y3, vx3, vy3 = w r12 = np.sqrt((x2 - x1)**2 + (y2 - y1)**2) r13 = np.sqrt((x3 - x1)**2 + (y3 - y1)**2) r23 = np.sqrt((x3 - x2)**2 + (y3 - y2)**2) dx1dt = vx1 dy1dt = vy1 dvx1dt = G * m2 * (x2 - x1) / r12**3 + G * m3 * (x3 - x1) / r13**3 dvy1dt = G * m2 * (y2 - y1) / r12**3 + G * m3 * (y3 - y1) / r13**3 dx2dt = vx2 dy2dt = vy2 dvx2dt = G * m1 * (x1 - x2) / r12**3 + G * m3 * (x3 - x2) / r23**3 dvy2dt = G * m1 * (y1 - y2) / r12**3 + G * m3 * (y3 - y2) / r23**3 dx3dt = vx3 dy3dt = vy3 dvx3dt = G * m1 * (x1 - x3) / r13**3 + G * m2 * (x2 - x3) / r23**3 dvy3dt = G * m1 * (y1 - y3) / r13**3 + G * m2 * (y2 - y3) / r23**3 return dx1dt, dy1dt, dvx1dt, dvy1dt, dx2dt, dy2dt, dvx2dt, dvy2dt, dx3dt, dy3dt, dvx3dt, dvy3dt # 定义初始状态和参数 w0 = [1, 0, 0, 6, -1, 0, 0, -6, 0, 0, 0, 0] t = np.linspace(0, 10, 1000) G = 1 m1 = 1 m2 = 1 m3 = 1 # 求解微分方程组 wsol = odeint(three_body_equations, w0, t, args=(G, m1, m2, m3)) # 绘制轨迹图 import matplotlib.pyplot as plt plt.plot(wsol[:,0], wsol[:,1], label='Body 1') plt.plot(wsol[:,4], wsol[:,5], label='Body 2') plt.plot(wsol[:,8], wsol[:,9], label='Body 3') plt.legend() plt.show() ```
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值