coreBugZJ

此 blog 已棄。

EOJ 1056 線性同余方程

/*
2

EOJ 1056 線性同余方程
3

----問題描述：
6

形如ax≡b(mod m) 的方程，稱為線性同余方程。編寫程序求解線性同余方程（基于歐幾里德算法）。
8

----輸入：
11

測試包含多組測試數據.
13

每組測試數據只含一行,每行有三個整數a,b,m(0<a,b,m<1000000)
14

----輸出：
17

每組測試數據只輸出一行.如果有解,則按解的大小,從小到大輸出,兩兩之間用空格分開.如果沒有解,則輸出"No Answer."(不含引號).
19

----樣例輸入：
22

12 54 34
24

4 2 4
25

----樣例輸出：
28

13 30
30

No Answer.
31

----分析：
34

[1]
36

擴展歐幾里得算法：
37

a * x + b * y = d (1)
39

d = gcd( a, b )
40

由
42

b * x1 + (a % b) * y1 = d
44

==>
46

b * x1 + (a - (a / b) * b) * y1 = d
48

==>
50

b * x1 + a * y1 - (a / b) * b * y1 = d
52

==>
54

a * y1 + b * (x1 - (a / b) * y1) = d (2)
56

比較 (1) 和 (2) 得
59

x = y1
61

y = x1 - (a / b) * y1
62

[2]
66

方程 ax≡b(mod m)
67

==>
69

ax + my = b
71

==>
73

當且僅當 b % gcd(a, m) == 0 時有解
75

*/
78

#include <iostream>
81

#include <cstdio>
82

#include <algorithm>
83

using namespace std;
85

#define N 1000009
87

typedef __int64 Lint;
89

int gcd( int a, int b ) {
91

int t;
92

while ( 0 != b ) {
93

t = a;
94

a = b;
95

b = t % b;
96

}
97

return a;
98

}
99

100

int gcd_ex( int a, int b, int &x, int &y ) {
101

if ( b == 0 ) {
102

x = 1;
103

y = 0;
104

return a;
105

}
106

int d = gcd_ex( b, a % b, x, y );
107

int t = x;
108

x = y;
109

y = (int)(t - ( a / b ) * ((Lint)(y)));
110

return d;
111

}
112

113

int a, b, m;
114

115

int n, x[ N ];
116

117

int solve() {
118

int i, d, tx, ty;
119

Lint tmp;
120

d = gcd( a, m );
121

if ( b % d != 0 ) {
122

return 0;
123

}
124

gcd_ex( a / d, m / d, tx, ty );
125

tx *= b / d;
126

n = d;
127

for ( i = 0; i < n; ++i ) {
128

tmp = tx + m / d * ((Lint)(i));
129

if ( 0 > tmp ) {
130

x[ i ] = m - ( (-tmp) % m );
131

}
132

else {
133

x[ i ] = tmp % m;
134

}
135

}
136

sort( x, x + n );
137

return 1;
138

}
139

140

int main() {
141

int i;
142

while ( scanf( "%d%d%d", &a, &b, &m ) == 3 ) {
143

if ( solve() ) {
144

printf( "%d", x[ 0 ] );
145

for ( i = 1; i < n; ++i ) {
146

printf( " %d", x[ i ] );
147

}
148

printf( "\n" );
149

}
150

else {
151

puts( "No Answer." );
152

}
153

}
154

return 0;
155

}
156

posted on 2012-06-01 21:26 coreBugZJ 閱讀(798) 評論(0) 編輯收藏引用所屬分類: ACM 、Algorithm 、Mathematics 、課內作業

只有注冊用戶登錄后才能發表評論。
【推薦】100%開源！大型工業跨平臺軟件C++源碼提供，建模，組態！

相關文章: 微軟2014實習生及秋令營技術類職位在線測試 TopCoder SRM 593 DIV2 第三題生成全排列的非回溯方法（TopCoder SRM 591 DIV 2） A* 算法求解八數碼問題，POJ 1077 Eight POJ 1067 取石子游戲 POJ 2068 Nim POJ 2975 Nim POJ 3696 The Luckiest number POJ 3604 Professor Ben EOJ 1117 剩余定理

網站導航: 博客園 IT新聞 BlogJava 博問 Chat2DB 管理

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

coreBugZJ

My Links

Blog Stats

常用鏈接

留言簿(10)

隨筆分類(458)

隨筆檔案(268)

相冊

ACM

AI

LaTeX

安全

編程語言

好有道理

技術

開源

科學

數學

圖形圖像

文化

問題（練習＆有趣）

資源

最新隨筆

搜索

最新評論

閱讀排行榜

評論排行榜

EOJ 1056 線性同余方程