求root(N,k)

题目描述
N<k时,root(N,k) = N,否则,root(N,k) = root(N’,k)。N’为N的k进制表示的各位数字之和。输入x,y,k,输出root(x^y,k)的值 (这里^为乘方,不是异或),2=<k<=16,0<x,y<2000000000,有一半的测试点里 x^y 会溢出int的范围(>=2000000000)

输入描述:
每组测试数据包括一行,x(0<x<2000000000), y(0<y<2000000000), k(2<=k<=16)

输出描述:
输入可能有多组数据,对于每一组数据,root(x^y, k)的值

示例1
输入
4 4 10
输出
4

题目解析:因为x^y会溢出,所以直接使用x的y次方时得到的N会有误差,所以不能直接求N进行%k,所以使用二分求幂法
关键在于等式 root(x*y,k) = root(root(x,k)*root(y,k),k)
所以有如下例 root(x^111,k) = root(root(x100,k)*root(x10,k)*root(x,k),k)
其中 root(x^100,k) = root(root(x10)*root(x10),k)
root(x^10,k) = root(root(x,k)*root(x,k),k)
(次数y为2进制)

代码:

#include<iostream>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<stdio.h>
#include<vector>
#include<map>

using namespace std;
//计算N值使用二分求幂 
int root(int N,int k){
	if(N < k){
		return N;
	}else{
		int sum = 0;
		while(N){
			sum += N % k;
			N /= k;
		}
		return root(sum , k);
	} 
}         
int computeroot(int x,int y,int k){  // root(x^y,k) = root(root(x,k)*root(y,k),k)
	int result = 1;
	int temp = root(x , k);   
	while(y > 0){  
		if(y % 2 == 1){  
			result = root(result * temp , k);
		}
		temp = root(temp * temp , k);
		y /= 2;
	}
	return result;
}
int main(){
	int x,y,k;
	while(cin >> x >> y >> k){
		cout << computeroot(x,y,k) <<endl;
	}
	return 0;
}
### 使用 MATLAB 的 `root` 函数解符号表达式的根 MATLAB 提供了多种方法用于解决代数方程的根问题,其中一种较为特殊的方式是利用 `root` 函数。尽管该函数通常与数值计算相关联,但在某些场景下也可以结合符号变量一起工作。 #### 符号表达式中的多项式根 如果目标是对一个由符号变量构成的一元多项式寻找所有的根,则推荐优先考虑 `solve` 或者 `vpasolve` 函数[^1]。然而,在特定条件下,比如希望返回显式形式的结果或者处理高次不可解析的情况时,可以采用如下策略间接调用 `root` 表达这些根: ```matlab syms x; p = x^3 - 6*x^2 + 11*x - 6; % 定义一个三次多项式作为示例 r = arrayfun(@(k) root(p,x,k), 1:degree(p)); % 假设 degree 是自定义获取最高幂次数目的辅助函数 disp(r); % 显示三个不同的根对象 ``` 需要注意的是,上述代码片段里的 `root` 并非直接来自基础库而是模拟实现的概念化写法,因为原生环境下的确不存在这样的通用接口可以直接作用于任意阶数的符号多项式上[^2]。 #### 数值近似替代方案 对于那些无法得到闭合解的情形,可转而助于数值手段得出具体答案。此时便能充分发挥内置版本的优势——即针对已知系数数组执行操作: ```matlab c = [1,-6,11,-6]; % 多项式展开后的各项系数组成向量 rn = roots(c); % 应用roots命令取得全部实部虚部组合而成的零点集合 disp(rn); ``` 这段脚本清晰明了地展示了如何便捷高效地定位给定多项式的各个孤立点位置[^3]。 #### 自定义封装实现隐含关系描述 假如确实有必要构建一套适用于广泛情况的新机制的话,那么可以从头设计满足需的功能模块。下面给出了一种初步设想框架图景: ```matlab function res = mysymbolicroot(expr,varargin) % MYSYMBOLICROOT Attempts to find symbolic representations of polynomial roots. % % RES = MYSYMBOLICROOT(EXPR,VARNAMES...) tries finding exact forms for all distinct zeros associated with EXPR, % assuming it represents a univariate polynomial over its first free variable unless otherwise specified via optional inputs. if nargin<1 || ~isa(expr,'sym') error('First argument must be an existing SYM instance.'); end varsyms = findsym(expr); if isempty(varsyms)||length(varsyms)>1 warning('Ambiguous or non-polynomial input detected!'); return; else targetvar = varsyms{1}; end deg = feval(symengine,'degree',expr,targetvar); res = cell(1,deg); for idx=1:deg try res{idx} = compose('RootOf(%s,%s)',char(expr),char(targetvar),num2str(idx)); catch ME fprintf(stderr,['Error occurred during construction:\n',ME.message]); end end end ``` 此 M 文件定义了一个新的命令叫做 `mysymbolicroot`,它接受至少一个参数 (代表待分析的目标公式),并允许额外附加选项指示哪个字母应该充当主导角色参与整个流程当中去[^4]。 --- ###
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值