背景:
今天我才发现快速乘有
Θ
(
1
)
\Theta(1)
Θ(1)的做法。
正题:
快速乘解决的是
x
∗
y
m
o
d
  
p
x*y \mod p
x∗ymodp的问题。
有时
x
,
y
,
p
x,y,p
x,y,p特别大,
Θ ( log n ) \Theta(\log n) Θ(logn):
原理:
y
≡
0
(
m
o
d
  
2
)
:
x
∗
y
m
o
d
  
p
=
2
∗
x
∗
y
2
m
o
d
  
p
y≡0(\mod2):x*y \mod p=2*x*\frac{y}{2}\mod p
y≡0(mod2):x∗ymodp=2∗x∗2ymodp
y
≡
1
(
m
o
d
  
2
)
:
x
∗
y
m
o
d
  
p
=
(
2
∗
x
∗
⌊
y
2
⌋
+
y
)
m
o
d
  
p
y≡1(\mod2):x*y \mod p=(2*x*\lfloor\frac{y}{2}\rfloor+y)\mod p
y≡1(mod2):x∗ymodp=(2∗x∗⌊2y⌋+y)modp
写成递归或递推的形式,在每一个可以的位置上
m
o
d
  
p
\mod p
modp即可。
注意:
a
m
o
d
  
b
a \mod b
amodb的结果符号由
a
a
a确定,因此
x
∗
y
m
o
d
  
p
x*y\mod p
x∗ymodp的结果符号由
x
∗
y
x*y
x∗y确定。
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
using namespace std;
LL x,y,mod;
LL dg(LL x,LL y,LL mod)
{
if(y==1) return x;
LL op=dg(x,y>>1,mod);
if(y&1) return (op*2%mod+x)%mod; else return op*2%mod;
}
int main()
{
scanf("%lld %lld %lld",&x,&y,&mod);
int op=((x<0&&y>0)||(x>0&&y<0));
printf("%lld\n",(op?-1ll:1ll)*dg(abs(x),abs(y),abs(mod)));
}
Θ ( 1 ) \Theta(1) Θ(1):
原理:
x
∗
y
m
o
d
  
p
=
(
x
∗
y
−
⌊
x
∗
y
p
⌋
∗
p
)
m
o
d
  
p
x*y\mod p=(x*y-\lfloor\frac{x*y}{p}\rfloor*p)\mod p
x∗ymodp=(x∗y−⌊px∗y⌋∗p)modp
此时,我们的过程量用
long double
\text{long double}
long double来存储即可。
注意:
a
m
o
d
  
b
a \mod b
amodb的结果符号由
a
a
a确定,因此
x
∗
y
m
o
d
  
p
x*y\mod p
x∗ymodp的结果符号由
x
∗
y
x*y
x∗y确定。
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
typedef long double LD;
LL x,y,mod;
LL times(LL a,LL b,LL mod)
{
a%=mod,b%=mod;
return ((a*b-(LL)((LL)((LD)a/mod*b+1e-3)*mod))%mod+mod)%mod;
}
int main()
{
scanf("%lld %lld %lld",&x,&y,&mod);
int op=((x<0&&y>0)||(x>0&&y<0));
printf("%lld\n",(op?-1ll:1ll)*times(abs(x),abs(y),abs(mod)));
}