POJ-1006问题

maybe_madio

POJ-1006

/ O+ l0 e/ p" g$ R2 Q在本地VC6.0下测试，得出正确结果， 为什么提交后总是Wrong Answer

现将代码贴出，请帮忙看看：

2. #include <iostream>; u. d9 D$ m6 c) c$ t
   
4. using namespace std;. r& H2 E: K5 J( s
   
6. #define PH      0
   
7. 2 \\\" h+ u! V7 W* [
8. #define EM      1& H( L7 C4 C# r+ Y* v
   
10. #define IN      2
    
11. 0 [6 G$ h/ w  w1 p& o8 t\\" q0 ?
12. int get_Top(int days, int kind);4 v+ S4 R. s& b5 H, v# i, ~: v. M
    
14. int triple(int ptop, int etop, int itop, int cur);0 N: m) h* u5 h
    
16. bool is_Integer(float n);2 C! O7 q+ l- P8 H7 v  b% `9 k
    
18. int main(); Z1 ]5 U$ k/ n- f, Y% X& Q
    
20. {9 z0 {+ G# L: Q
    
22. int ptop,  etop,  itop;
    
23. ' `- E( `/ `8 I4 z
24. int pdays, edays, idays, curdays;
    
25. 7 D0 O( |/ }% t& Q
26. int count =0;$ i2 R  l) q4 Z  D
    
28. $ U3 z+ q\\" g1 n
    
30. while(cin >> pdays >> edays >> idays >> curdays)! @+ y0 _  h2 Z- q
    
32. {) e0 G' r6 R  Y: K) O
    
34.   if(pdays==edays&&pdays==idays&&pdays==curdays&&pdays==-1), R- n9 o( p\\" E# n& n
    
36.    break;
    
37. 8 h0 M6 @- u+ A( c4 ~8 L8 F
38.   count++;! g4 W4 d; o9 [) W  q+ h7 l5 \
    
40.   ptop = get_Top(pdays,PH);# \' p. r- j+ S/ P# w8 {\\" ~: P
    
42.      etop = get_Top(edays,EM);6 j1 P, r- l+ s# _! P, m$ J, W7 T
    
44.      itop = get_Top(idays,IN);# ?; Q2 t% p& o# ^2 h8 U6 V
    
46.   cout<<"Case "<<count<<": "<<"the next triple peak occurs in "<<triple(ptop, etop, itop, curdays)<<" days.";
    
47. ( ?! i( l' ^# j7 B) j' O
48. }
    
49. 5 Y  k' F\\" K  s
50. & i& g# S2 `+ K) u0 ^
    
52. return 0;
    
53. $ ~  ]\\" p3 T; X& T
54. }% S( E! I  T& V; r0 Q
    
56. int get_Top(int days, int kind)
    
57. : g, @, y( P% {. ?6 j! o. @& M\\" B
58. {3 o2 t7 f2 z: E2 F% q/ T' W! t0 o4 ~* r
    
60. switch(kind)9 O; V4 P4 i+ e- E, j0 K1 T
    
62. {  i. A* H3 r8 i3 G\\" L7 H% g' P
    
64. case 0:+ h( r0 ^) j# D0 N3 J& @
    
66.   return days%23;
    
67. / r& _) h7 |& H, p; p; Y+ X
68. case 1:
    
69. ! u+ @  D/ L% `; r
70.   return days%28;4 d/ D' V\\" k5 v$ N) |
    
72. case 2:
    
73. 9 c/ z- V* ?4 ~4 M6 r; X\\" H# D2 g
74.   return days%33;
    
75. ' o( Q3 m2 D; b& _1 i
76. }  p- c6 {1 {; y! P
    
78. return 0;
    
79. : v1 ]4 e0 A$ i: y- a  p* s
80. }9 H6 X) \7 f\\" S) o9 [
    
82. 4 k7 z/ N; }$ l6 j  r6 T
    
84. int triple(int ptop, int etop, int itop, int cur)/ A( A( t2 K/ `( o$ @& n! c! T( G
    
86. {1 a9 M9 Q( d  e9 i2 h3 O# R
    
88. float x, y, z;
    
89. & N( i  g, R5 i7 C
90. int m1   = etop-ptop-5,
    
91.   r  r3 y$ y1 P+ W1 [; Z
92.   m2   = itop-etop-5;9 ~  _0 U1 ^0 U1 f
    
94. int temp = 0;
    
95. 4 c% E4 N2 ~/ U* Z' j7 D
96. for(z=2;  (itop+(z-1)*33) <= 21252;  z++)  w2 J\\" H4 l+ h; X! P8 z) N
    
98. {
    
99. & f! F4 R' X/ u
100.   temp = m2 + 33*z;
     
101. 2 O; @! n/ M* l/ f, Q
102.   y = (float)temp/28;) d# k' N+ {. K; |% s# E# B% b
     
104.   x = (float)(m1+temp)/23;
     
105. \\" ^8 V. Q. z  I4 z( [4 u
106.   if(x>0 && y>0 && is_Integer(x) && is_Integer(y))
     
107. \\" ?, ]6 ^( }( ~
108.    return itop + (z-1)*33 -cur;& g! f4 |6 N: v
     
110. }) @7 s9 u. Y: C. O( r3 x( p0 o. U
     
112. return 0;: j' o% R9 _- W# z: ]: i- o! t
     
114. }
     
115. ) M# G0 j, Q; O5 O0 R5 l8 Q- c% T
116. 
     
117. - ]2 r9 v7 E+ s, t9 W
118. bool is_Integer(float n)
     
119. # [6 c- h6 a* R3 I
120. {& n6 b\\" k. P5 J' @% t; M! @
     
122. if( (int)n - n ==0)* q9 n% d* n+ `\\" y/ q
     
124.   return true;, J8 T- @/ e. a6 b  u  ~
     
126. else
     
127. 5 ]$ e\\" x& o/ t* [  B! N8 }, t+ }
128.   return false;
     
129. 0 ]4 ^+ }5 }8 R; }- N
130. }) Q2 c# d% K4 H2 D
     
132. 
     
133. 9 l; j' j; d* v/ h

#include <iostream><font class=\\"jammer\\">; u. d9 D$ m6 c) c$ t</font><br /> using namespace std;<font class=\\"jammer\\">. r& H2 E: K5 J( s</font><br /> #define PH      0<br /> <span style=\\"display:none\\">2 \\\" h+ u! V7 W* [</span> #define EM      1<font class=\\"jammer\\">& H( L7 C4 C# r+ Y* v</font><br /> #define IN      2<br /> <span style=\\"display:none\\">0 [6 G$ h/ w  w1 p& o8 t\\" q0 ?</span> int get_Top(int days, int kind);<font class=\\"jammer\\">4 v+ S4 R. s& b5 H, v# i, ~: v. M</font><br /> int triple(int ptop, int etop, int itop, int cur);<font class=\\"jammer\\">0 N: m) h* u5 h</font><br /> bool is_Integer(float n);<font class=\\"jammer\\">2 C! O7 q+ l- P8 H7 v  b% `9 k</font><br /> int main()<font class=\\"jammer\\">; Z1 ]5 U$ k/ n- f, Y% X& Q</font><br /> {<font class=\\"jammer\\">9 z0 {+ G# L: Q</font><br /> int ptop,  etop,  itop;<br /> <span style=\\"display:none\\">' `- E( `/ `8 I4 z</span> int pdays, edays, idays, curdays;<br /> <span style=\\"display:none\\">7 D0 O( |/ }% t& Q</span> int count =0;<font class=\\"jammer\\">$ i2 R  l) q4 Z  D</font><br /> <font class=\\"jammer\\">$ U3 z+ q\\" g1 n</font><br /> while(cin >> pdays >> edays >> idays >> curdays)<font class=\\"jammer\\">! @+ y0 _  h2 Z- q</font><br /> {<font class=\\"jammer\\">) e0 G' r6 R  Y: K) O</font><br />   if(pdays==edays&&pdays==idays&&pdays==curdays&&pdays==-1)<font class=\\"jammer\\">, R- n9 o( p\\" E# n& n</font><br />    break;<br /> <span style=\\"display:none\\">8 h0 M6 @- u+ A( c4 ~8 L8 F</span>   count++;<font class=\\"jammer\\">! g4 W4 d; o9 [) W  q+ h7 l5 \</font><br />   ptop = get_Top(pdays,PH);<font class=\\"jammer\\"># \' p. r- j+ S/ P# w8 {\\" ~: P</font><br />      etop = get_Top(edays,EM);<font class=\\"jammer\\">6 j1 P, r- l+ s# _! P, m$ J, W7 T</font><br />      itop = get_Top(idays,IN);<font class=\\"jammer\\"># ?; Q2 t% p& o# ^2 h8 U6 V</font><br />   cout<<"Case "<<count<<": "<<"the next triple peak occurs in "<<triple(ptop, etop, itop, curdays)<<" days.";<br /> <span style=\\"display:none\\">( ?! i( l' ^# j7 B) j' O</span> }<br /> <span style=\\"display:none\\">5 Y  k' F\\" K  s</span> <font class=\\"jammer\\">& i& g# S2 `+ K) u0 ^</font><br /> return 0;<br /> <span style=\\"display:none\\">$ ~  ]\\" p3 T; X& T</span> }<font class=\\"jammer\\">% S( E! I  T& V; r0 Q</font><br /> int get_Top(int days, int kind)<br /> <span style=\\"display:none\\">: g, @, y( P% {. ?6 j! o. @& M\\" B</span> {<font class=\\"jammer\\">3 o2 t7 f2 z: E2 F% q/ T' W! t0 o4 ~* r</font><br /> switch(kind)<font class=\\"jammer\\">9 O; V4 P4 i+ e- E, j0 K1 T</font><br /> {<font class=\\"jammer\\">  i. A* H3 r8 i3 G\\" L7 H% g' P</font><br /> case 0:<font class=\\"jammer\\">+ h( r0 ^) j# D0 N3 J& @</font><br />   return days%23;<br /> <span style=\\"display:none\\">/ r& _) h7 |& H, p; p; Y+ X</span> case 1:<br /> <span style=\\"display:none\\">! u+ @  D/ L% `; r</span>   return days%28;<font class=\\"jammer\\">4 d/ D' V\\" k5 v$ N) |</font><br /> case 2:<br /> <span style=\\"display:none\\">9 c/ z- V* ?4 ~4 M6 r; X\\" H# D2 g</span>   return days%33;<br /> <span style=\\"display:none\\">' o( Q3 m2 D; b& _1 i</span> }<font class=\\"jammer\\">  p- c6 {1 {; y! P</font><br /> return 0;<br /> <span style=\\"display:none\\">: v1 ]4 e0 A$ i: y- a  p* s</span> }<font class=\\"jammer\\">9 H6 X) \7 f\\" S) o9 [</font><br /> <font class=\\"jammer\\">4 k7 z/ N; }$ l6 j  r6 T</font><br /> int triple(int ptop, int etop, int itop, int cur)<font class=\\"jammer\\">/ A( A( t2 K/ `( o$ @& n! c! T( G</font><br /> {<font class=\\"jammer\\">1 a9 M9 Q( d  e9 i2 h3 O# R</font><br /> float x, y, z;<br /> <span style=\\"display:none\\">& N( i  g, R5 i7 C</span> int m1   = etop-ptop-5,<br /> <span style=\\"display:none\\">  r  r3 y$ y1 P+ W1 [; Z</span>   m2   = itop-etop-5;<font class=\\"jammer\\">9 ~  _0 U1 ^0 U1 f</font><br /> int temp = 0;<br /> <span style=\\"display:none\\">4 c% E4 N2 ~/ U* Z' j7 D</span> for(z=2;  (itop+(z-1)*33) <= 21252;  z++)<font class=\\"jammer\\">  w2 J\\" H4 l+ h; X! P8 z) N</font><br /> {<br /> <span style=\\"display:none\\">& f! F4 R' X/ u</span>   temp = m2 + 33*z;<br /> <span style=\\"display:none\\">2 O; @! n/ M* l/ f, Q</span>   y = (float)temp/28;<font class=\\"jammer\\">) d# k' N+ {. K; |% s# E# B% b</font><br />   x = (float)(m1+temp)/23;<br /> <span style=\\"display:none\\">\\" ^8 V. Q. z  I4 z( [4 u</span>   if(x>0 && y>0 && is_Integer(x) && is_Integer(y))<br /> <span style=\\"display:none\\">\\" ?, ]6 ^( }( ~</span>    return itop + (z-1)*33 -cur;<font class=\\"jammer\\">& g! f4 |6 N: v</font><br /> }<font class=\\"jammer\\">) @7 s9 u. Y: C. O( r3 x( p0 o. U</font><br /> return 0;<font class=\\"jammer\\">: j' o% R9 _- W# z: ]: i- o! t</font><br /> }<br /> <span style=\\"display:none\\">) M# G0 j, Q; O5 O0 R5 l8 Q- c% T</span> <br /> <span style=\\"display:none\\">- ]2 r9 v7 E+ s, t9 W</span> bool is_Integer(float n)<br /> <span style=\\"display:none\\"># [6 c- h6 a* R3 I</span> {<font class=\\"jammer\\">& n6 b\\" k. P5 J' @% t; M! @</font><br /> if( (int)n - n ==0)<font class=\\"jammer\\">* q9 n% d* n+ `\\" y/ q</font><br />   return true;<font class=\\"jammer\\">, J8 T- @/ e. a6 b  u  ~</font><br /> else<br /> <span style=\\"display:none\\">5 ]$ e\\" x& o/ t* [  B! N8 }, t+ }</span>   return false;<br /> <span style=\\"display:none\\">0 ]4 ^+ }5 }8 R; }- N</span> }<font class=\\"jammer\\">) Q2 c# d% K4 H2 D</font><br /> <br /> <span style=\\"display:none\\">9 l; j' j; d* v/ h</span>

/ {. j3 E$ n0 b+ u+ b( W- D: i

ultra1989

( W" x! Q; T+ T# u" l
用float可能会有精度问题用中国剩余定理求解
$ i8 Y& c' [' X1 P& S. h参考程序：

#include<stdio.h>
int main(){
5 \1 t$ h3 @% \" p' {% V, g   int a,b,c,m,n,i,d=1;
   while(scanf("%d%d%d%d",&a,&b,&c,&m)){
1 E5 @" p1 W+ B' N5 t
if(a==-1&&b==-1&&c==-1&&m==-1) break;
+ M# x8 O5 A( D3 t& J  e8 D+ {
' C6 x9 K1 ^2 c" un=0;
" M. B6 d) P6 _$ Y
i=(5544*a+14421*b+1288*c-m+21252)%21252;
* ?% B3 I, y0 z/ l/ t* t
& n; n/ x+ f4 U7 A' x' tif(!i)

i=21252;
$ c" b' w  Z; @* K
, [+ m: k2 M, q# ^- d4 Hprintf("Case %d: the next triple peak occurs in %d days.\n",d++,i);
$ _' x8 r, n' d. ~2 s& Z: J
3 Z# s' F! _. w2 K3 V9 N   }
; ?  Q" o: ^6 ?) K
return 0;
, c4 z9 @3 P: s. r}

+ V9 p# g" W) Z! e0 y

qbist

小草远在天涯

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

maybe_madio

回复 ultra1989 的帖子

& h. I5 i5 r* i& x6 L
    确实这道题使用剩余定理解效率太高了！ 哈哈！ 学习了！  我还没有看过数论方面的书！ 以后肯定要看看

两个程序一对比，一下就体现出数学的强大力量！

39133120

好复杂啊，得仔细

鸿光满面

好复杂啊