POJ-1006问题

maybe_madio

POJ-1006
$ U1 v8 C' |0 ]/ d* ^* o
+ }2 d$ Y' y* |5 ]在本地VC6.0下测试，得出正确结果， 为什么提交后总是Wrong Answer
0 x7 [: L& F( X  v+ _0 `
# d! G! G9 b: X  w现将代码贴出，请帮忙看看：

2. #include <iostream>, O2 m5 @0 t: ?
   
4. using namespace std;
   
5. 0 @! Z! I, _3 E
6. #define PH      0
   
7. # \' k- H* R( I% _6 s) o) [( ?
8. #define EM      1$ s9 F- J6 X2 y' C\\" X. ]
   
10. #define IN      2  ~# w! y$ D' {4 y
    
12. int get_Top(int days, int kind);
    
13. ; F6 W( K: m1 y! W, V  H: f
14. int triple(int ptop, int etop, int itop, int cur);
    
15. . s: n' |$ G/ y$ W4 f) c
16. bool is_Integer(float n);
    
17. : c  W/ U0 t5 _7 N# x6 S* M
18. int main()
    
19. ! F. s: z- W9 r: ?
20. {
    
21. \\" t- H! J6 ^1 m  ]! D8 X
22. int ptop,  etop,  itop;
    
23. # T, ?. k9 j) y) Q
24. int pdays, edays, idays, curdays;' ]( }( T8 U/ s$ G9 k\\" }
    
26. int count =0;
    
27. & T! G. ?1 R3 n2 M  U$ j
28. 
    
29. : f/ U0 ^+ N  ?! @  c4 V
30. while(cin >> pdays >> edays >> idays >> curdays)8 r, h, z6 M\\" _2 b
    
32. {/ }  a. W* c) D\\" e) o+ W
    
34.   if(pdays==edays&&pdays==idays&&pdays==curdays&&pdays==-1)
    
35. 3 J! J7 b: [7 z: z0 @$ i
36.    break;
    
37. 0 {# ?' o* s5 F5 n3 Z
38.   count++;
    
39. - R3 H* \4 }! |5 ?\\" ?5 G- e+ f
40.   ptop = get_Top(pdays,PH);8 c2 _' U% H: f3 @# j8 A
    
42.      etop = get_Top(edays,EM);% z\\" j7 f' @4 ]2 O. b$ r
    
44.      itop = get_Top(idays,IN);3 s' q6 H5 ]# s+ X) j\\" J
    
46.   cout<<"Case "<<count<<": "<<"the next triple peak occurs in "<<triple(ptop, etop, itop, curdays)<<" days.";4 X: r0 C4 A5 ~4 G
    
48. }
    
49. , q- W  h+ j! E$ q/ P3 n! Z
50. 
    
51. + O7 x3 y) ~+ [
52. return 0;7 j$ Z( a4 |3 M! {
    
54. }
    
55. . U- [7 N- v4 K- h\\" W% ^\\" R
56. int get_Top(int days, int kind)
    
57. \\" l1 z4 t) ]8 k- D( k7 ~% u4 O
58. {9 r, A$ S* L  K% E\\" e
    
60. switch(kind): [1 N6 P* m\\" A2 n3 R% n
    
62. {8 j# F. z  n5 a
    
64. case 0:
    
65. / i6 Q% R9 v! l$ e# Q$ b
66.   return days%23;
    
67. 6 @7 ^& F8 y) Q/ _& `5 E
68. case 1:8 T' Z5 a- v& J! ?/ U) z
    
70.   return days%28;\\" j3 R% [7 F$ D
    
72. case 2:  {; n9 }- P+ A1 t1 T  x
    
74.   return days%33;+ j; H; ]9 f# M6 W% H4 K: c/ q
    
76. }  R6 ?% H+ {. G\\" T
    
78. return 0;
    
79. ; C' b6 l8 B\\" t+ e$ `; c' c% S3 B
80. }! E\\" H: \1 x' [7 @) l  o
    
82. 
    
