标签:mod P20 using return name ble strong style turn
题目链接:
Luogu:https://www.luogu.org/problem/P2054
bzoj:https://www.lydsy.com/JudgeOnline/problem.php?id=1965
找到规律:
上一次在x位置,下一次就会在 x*2%(n+1) 位置
那么就是要求:
x * (2^m) = L(mod n+1)
第一反应exgcd,但是突然感觉有点不太优秀。
因为n是偶数(题目说了)
所以2和n+1为互质
所以设 z * 2 =1(mod n+1)
z就为2在mod n+1下的逆元
z * 2 +(n+1)*y = 1
当 y=-1 时得出 z 最小为 n/2+1
所以 x=(n/2+1)^ m *L (mod n+1)
用一个快速幂和一个快速乘就可以切了。
不用快速乘中间会爆成负数。
代码如下:
#include<bits/stdc++.h> #define ll long long using namespace std; const ll maxn=1e10+10; ll n,m,l; inline ll qm(ll x,ll y){ ll mod=(n+1),ans=0; x%=mod; while(y){ if(y&1) ans=(ans+x)%mod; x=(x+x)%mod; y>>=1; } return ans%mod; } inline ll qp(ll x,ll y){ ll mod=(n+1),ans=1; x%=mod; while(y){ if(y&1) ans=qm(ans,x)%mod; x=qm(x,x)%mod; y>>=1; } return ans%mod; } int main() { scanf("%lld%lld%lld",&n,&m,&l); ll ans=qp(n/2+1,m); ans=qm(ans,l); printf("%lld\n",ans); //system("pause"); return 0; }
luogu_P2054 bzoj 1965 洗牌 【题解】 快速幂 快速乘
标签:mod P20 using return name ble strong style turn
原文地址:https://www.cnblogs.com/ChrisKKK/p/11442602.html