POJ-1006问题

maybe_madio

POJ-1006
; A  |' ?0 j( @: E
在本地VC6.0下测试，得出正确结果， 为什么提交后总是Wrong Answer
; ?  t6 _5 u. l/ {6 m
现将代码贴出，请帮忙看看：

2. #include <iostream># S) Z1 H4 l& t* y/ a\\" e: ]$ F
   
4. using namespace std;) Z4 P3 p7 j\\" d  {; ?
   
6. #define PH      05 v9 w5 Z$ l! `) X  v2 ^
   
8. #define EM      12 P- D* b1 Q' S* X
   
10. #define IN      2
    
11. - o$ p- t# Y$ }. g
12. int get_Top(int days, int kind);
    
13. ( b1 e7 ~\\" \! S4 B
14. int triple(int ptop, int etop, int itop, int cur);# C4 ^( q/ ~8 N2 j! M/ [
    
16. bool is_Integer(float n);
    
17. : R! B5 E; ]  }/ [
18. int main()
    
19. * n, T% T  h) y
20. {
    
21. 6 S+ J# k; g1 b* i/ h; c' g
22. int ptop,  etop,  itop;
    
23. # J* h3 h2 o) U2 Z! I5 S
24. int pdays, edays, idays, curdays;( l3 a9 |/ T# Q! x* a# z: m
    
26. int count =0;/ c9 I5 W% K  l* j1 H+ E) D! d# n
    
28. 
    
29. $ M% \( L: ~* I7 h, a- N2 M! u
30. while(cin >> pdays >> edays >> idays >> curdays)' x2 ]' {6 z) @6 p; i+ G
    
32. {
    
33. + `  X* E0 e' E5 {* M, r
34.   if(pdays==edays&&pdays==idays&&pdays==curdays&&pdays==-1)
    
35. - V2 o' u, n( _
36.    break;
    
37. / y6 o* p: M( P6 K6 U
38.   count++;; l1 o5 |( ?5 t4 |5 G
    
40.   ptop = get_Top(pdays,PH);' d$ q# Z. }4 f6 z2 F7 p
    
42.      etop = get_Top(edays,EM);; {1 [2 _\\" E3 d+ }% `7 w% }: m: X2 I: D
    
44.      itop = get_Top(idays,IN);- k0 B/ [) W0 P3 ?2 B7 R
    
46.   cout<<"Case "<<count<<": "<<"the next triple peak occurs in "<<triple(ptop, etop, itop, curdays)<<" days.";6 ?1 p) Q4 N  V- G+ [
    
48. }  p7 M\\" V( m( u5 i/ g$ p3 d
    
50. - w, A# m( M( V( ~, `# _! b
    
52. return 0;\\" ]7 [4 t) ?6 P0 N$ I
    
54. }
    
55. ( @$ c0 }% C3 U/ e
56. int get_Top(int days, int kind)
    
57. 3 ^; I6 T  I/ ~5 t1 Y
58. {  F- J! R' q  k; ~8 g6 k
    
60. switch(kind)1 a  Z0 W* d  R4 x+ t. t! J
    
62. {7 ~( u( ~. h% t1 Y7 N/ w) a' t+ F
    
64. case 0:
    
65. 3 [1 g5 t7 K3 f% h
66.   return days%23;6 E8 a, j: ]  E: U, j' T3 g. B
    
68. case 1:! n5 Y1 p/ F8 v
    
70.   return days%28;\\" K$ [+ h, L: O7 B
    
72. case 2:  ?5 d, F- G2 P  d& V+ k0 S2 \
    
74.   return days%33;
    
75. 4 |) O0 Y+ Z) T5 L* `
76. }5 g, @3 }* V( ~/ n4 F9 g- ]5 H) [' ]8 B
    
78. return 0;
    
79. % N, {9 C8 J+ q8 y
80. }# Q5 n; i7 `$ F7 U, U5 P7 k4 v2 b
    
82. 5 Z& ]. _\\" m7 g5 L/ C
    
84. int triple(int ptop, int etop, int itop, int cur)8 f. V% y/ R' t
    
