为什么在Julia中，负1 (-1)被提升到一个偶数的幂，返回的结果总是相同的？
在Julia中，负1 (-1)被提升到一个偶数的幂，返回的结果总是相同的。这是因为Julia中的幂运算符（^）对于复数的处理方式是通过使用复数的极坐标形式来计算的。
复数可以表示为r * exp(iθ)，其中r是模长，θ是幅角。对于负1 (-1)，它的模长r为1，幅角θ为π。
当负1被提升到一个偶数的幂时，Julia会将其转换为极坐标形式，并根据幂运算的性质进行计算。由于幂运算的性质中有一个规则是：对于任意复数z，z^a * z^b = z^(a+b)，其中a和b是实数。
因此，当负1被提升到一个偶数的幂时，Julia会将其转换为极坐标形式的负1，即1 * exp(iπ)。由于指数函数exp(iπ)的周期性质，它的幂次方结果总是相同的。
总结起来，负1 (-1)被提升到一个偶数的幂，在Julia中返回的结果总是相同的，是因为Julia使用复数的极坐标形式进行计算，并利用幂运算的性质进行简化。
相关·内容
1回答
为什么常规重复乘法的幂算法是不有效的？
python、algorithm、numbers、exponential
作为一个问题，我不明白为什么我们需要一些算法像指数的平方或模指数计算一个数字的幂。
例如，我有这样一个重复的乘法算法：
def expt_mul(a, n):
r = 1
for i in xrange(n):
r *= a
return r
A乘以n倍，所以它的复杂度是O(n)，为什么它没有效率？
浏览 3提问于2015-02-04得票数 0
回答已采纳
2回答
有人能给我解释下指数函数背后的数学吗？
c++、exponential
我看到了指数函数的实现，如下图所示，当时我正在阅读有关constexpr的内容。我试着理解逻辑，但无法理解。
我试过Youtube和Google来寻找解释，但没有成功。所以，有人能提供一个我可以阅读或解释的消息来源吗？
constexpr float exp(float x, int n)
{
return n == 0 ? 1 :
n % 2 == 0 ? exp(x * x, n / 2) :
exp(x * x, (n - 1) / 2) * x;
};
浏览 0提问于2019-07-21得票数 0
回答已采纳
2回答
计算求和的问题
python
我正在尝试计算一个求和，作为类赋值的一部分，以开发梯度方法的代码。求和是我需要找出最小值的目标函数。计算求和时，我得到了这个错误"RuntimeWarning:在double_scalars y = y + (-0.01*i + exp(-(1/x1)*((u - x2)**x3))**2)中遇到的无效值“。 为什么会发生这种情况？我怎么才能修复它？ from sympy import *
from numpy import *
x1 = 1
x2 = 50
x3 = 1
def f(x1,x2,x3):
y = 0
for i in range(1,100):
浏览 14提问于2019-05-03得票数 0
1回答
极坐标中的polar网格
python、matplotlib
我正在尝试用matplotlib绘制H原子的波函数的概率密度。我设法在笛卡尔坐标中做到了这一点，但对于以后的计算，如果我在极坐标中指定psi，效果会更好。现在，我试图使绘图工作，但它给出了错误的结果(虽然绘图的轴线应该是笛卡尔坐标)。你知道怎么解决这个问题吗？
import numpy as np
from matplotlib import pyplot as plt
import matplotlib.cm as cm
from scipy import integrate
Z = 1
a_0 = 1
pi = np.pi
n = 300
r = np.linspace(-10
浏览 0提问于2016-03-20得票数 0
2回答
如何在c中得到(3623752876 * 3623752876) % 4294434817
c、long-integer
我实现了一个快速幂算法来执行(x^y) % m操作.
由于m非常大(4294434817)，所以我使用long long来存储结果。然而，在操作过程中，long long似乎仍然不够。例如，我得到了一个(3623752876 * 3623752876) % 4294434817的负数。
有办法弄清楚吗？
浏览 3提问于2015-09-11得票数 2
回答已采纳
3回答
将矩阵提升为复数幂
c、matrix、complex-numbers、gsl、exponentiation
我正在实现一个利用GSL进行矩阵运算的库。我现在需要将事物提升到任何力量，包括想象中的事物。我已经有了处理矩阵负幂的代码，但现在我需要处理虚幂，因为我的库中的所有数字都是复数。
GSL已经有这样做的工具了吗，或者我要进入循环地狱，试图为这个创建一个算法？我需要不仅能提升到虚数，还能提升到复数，比如3+2i。由于整体上对矩阵的经验有限，我甚至不确定手动完成此操作的过程，更不用说使用计算机了。
浏览 0提问于2015-06-08得票数 3
1回答
如何找到大值的对对数？
python、pow、modular-arithmetic、exp
我想知道如何找到浮子的对对数。我的第一种方法是使用Python & C中的exp()、pow()等内置函数，但它们超出了范围错误。
然后我试着把它分解成两个部分，一个整数和其他浮点数，然后分别为它们计算10提升幂，然后乘以它们，得到结果。因此，当我试图用Python计算(a*b)时，它说长int太大，无法转换为浮动。
我最初的任务是计算antilog(x)%m &我将其转换为(a*b)%m，其中a是非常大的整数&b是浮点数。
有人能帮我解决这个问题吗？是否有适用于浮点数的模块化属性？或者是否有任何“快速”和“有效”的方法来计算反日志(X)？
浏览 4提问于2014-11-16得票数 0
4回答
我不知道其他R用户是否发现了以下问题。
在R中，我执行以下操作：
> (3/-2)^(1/3)
[1] NaN
我得到了一个NaN结果。
我以类似的方式设置：
> w<-(3/-2)
> g<-1/3
> w^g
[1] NaN
但是，如果我这样做了：
> 3/-2
[1] -1.5
> -1.5^(1/3)
[1] -1.144714
有没有人可以解释这个矛盾？
浏览 0提问于2014-05-23得票数 0
2回答
用线性同余方程在java中的实现
java、math、random
我在Java中的Random类下看到了一个LCG实现，如下所示：
/*
* This is a linear congruential pseudorandom number generator, as
* defined by D. H. Lehmer and described by Donald E. Knuth in
* <i>The Art of Computer Programming,</i> Volume 3:
* <i>Seminumerical Algorithms</i>, section 3.2.1.
*
*
浏览 1提问于2014-02-20得票数 1
2回答
复数的指数形式
java、math
我有以下问题：
我需要扩展下面的类Complex (表示复数)，并添加具有指数形式的字段，以及从complex重载方法：
class Complex{
public double re;
public double im;
public Complex(double r, double i){
re =r;
im = i;
}
public Complex add(Complex a,Complex b){
return new Complex(a.re + b.re,a.im + b.im );
}
public Complex sub(Complex a,Co
浏览 1提问于2012-11-21得票数 5
回答已采纳
1回答
用公钥存储R^2有什么用？
rsa、public-key
我认为我已经成功地反向设计了三星RSA公钥这里。然而，公钥似乎主要由模组成，但它也包含一个32位整数-1 / n[0] mod 2^32，即模的前32位字的逆以及R^2 (可能是mod n?)。
有人能解释为什么RSA公钥包含这些值吗？这些价值观能做些什么？我首先想到的也许是防范侧通道攻击，但这对公钥来说是没有意义的。
在源代码中找到了更多的信息：
/* montgomery c[] += a * b[] / R % mod */
static void montMulAdd(const RSAPublicKey *key,
uint32_t* c
浏览 0提问于2021-08-02得票数 3
回答已采纳
4回答
Python模块化结果与wolfram不同吗？
python、cryptography、rsa、modulo、wolframalpha
当我运行python 3程序时：
exp = 211
p = 199
q = 337
d = (exp ** (-1)) % ((p - 1)*(q - 1))
结果211^(-1)。
但是当我运行时，我得到了我期望的结果。
我做了一些测试输出，程序中的变量exp、p和q都是我在wolfram中使用的整数值。
我的目标是从(弱)加密整数派生私钥。如果我测试我的wolfram结果，我可以正确地解密加密的消息。
浏览 1提问于2015-11-23得票数 8
回答已采纳
3回答
如何知道复杂的结果有一个或多个真实的根源
c#、math、mathematical-expressions
我想用复数来计算减号。
因为使用Math.Pow的$(-10)^(1/3)$将得到NAN。
但是当我使用Complex.Pow时，我可以得到结果。
Complex c = Complex.Pow(-10, 1.0/3);
var r = c.Magnitude;//2.1544346900318843
但真正的结果是-2.1544346900318843。
我不知道结果是肯定的，因为震级是Abs。
我不能在所有的结果中添加-1*，因为pow可以使结果为正。
(-10)^(2/3)应得到阳性结果。
我如何判断结果是肯定的，并确定是否有真正的解决办法？
请参阅：
-
浏览 3提问于2017-04-25得票数 0
回答已采纳
2回答
遇到溢出: exp()乘以erfc()
python、numpy、scipy
我想评估以下公式：
c = exp{x}*erfc{y}
(请参阅下面代码中x和y的定义。)
问题是x和y变得很大，我得到了非常大的exp{x}值和非常小的erfc(y)值。
import numpy as np
import scipy as sci
k = 5.7e-3
D = 1.53e-8
R = 1.5e-5
r = 1e-6
t = np.linspace(0.0,12,10)
x = (r/R) + (D/(R*R) - k)*t
y = (r/(2*np.sqrt(D*t))) + np.sqrt(D*t)/R
exp_x = np.exp(x)
erfc_y = s
浏览 1提问于2019-01-15得票数 0
4回答
为什么pow(a，d，n)比**d%n快这么多？
python、performance、pypy
我正在尝试实现一个，我很困惑为什么中型数字(~7位)需要这么长时间(> 20秒)。我最终发现下面这行代码是问题的根源：
x = a**d % n
(其中a、d和n都是相似但不相等的中型数字，**是求幂运算符，%是模运算符)
然后，我尝试用以下代码替换它：
x = pow(a, d, n)
相比之下，它几乎是瞬间发生的。
对于上下文，以下是原始函数：
from random import randint
def primalityTest(n, k):
if n < 2:
return False
if n % 2 == 0:
ret
浏览 2提问于2013-01-03得票数 117
回答已采纳
2回答
在OpenCV中实现复值矩阵方程
c++、opencv、matrix
我正在尝试用OpenCV实现一个复值矩阵方程。我已经在MATLAB中建立了原型，它工作得很好。从方程开始(精确的MATLAB实现)：
kernel = exp(1i .* k .* Circ3D) .* z ./ (1i .* lambda .* Circ3D .* Circ3D)
其中
1i = complex number
k = constant (float)
Circ3D = real-valued matrix of known size
lambda = constant (float)
.* = element-wise multiplication
./ = element-
浏览 0提问于2012-11-28得票数 1
回答已采纳
1回答
如何编写灵活的多重指数拟合
scipy、curve-fitting
我想为一般函数编写一个或多或少的通用拟合函数。
$f_i = \sum_i a_i exp(-t/τ_i)$
我有一些数据。
下面是一个双指数函数的示例代码，但我希望能够用最小的代码适配来拟合单指数函数或三指数函数。
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
import numpy as np
from scipy.optimize import curve_fit
import matplotlib.pyplot as plt
t = np.linspace(0, 10, 100)
a_1 = 1
a_2 = 1
tau_1 =
浏览 0提问于2018-07-11得票数 2
回答已采纳
2回答
计算幂的最有效的算法是什么？
python-3.x、exponentiation
我正在通过Python 3模拟公钥和私钥设置的RSA协议，我必须处理巨大的指数。由于pow(base,exp)似乎没有在合理的时间内运行，所以我一直在尝试使用不同的算法，但目前似乎没有一种算法有效。
到目前为止，最有效的算法是哪一个？
浏览 1提问于2019-09-20得票数 2
回答已采纳
1回答
为什么我不能从R中的分布中积分到零附近？
r、integration、distribution
我有一个分布函数，需要用一个系数来调整它。这种风险调整需要标准化，以确保结果也是PDF格式。作为一个简单的例子，假设我们的第一个PDF是一个普通的：
Qpdf <- function(r) dnorm(r , 0 ,1)
然后我有这个风险调整：
RRA <- 0.7
norINTEG <- function(r) Qpdf(r) * r^(RRA)
a <- integrate( norINTEG , -Inf , Inf )[[1]]
realPDF <- function(r) norINTEG(r) / a
但是积分a和结果是realPDF不能处
浏览 0提问于2018-05-04得票数 0
1回答
如何用不动点数实现指数？
verilog、hdl
如何在verilog中实现一个代码，该代码解析一个指数方程，其中的数字必须表示为不动点。
例如，我在C++上有这个方程，并希望转换为Verilog或VHDL：
double y = 0.1+0.75*(1.0/(1.0+exp((x[i]+40.5)/6.0)));
其中'y‘和'x’必须是不动点数。“x”也是一个向量。
我查找了有定点的模块和库，但没有一个模块和库有指数函数。
浏览 1提问于2014-08-22得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
即时通信 IM
ICP备案
对象存储
实时音视频活动推荐
腾讯云开发者