83. 2 m+ b, _5 s6 L& ~4 T- C
84. int triple(int ptop, int etop, int itop, int cur)9 i; D% Q( Y8 G5 v: E\\" V2 J
    
86. {0 T! R. S; P0 K/ R; C
    
88. float x, y, z;
    
89. ) H. X$ u. B- K1 a% I( e
90. int m1   = etop-ptop-5,! n\\" |; p9 M& C- K) Y. |
    
92.   m2   = itop-etop-5;. q# |/ O: C. f, ?7 q\\" R
    
94. int temp = 0;* Z: y2 @& m& O5 ]. }* u
    
96. for(z=2;  (itop+(z-1)*33) <= 21252;  z++)! g3 t& Y/ f\\" r: U
    
98. {6 B5 i; P\\" H4 p' b- G( Y
    
100.   temp = m2 + 33*z;0 ?, Y9 u8 I5 [
     
102.   y = (float)temp/28;
     
103. 8 _* G5 j* y' X% {
104.   x = (float)(m1+temp)/23;9 A. K, T1 ^# @
     
106.   if(x>0 && y>0 && is_Integer(x) && is_Integer(y))7 z' u1 c! [\\" F$ ~
     
108.    return itop + (z-1)*33 -cur;5 \4 L; v) _6 M/ O( l* g
     
110. }
     
111. + W' [, U4 L6 I5 _  x1 F' A7 ]) p
112. return 0;
     
113. - r5 ?2 |! j. ?( b- c
114. }
     
115. . ^& a0 z0 d9 w$ N3 F
116. 
     
117. 3 E$ I4 I\\" d6 C4 j+ ?\\" Q+ b% v7 `
118. bool is_Integer(float n)8 x. C\\" b, Q( G  {% r: q' r
     
120. {' E& o7 P9 X0 T, [/ u( R- l\\" Q# _; O% Q
     
122. if( (int)n - n ==0)
     
123. , a. Q9 z+ L+ u$ G8 T
124.   return true;
     
125. ' w) B2 f- ]* ?7 h  t\\" I: L* e
126. else
     
127. 5 @  h; d! l+ E5 Y3 O3 x' Q2 t
128.   return false;+ |0 a& g0 M1 {' {3 `2 X' I- v\\" C
     
130. }
     
131. $ G. o: @1 l5 ?+ j
132. 
     
133. 7 ?& e\\" W3 p$ u