86. {
    
87. $ d& Z* v9 e! o
88. float x, y, z;
    
89. \\" V1 }' R& y2 H# f) v
90. int m1   = etop-ptop-5,
    
91. + Q$ {+ M$ {2 y  L4 T5 j4 \, X1 _
92.   m2   = itop-etop-5;\\" X0 [3 f& d5 h( h* s
    
94. int temp = 0;
    
95. 6 |! R( ~, M1 E$ c6 g) f9 S  G
96. for(z=2;  (itop+(z-1)*33) <= 21252;  z++)9 c. `3 o; B, P6 }7 W  F; r
    
98. {3 s/ Q\\" k1 x6 d5 p: t6 h
    
100.   temp = m2 + 33*z;
     
101. ) q\\" T- ?1 q; T  [/ Q\\" j
102.   y = (float)temp/28;
     
103. ) D) Y7 n. i3 U7 m# P
104.   x = (float)(m1+temp)/23;
     
105. 2 C. a. `8 f# J- m8 K) V) d# _
106.   if(x>0 && y>0 && is_Integer(x) && is_Integer(y))0 ?, D5 p/ F2 c2 Y5 J
     
108.    return itop + (z-1)*33 -cur;5 S' y- ?0 p1 A7 e- H6 }$ B- v
     
110. }) |$ v/ X( m6 ^5 K' s
     
112. return 0;# [0 {5 b$ J# C\\" n4 x: \/ i
     
114. }& S3 j) r$ N) [7 ]. H' y
     
116. 
     
117. $ B' N- A* q2 D
118. bool is_Integer(float n)
     
