6.A+B(Ⅰ)

题目

Description

学过《计算机科学导论》的你应该熟悉“二进制运算”吧?

和十进制不同的是:二进制运算“逢二进一”。下面举一个二进制加法的运算实例:

   11101
+    110
--------
  100011

下面请你模拟这个过程。

Input

第一行输入一个正整数T,表示接下来有T组数据;

接下来T行,每组输入两个二进制串a和b,中间用空格隔开,保证它们长度并且没有前导0。

Output

对于每组数据,请按模拟二进制加法,按题目描述的格式输出正确的运算结果,注意换行,没有多余的空格和换行。

思路

思路:本想二进制转十进制,相加再转二进制,但是没注意到数据的大小(不是10^5, 而是2^10^5)。故应该用两个数组倒序存储数据,之后相加。注意char与int相加的区别。

小坑:注意数据大小,相加的情况要考虑全,打印的格式问题。

改进:可以使用函数简写正存、倒叙、打印部分。

代码

#include<stdio.h>
#include<string.h>
#include<math.h>
char arr1[1000010]={0};
char arr2[1000010]={0};
char arr3[1000010]={0};
int main()
{
	int t = 0; 
	int wei1,wei2,wei3;
	scanf("%d",&t) ;
	char ch; ch = getchar();
	while(t--)//t组用例 
	{
		for(wei1 = 0;wei1 < 1000010; wei1++)//正存1 
		{
			ch = getchar();
			if(ch ==' ') break;
			arr1[wei1] = ch;
		}
		for(wei2 = 0;wei2 < 1000010; wei2++)//正存2 
		{
			ch = getchar();
			if(ch == '\n') break;
			arr2[wei2] = ch;
		}
		for(int i = 0; i < wei1/2; i++)//倒存1 
		{
			char m = arr1[wei1-1-i];
			arr1[wei1-1-i] = arr1[i] ;
			arr1[i] = m;
		}
		for(int i = 0; i < wei2/2; i++)//倒存2 
		{
			char m = arr2[wei2-1-i];
			arr2[wei2-1-i] = arr2[i];
			arr2[i] = m;
		}
		
		(wei1>wei2)?(wei3=wei1):(wei3=wei2);//算 
		for(int i = 0; i < wei3+1; i++)
		{
			char c = arr1[i]+arr2[i]+arr3[i];
			if(c==96||c==48) arr3[i]='0'; 
			else if(c==97||c==49) arr3[i]='1';
			else if(c==98||c=='2'||c==-110) 
			{
				arr3[i] = '0';arr3[i+1] = '1';
			}
			else if(c==99||c==-109) 
			{
				arr3[i] = '1';arr3[i+1] = '1';
			}
		}
			if(arr3[wei3] == 49) wei3++;
			
	   		int j = 0;//打印 
			for(int i = 0; i < wei3+2; i++)
			{	
				if(i<wei3+2-wei1) printf(" ");
				else  
				{   
					printf("%c",arr1[wei1-1-j]); j++;
				}
			 } printf("\n");  j = 0;
			for(int i = 0; i < wei3+2; i++)
			{	
				if(i == 0) printf("+");
				else if(i < wei3+2-wei2) printf(" ");
				else  
				{
					printf("%c",arr2[wei2-1-j]); j++;
				}
			 } printf("\n");
			for(int i = 0; i < wei3+2 ;i++) 
			{
				printf("-");
			}printf("\n");printf(" ");printf(" ");
			for(int i = 0; i < wei3; i++) 
			{
				printf("%c",arr3[wei3-1-i]);
			}printf("\n");
	for(int i = 0; i < 1000010; i++)//重置 
	{
		arr1[i] = 0;arr2[i] = 0;arr3[i] = 0;
	}	
	}
	return 0;
}

% 1. 读取数据(假设数据已保存为Excel文件,路径需替换) % 注意:Excel中第一行为横向坐标x(0.00,0.02,...,4.00),第一列为纵向坐标y(0.00,0.02,...,1.80) data = xlsread('附件.xlsx', 'Sheet1'); % 读取整个数据区域 % 2. 提取区域Ⅰ的x、y、z数据(假设区域Ⅰ:x∈[0,2.00] NM,y∈[0,1.00] NM) % 2.1 提取x坐标(横向,由西向东)和y坐标(纵向,由南向北) x_all = data(1, 3:end); % 第一行第3列起为横向坐标(0.00,0.02,...) y_all = data(3:end, 2); % 第三行起第2列为纵向坐标(0.00,0.02,...) % 2.2 确定区域Ⅰ的索引(x≤2.00,y≤1.00) x_idx = x_all <= 2.00; % x在区域Ⅰ内的索引 y_idx = y_all <= 1.00; % y在区域Ⅰ内的索引 % 2.3 提取区域Ⅰ的x、y向量 x = x_all(x_idx); y = y_all(y_idx); n_x = length(x); % x方向数据量 n_y = length(y); % y方向数据量 % 2.4 提取区域Ⅰ的深度数据z(对应x和y的网格) z_grid = data(3:end, 3:end); % 深度数据网格(行对应y,列对应x) z_region = z_grid(y_idx, x_idx); % 区域Ⅰ的深度网格 % 2.5 转换为向量(用于最小二乘拟合) x_vec = repmat(x, n_y, 1); % x向量(n_y×n_x,与z对应) y_vec = repmat(y', 1, n_x); % y向量(n_y×n_x,与z对应) z_vec = z_region(:); % z向量(列向量) x_vec = x_vec(:); y_vec = y_vec(:); % 3. 最小二乘法平面拟合 % 构造设计矩阵X = [1, x, y] X = [ones(length(z_vec), 1), x_vec, y_vec]; % 求解系数a,b,c(z = a + b*x + c*y) coeff = (X' * X) \ (X' * z_vec); % 最小二乘解 a = coeff(1); b = coeff(2); c = coeff(3); % 4. 输出拟合结果 fprintf('平面拟合方程:z = %.4f + %.4f*x + %.4f*y\n', a, b, c); % 5. 可视化拟合结果(修正linspace范围) figure; % 绘制原始数据散点 scatter3(x_vec, y_vec, z_vec, 10, 'b', 'filled'); hold on; % 绘制拟合平面:基于区域Ⅰ的x、y实际范围生成网格(修正处) [X_grid, Y_grid] = meshgrid(linspace(min(x), max(x), 50), linspace(min(y), max(y), 50)); Z_fit = a + b*X_grid + c*Y_grid; surf(X_grid, Y_grid, Z_fit, 'FaceAlpha', 0.5); % 半透明平面,避免遮挡原始数据 xlabel('横向坐标x(NM)'); ylabel('纵向坐标y(NM)'); zlabel('海水深度z(m)'); title('区域Ⅰ海水深度平面拟合(最小二乘法)'); legend('原始数据', '拟合平面'); grid on; hold off;修改该段代码使其能够正确运行
最新发布
07-21
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值