#include <iostream><font class=\\"jammer\\">, O2 m5 @0 t: ?</font><br /> using namespace std;<br /> <span style=\\"display:none\\">0 @! Z! I, _3 E</span> #define PH      0<br /> <span style=\\"display:none\\"># \' k- H* R( I% _6 s) o) [( ?</span> #define EM      1<font class=\\"jammer\\">$ s9 F- J6 X2 y' C\\" X. ]</font><br /> #define IN      2<font class=\\"jammer\\">  ~# w! y$ D' {4 y</font><br /> int get_Top(int days, int kind);<br /> <span style=\\"display:none\\">; F6 W( K: m1 y! W, V  H: f</span> int triple(int ptop, int etop, int itop, int cur);<br /> <span style=\\"display:none\\">. s: n' |$ G/ y$ W4 f) c</span> bool is_Integer(float n);<br /> <span style=\\"display:none\\">: c  W/ U0 t5 _7 N# x6 S* M</span> int main()<br /> <span style=\\"display:none\\">! F. s: z- W9 r: ?</span> {<br /> <span style=\\"display:none\\">\\" t- H! J6 ^1 m  ]! D8 X</span> int ptop,  etop,  itop;<br /> <span style=\\"display:none\\"># T, ?. k9 j) y) Q</span> int pdays, edays, idays, curdays;<font class=\\"jammer\\">' ]( }( T8 U/ s$ G9 k\\" }</font><br /> int count =0;<br /> <span style=\\"display:none\\">& T! G. ?1 R3 n2 M  U$ j</span> <br /> <span style=\\"display:none\\">: f/ U0 ^+ N  ?! @  c4 V</span> while(cin >> pdays >> edays >> idays >> curdays)<font class=\\"jammer\\">8 r, h, z6 M\\" _2 b</font><br /> {<font class=\\"jammer\\">/ }  a. W* c) D\\" e) o+ W</font><br />   if(pdays==edays&&pdays==idays&&pdays==curdays&&pdays==-1)<br /> <span style=\\"display:none\\">3 J! J7 b: [7 z: z0 @$ i</span>    break;<br /> <span style=\\"display:none\\">0 {# ?' o* s5 F5 n3 Z</span>   count++;<br /> <span style=\\"display:none\\">- R3 H* \4 }! |5 ?\\" ?5 G- e+ f</span>   ptop = get_Top(pdays,PH);<font class=\\"jammer\\">8 c2 _' U% H: f3 @# j8 A</font><br />      etop = get_Top(edays,EM);<font class=\\"jammer\\">% z\\" j7 f' @4 ]2 O. b$ r</font><br />      itop = get_Top(idays,IN);<font class=\\"jammer\\">3 s' q6 H5 ]# s+ X) j\\" J</font><br />   cout<<"Case "<<count<<": "<<"the next triple peak occurs in "<<triple(ptop, etop, itop, curdays)<<" days.";<font class=\\"jammer\\">4 X: r0 C4 A5 ~4 G</font><br /> }<br /> <span style=\\"display:none\\">, q- W  h+ j! E$ q/ P3 n! Z</span> <br /> <span style=\\"display:none\\">+ O7 x3 y) ~+ [</span> return 0;<font class=\\"jammer\\">7 j$ Z( a4 |3 M! {</font><br /> }<br /> <span style=\\"display:none\\">. U- [7 N- v4 K- h\\" W% ^\\" R</span> int get_Top(int days, int kind)<br /> <span style=\\"display:none\\">\\" l1 z4 t) ]8 k- D( k7 ~% u4 O</span> {<font class=\\"jammer\\">9 r, A$ S* L  K% E\\" e</font><br /> switch(kind)<font class=\\"jammer\\">: [1 N6 P* m\\" A2 n3 R% n</font><br /> {<font class=\\"jammer\\">8 j# F. z  n5 a</font><br /> case 0:<br /> <span style=\\"display:none\\">/ i6 Q% R9 v! l$ e# Q$ b</span>   return days%23;<br /> <span style=\\"display:none\\">6 @7 ^& F8 y) Q/ _& `5 E</span> case 1:<font class=\\"jammer\\">8 T' Z5 a- v& J! ?/ U) z</font><br />   return days%28;<font class=\\"jammer\\">\\" j3 R% [7 F$ D</font><br /> case 2:<font class=\\"jammer\\">  {; n9 }- P+ A1 t1 T  x</font><br />   return days%33;<font class=\\"jammer\\">+ j; H; ]9 f# M6 W% H4 K: c/ q</font><br /> }<font class=\\"jammer\\">  R6 ?% H+ {. G\\" T</font><br /> return 0;<br /> <span style=\\"display:none\\">; C' b6 l8 B\\" t+ e$ `; c' c% S3 B</span> }<font class=\\"jammer\\">! E\\" H: \1 x' [7 @) l  o</font><br /> <br /> <span style=\\"display:none\\">2 m+ b, _5 s6 L& ~4 T- C</span> int triple(int ptop, int etop, int itop, int cur)<font class=\\"jammer\\">9 i; D% Q( Y8 G5 v: E\\" V2 J</font><br /> {<font class=\\"jammer\\">0 T! R. S; P0 K/ R; C</font><br /> float x, y, z;<br /> <span style=\\"display:none\\">) H. X$ u. B- K1 a% I( e</span> int m1   = etop-ptop-5,<font class=\\"jammer\\">! n\\" |; p9 M& C- K) Y. |</font><br />   m2   = itop-etop-5;<font class=\\"jammer\\">. q# |/ O: C. f, ?7 q\\" R</font><br /> int temp = 0;<font class=\\"jammer\\">* Z: y2 @& m& O5 ]. }* u</font><br /> for(z=2;  (itop+(z-1)*33) <= 21252;  z++)<font class=\\"jammer\\">! g3 t& Y/ f\\" r: U</font><br /> {<font class=\\"jammer\\">6 B5 i; P\\" H4 p' b- G( Y</font><br />   temp = m2 + 33*z;<font class=\\"jammer\\">0 ?, Y9 u8 I5 [</font><br />   y = (float)temp/28;<br /> <span style=\\"display:none\\">8 _* G5 j* y' X% {</span>   x = (float)(m1+temp)/23;<font class=\\"jammer\\">9 A. K, T1 ^# @</font><br />   if(x>0 && y>0 && is_Integer(x) && is_Integer(y))<font class=\\"jammer\\">7 z' u1 c! [\\" F$ ~</font><br />    return itop + (z-1)*33 -cur;<font class=\\"jammer\\">5 \4 L; v) _6 M/ O( l* g</font><br /> }<br /> <span style=\\"display:none\\">+ W' [, U4 L6 I5 _  x1 F' A7 ]) p</span> return 0;<br /> <span style=\\"display:none\\">- r5 ?2 |! j. ?( b- c</span> }<br /> <span style=\\"display:none\\">. ^& a0 z0 d9 w$ N3 F</span> <br /> <span style=\\"display:none\\">3 E$ I4 I\\" d6 C4 j+ ?\\" Q+ b% v7 `</span> bool is_Integer(float n)<font class=\\"jammer\\">8 x. C\\" b, Q( G  {% r: q' r</font><br /> {<font class=\\"jammer\\">' E& o7 P9 X0 T, [/ u( R- l\\" Q# _; O% Q</font><br /> if( (int)n - n ==0)<br /> <span style=\\"display:none\\">, a. Q9 z+ L+ u$ G8 T</span>   return true;<br /> <span style=\\"display:none\\">' w) B2 f- ]* ?7 h  t\\" I: L* e</span> else<br /> <span style=\\"display:none\\">5 @  h; d! l+ E5 Y3 O3 x' Q2 t</span>   return false;<font class=\\"jammer\\">+ |0 a& g0 M1 {' {3 `2 X' I- v\\" C</font><br /> }<br /> <span style=\\"display:none\\">$ G. o: @1 l5 ?+ j</span> <br /> <span style=\\"display:none\\">7 ?& e\\" W3 p$ u</span>