119. 6 r/ X, n; D+ @( g5 K/ n
120. {
     
121. # z: X, D( U' x5 _$ c/ z
122. if( (int)n - n ==0); C7 G! G4 B( g) ~3 b+ n$ J+ d9 O
     
124.   return true;
     
125. 9 W' c; s2 F. Z\\" d2 q* C- ]\\" V1 I
126. else2 x4 V! B* j- ^6 n) {4 f: x4 O4 M
     
128.   return false;
     
129. 2 S& D% V- \! t7 M6 ^
130. }
     
131. - P$ j- Q# ]* }\\" V
132. 3 e. k4 M8 A/ v1 H, x% ~7 ?5 [
     

#include <iostream><font class=\\"jammer\\"># S) Z1 H4 l& t* y/ a\\" e: ]$ F</font><br /> using namespace std;<font class=\\"jammer\\">) Z4 P3 p7 j\\" d  {; ?</font><br /> #define PH      0<font class=\\"jammer\\">5 v9 w5 Z$ l! `) X  v2 ^</font><br /> #define EM      1<font class=\\"jammer\\">2 P- D* b1 Q' S* X</font><br /> #define IN      2<br /> <span style=\\"display:none\\">- o$ p- t# Y$ }. g</span> int get_Top(int days, int kind);<br /> <span style=\\"display:none\\">( b1 e7 ~\\" \! S4 B</span> int triple(int ptop, int etop, int itop, int cur);<font class=\\"jammer\\"># C4 ^( q/ ~8 N2 j! M/ [</font><br /> bool is_Integer(float n);<br /> <span style=\\"display:none\\">: R! B5 E; ]  }/ [</span> int main()<br /> <span style=\\"display:none\\">* n, T% T  h) y</span> {<br /> <span style=\\"display:none\\">6 S+ J# k; g1 b* i/ h; c' g</span> int ptop,  etop,  itop;<br /> <span style=\\"display:none\\"># J* h3 h2 o) U2 Z! I5 S</span> int pdays, edays, idays, curdays;<font class=\\"jammer\\">( l3 a9 |/ T# Q! x* a# z: m</font><br /> int count =0;<font class=\\"jammer\\">/ c9 I5 W% K  l* j1 H+ E) D! d# n</font><br /> <br /> <span style=\\"display:none\\">$ M% \( L: ~* I7 h, a- N2 M! u</span> while(cin >> pdays >> edays >> idays >> curdays)<font class=\\"jammer\\">' x2 ]' {6 z) @6 p; i+ G</font><br /> {<br /> <span style=\\"display:none\\">+ `  X* E0 e' E5 {* M, r</span>   if(pdays==edays&&pdays==idays&&pdays==curdays&&pdays==-1)<br /> <span style=\\"display:none\\">- V2 o' u, n( _</span>    break;<br /> <span style=\\"display:none\\">/ y6 o* p: M( P6 K6 U</span>   count++;<font class=\\"jammer\\">; l1 o5 |( ?5 t4 |5 G</font><br />   ptop = get_Top(pdays,PH);<font class=\\"jammer\\">' d$ q# Z. }4 f6 z2 F7 p</font><br />      etop = get_Top(edays,EM);<font class=\\"jammer\\">; {1 [2 _\\" E3 d+ }% `7 w% }: m: X2 I: D</font><br />      itop = get_Top(idays,IN);<font class=\\"jammer\\">- k0 B/ [) W0 P3 ?2 B7 R</font><br />   cout<<"Case "<<count<<": "<<"the next triple peak occurs in "<<triple(ptop, etop, itop, curdays)<<" days.";<font class=\\"jammer\\">6 ?1 p) Q4 N  V- G+ [</font><br /> }<font class=\\"jammer\\">  p7 M\\" V( m( u5 i/ g$ p3 d</font><br /> <font class=\\"jammer\\">- w, A# m( M( V( ~, `# _! b</font><br /> return 0;<font class=\\"jammer\\">\\" ]7 [4 t) ?6 P0 N$ I</font><br /> }<br /> <span style=\\"display:none\\">( @$ c0 }% C3 U/ e</span> int get_Top(int days, int kind)<br /> <span style=\\"display:none\\">3 ^; I6 T  I/ ~5 t1 Y</span> {<font class=\\"jammer\\">  F- J! R' q  k; ~8 g6 k</font><br /> switch(kind)<font class=\\"jammer\\">1 a  Z0 W* d  R4 x+ t. t! J</font><br /> {<font class=\\"jammer\\">7 ~( u( ~. h% t1 Y7 N/ w) a' t+ F</font><br /> case 0:<br /> <span style=\\"display:none\\">3 [1 g5 t7 K3 f% h</span>   return days%23;<font class=\\"jammer\\">6 E8 a, j: ]  E: U, j' T3 g. B</font><br /> case 1:<font class=\\"jammer\\">! n5 Y1 p/ F8 v</font><br />   return days%28;<font class=\\"jammer\\">\\" K$ [+ h, L: O7 B</font><br /> case 2:<font class=\\"jammer\\">  ?5 d, F- G2 P  d& V+ k0 S2 \</font><br />   return days%33;<br /> <span style=\\"display:none\\">4 |) O0 Y+ Z) T5 L* `</span> }<font class=\\"jammer\\">5 g, @3 }* V( ~/ n4 F9 g- ]5 H) [' ]8 B</font><br /> return 0;<br /> <span style=\\"display:none\\">% N, {9 C8 J+ q8 y</span> }<font class=\\"jammer\\"># Q5 n; i7 `$ F7 U, U5 P7 k4 v2 b</font><br /> <font class=\\"jammer\\">5 Z& ]. _\\" m7 g5 L/ C</font><br /> int triple(int ptop, int etop, int itop, int cur)<font class=\\"jammer\\">8 f. V% y/ R' t</font><br /> {<br /> <span style=\\"display:none\\">$ d& Z* v9 e! o</span> float x, y, z;<br /> <span style=\\"display:none\\">\\" V1 }' R& y2 H# f) v</span> int m1   = etop-ptop-5,<br /> <span style=\\"display:none\\">+ Q$ {+ M$ {2 y  L4 T5 j4 \, X1 _</span>   m2   = itop-etop-5;<font class=\\"jammer\\">\\" X0 [3 f& d5 h( h* s</font><br /> int temp = 0;<br /> <span style=\\"display:none\\">6 |! R( ~, M1 E$ c6 g) f9 S  G</span> for(z=2;  (itop+(z-1)*33) <= 21252;  z++)<font class=\\"jammer\\">9 c. `3 o; B, P6 }7 W  F; r</font><br /> {<font class=\\"jammer\\">3 s/ Q\\" k1 x6 d5 p: t6 h</font><br />   temp = m2 + 33*z;<br /> <span style=\\"display:none\\">) q\\" T- ?1 q; T  [/ Q\\" j</span>   y = (float)temp/28;<br /> <span style=\\"display:none\\">) D) Y7 n. i3 U7 m# P</span>   x = (float)(m1+temp)/23;<br /> <span style=\\"display:none\\">2 C. a. `8 f# J- m8 K) V) d# _</span>   if(x>0 && y>0 && is_Integer(x) && is_Integer(y))<font class=\\"jammer\\">0 ?, D5 p/ F2 c2 Y5 J</font><br />    return itop + (z-1)*33 -cur;<font class=\\"jammer\\">5 S' y- ?0 p1 A7 e- H6 }$ B- v</font><br /> }<font class=\\"jammer\\">) |$ v/ X( m6 ^5 K' s</font><br /> return 0;<font class=\\"jammer\\"># [0 {5 b$ J# C\\" n4 x: \/ i</font><br /> }<font class=\\"jammer\\">& S3 j) r$ N) [7 ]. H' y</font><br /> <br /> <span style=\\"display:none\\">$ B' N- A* q2 D</span> bool is_Integer(float n)<br /> <span style=\\"display:none\\">6 r/ X, n; D+ @( g5 K/ n</span> {<br /> <span style=\\"display:none\\"># z: X, D( U' x5 _$ c/ z</span> if( (int)n - n ==0)<font class=\\"jammer\\">; C7 G! G4 B( g) ~3 b+ n$ J+ d9 O</font><br />   return true;<br /> <span style=\\"display:none\\">9 W' c; s2 F. Z\\" d2 q* C- ]\\" V1 I</span> else<font class=\\"jammer\\">2 x4 V! B* j- ^6 n) {4 f: x4 O4 M</font><br />   return false;<br /> <span style=\\"display:none\\">2 S& D% V- \! t7 M6 ^</span> }<br /> <span style=\\"display:none\\">- P$ j- Q# ]* }\\" V</span> <font class=\\"jammer\\">3 e. k4 M8 A/ v1 H, x% ~7 ?5 [</font><br />

  ]8 S* [. T" y- e: d' q

鸿光满面

好复杂啊         

39133120

好复杂啊，得仔细

maybe_madio

回复 ultra1989 的帖子
5 ?* F3 S* h* h- e: U. i
1 \$ ^3 V: g: c+ _" O8 }( ?$ w
    确实这道题使用剩余定理解效率太高了！ 哈哈！ 学习了！  我还没有看过数论方面的书！ 以后肯定要看看

( l) R" N4 s+ m% m1 K) Q两个程序一对比，一下就体现出数学的强大力量！

小草远在天涯

好复杂啊，我是初级的。。。路过。。。

qbist

ultra1989

) _, y1 N) U! K$ ]
用float可能会有精度问题用中国剩余定理求解
! w* P- M2 ]1 v( v% E2 ?2 O参考程序：

. r+ j0 R. b7 {#include<stdio.h>
int main(){
7 X4 q: L5 J' W) t. F% L   int a,b,c,m,n,i,d=1;
! H4 j/ H* I) b   while(scanf("%d%d%d%d",&a,&b,&c,&m)){
( |$ f8 q) C. P  _
if(a==-1&&b==-1&&c==-1&&m==-1) break;

5 X/ p' N/ d6 T: P2 Bn=0;
5 o  v1 O  m/ H' m7 L4 d3 h1 S6 S, a
i=(5544*a+14421*b+1288*c-m+21252)%21252;
+ W) j6 t$ e9 G* ^- c2 j8 `
$ b+ K4 d8 g: z4 u$ [) H. Pif(!i)
! x; M# N6 d1 `. v$ X
i=21252;

printf("Case %d: the next triple peak occurs in %d days.\n",d++,i);

$ a4 J( d# Y9 y  }3 j* y   }

return 0;
}
5 X4 ]0 g9 e) a$ C