算法总结之求解模线性方程组1)求解模线性方程 ax = b(mod n) 方程ax = b(mod n) -> ax = b + ny ->ax - ny = b -> ax + n (-y) =b 其中a,n,b已知。 可用扩展欧几里得来求解该方程的一组特解。 这里给出下列几个定理用来求解方...
分类:
编程语言 时间:
2014-10-30 20:51:32
阅读次数:
297
题目大意:给定n个形如xi=ki*x_pi+bi mod p的同余方程组 支持修改操作和求解操作
确实好题 感谢此题作者 顺便吐槽一下作者的Splay不加空节点太蛋疼了0.0
将每个点i的父亲设为pi 我们将会得到一座基环树林 将环上的一条边拆掉,在边的起始节点新开个域special_father记录这条边(P.S:好浪费 但是没办法)
于是我们得到了一座森林 显然可以用LCT来维护 每个节...
分类:
其他好文 时间:
2014-10-25 09:20:13
阅读次数:
284
题目大意:……简洁明了自己看
第一问快速幂
第二问扩展欧几里得
第三问BSGS
顺便一开始没看到p是质数0.0 去弄了EXBSGS的模板0.0 懒得改了
#include
#include
#include
#include
#include
#define M 1001001
using namespace std;
typedef long long ll;
typedef pa...
分类:
其他好文 时间:
2014-10-23 14:28:12
阅读次数:
239
题目大意:给定一个数列X(i+1)=(a*Xi+b)%p 求最小的i>0,使Xi=t
0.0 此题能1A真是太好了
首先讨论特殊情况
若X1=t ans=1
若a=0 ans=b==t?2:-1
若a=1 X1+b*(ans-1)==t (%p) 扩展欧几里得
令
temp=b/(a-1)
则有
(X(i+1)+temp)=a*(Xi+temp)
Xans=(X1+temp)*...
分类:
其他好文 时间:
2014-10-21 17:52:32
阅读次数:
197
扩展欧几里得,给两个点,就可以求出直线方程为 (yy-y)*x0 + (x-xx)*y0 =x*yy - y*xx,求的是在线段上的整点个数。所以就是(yy-y)*10*x0 + (x-xx)*10*y0 = x*yy - y*xx满足条件的解的个数。用exgcd搞之后求出一个解,再求出在线段上第一...
分类:
其他好文 时间:
2014-10-21 02:15:40
阅读次数:
197
青蛙的约会
Time Limit: 1000MS
Memory Limit: 10000K
Total Submissions: 91753
Accepted: 16849
Description
两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面。它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西...
分类:
其他好文 时间:
2014-10-17 20:33:41
阅读次数:
211
原题链接:hdu 1576 A/B同样是用扩展的欧几里得算法。$A = 9973k+n = xB$从而转化为:$xB-9973k=n$求解$x$即可。具体扩展欧几里得算法请参考:hdu 2669 Romantic代码如下: 1 #include 2 #include 3 #include 4 ...
分类:
其他好文 时间:
2014-10-06 16:16:00
阅读次数:
214
题目链接:HDU 1576 A/B
中文题,
思路:设X=(A/B)%9973。A/B=k_1*9973+X。A=B*k_1*9973+x*B。n=A%9973,A=k_2*9973+n。k_2*9973+n=B*k_1*9973+x*B
B*X ≡ n mod 9973 就是转化为 求B关于n模9973 的逆元。gcd(B,9973) = 1 得知一定有解。
AC代码:
#in...
分类:
其他好文 时间:
2014-09-25 18:46:57
阅读次数:
194
听说这是数论中比较重要的部分了,一点点的总结吧。。一.线性同余方程与不定方程:单个一元线性方程求解方法:扩展欧几里得 exgcd模板:long long exgcd(long long a,long long b,long long &x,long long &y){ if(!b) { ...
分类:
其他好文 时间:
2014-09-22 20:09:43
阅读次数:
166
整数性质
时间限制:500 ms | 内存限制:65535 KB
难度:1
描述
我们知道,在数学中,对于任意两个正整数a和b,必定存在一对整数s、t使得sa+tb=gcd(a,b)。
输入多组测试数据。
每组数据输入两个非负整数a和b且a+b>0且a不等于b。
其中0
输出输出满足条件的 s 和 t 。
样例输入
2 4
3 8
737 6...
分类:
其他好文 时间:
2014-09-18 22:18:44
阅读次数:
165