鸿光满面

好复杂啊         

39133120

好复杂啊，得仔细

maybe_madio

回复 ultra1989 的帖子
5 a  k: P! y* u3 c4 U/ r* G
$ l; }  l( k) R! @( S, V, o
    确实这道题使用剩余定理解效率太高了！ 哈哈！ 学习了！  我还没有看过数论方面的书！ 以后肯定要看看

; O8 W! W# T& O+ x1 L: C, H/ w, o两个程序一对比，一下就体现出数学的强大力量！

小草远在天涯

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

qbist

ultra1989

用float可能会有精度问题用中国剩余定理求解
参考程序：

#include<stdio.h>
int main(){
+ w6 E4 Y7 m6 }/ M8 d   int a,b,c,m,n,i,d=1;
   while(scanf("%d%d%d%d",&a,&b,&c,&m)){

if(a==-1&&b==-1&&c==-1&&m==-1) break;
. h, h8 T6 m* V2 q' e4 s& {3 F9 }! f
n=0;
1 i7 I7 I0 h& ~9 A7 h1 H. |- V
i=(5544*a+14421*b+1288*c-m+21252)%21252;
) x1 U9 K/ E8 t7 s" c; W* Q3 J% N/ C
' y6 D- i; P: S* O- \: r, _0 P( Q, _if(!i)

4 d3 O1 r4 ?$ g- ?i=21252;

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

! N% ~" B: W' c' Q( M   }

* t* B$ Y8 w; Z% x$ Kreturn 0;
}
$ z7 S& u) b1 Q6 C
' D; D1 v, A0 H( p* c) ]