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

隨筆 - 97, 文章 - 22, 評論 - 81, 引用 - 0
數據加載中……

SGU 110 - 119 解題報告

110 Dungeon                                                計算幾何:射線和球體相交
111 Very simple problem                             二分枚舉
112 a^b-b^a                                                二分求冪 
113 Nearly prime numbers                          數論:素數篩選
114 Telecasting station                               枚舉
115 Calendar                                                簡單模擬
116 Index of super-prime                            動態規劃:記憶化搜索
117 Counting                                                二分求冪
118 Digital Root                                           數學題
119 Magic Pairs                                            數論:擴展歐幾里得


110 Dungeon

      計算幾何:射線和球體相交

       題意:給定一條3D射線 N個球體,問射線的各種反彈經過的球體編號,如果反射超過10次,只需要輸出前十次。

       題解:很明顯,首先需要枚舉每個球體和當前射線的相交情況,如果和多個球體相交肯定是取距離最近的球體進行相交計算,然后計算出空間反射射線,重復以上操作11次,如果某次的反射射線已經不能和任何球體相交那么直接跳出這個過程。

 


圖1
       1) 射線和球體相交

 

       如圖1,表示射線AB和球O相交于B點,反射后的射線為BCBD為球心O到交點的延長線,那么必定滿足以下幾個條件:

       已知A點坐標(xa, ya, za),單位向量tab

       B點坐標 (xb, yb, zb)  =  (xa, ya, za) + |AB| * tab   (1)

         B點在球O上,所以B點坐標滿足

       | (xb, yb, zb)  - (xo, yo, zo) | = OB = Ro              (2)

       (1)代入(2),可計算出|AB|的長度(由于是射線和球體求交,這兩個方程求出來的是直線和球體的交點,所以還需要計算方向判斷交點的可行性,如果交點有兩個,則取距離小的那個,如果只有一個交點,說明是相切),然后利用(1)式計算出B點坐標。

       2) 射線經過球體的反射射線

       假設C點在A經過B點反射后的反射射線上,并且AB = BC,那么AC必定和OB延長線相交,假設交于D點;

       單位向量cos(ABD) = Iba ­­* I­­bd

BD = AB * cos(ABD)

向量AD = 向量AB + 向量BD

向量AC = 2 * AD

C (x­c, yc, zc) = 向量AC + (xa, ya, za)

向量BC求出后返回1) 繼續判斷和別的球的相交情況;

111 Very simple problem

       二分答案 + 大數模擬

       題意:給定X( X <= 10100),求X開方后的下取整。

       題解:數組模擬大數。二分答案A,找滿足 A*A <= X 最大的A即為答案。

 
112 
a^b-b^a

       二分求冪 + 大數模擬

       題意:求 ab - ba ( 0 < a, b < 100)

       題解:數組模擬大數。二分求 AB

       B == 0  則 AB = 1; 否則 AB = (A2)(B/2) * ( (B mod 2) ? A : 1 ); 遞歸求解。

 

113 Nearly prime numbers

       素數篩選 + 素數判定

       題意:判斷一個數是否為兩個素數的積。

       題解:利用傳統的素數篩選將 1-31623 的素數篩選出來,31623為平方大于等于109的最小整數,因為一個數的開方內找不到一個(非1或它本身)因子的話它本身就是素數,所以素數只需要篩選到 109 的開方 即可。

       對于每個數,遍歷素數數組,如果能被某個素數整除,判斷商是不是一個素數,如果商也是素數結果為Yes,否則為No


114 Telecasting station

      枚舉 + 統計


       題意:在X軸上規定一些有權值a­i的點 (x1,a1), (x1,a1), (x2,a2) .. (xn,an)要求在x坐標上找到一個點X使得 M = |X-x1|*a1 + ... |X-xn|*an 的值最小。

       題解:細心觀察可以發現,M是關于X的一次函數,所以可以確定X必定為 x1...xn中的某個點。

       由于原式中存在絕對值,為了將絕對值化簡,可以將X的區間分段,比如當X xj時可以簡化為 M = (X-x1)*a1 + ...(X-xj)*aj + (xj+1-X)*aj+1 + ... + (xn-X)*an

            = (X * presum[j] - preproduct[j]) + (postproduct[j+1] - X * postsum[j+1])

 
圖2

      這些輔助數組可以通過四次線性掃描得出,復雜度O(n)然后,先將xi遞增排序,枚舉每個xi,計算最小值更新M,就可以在O(n)的時間內求解了。總體算法復雜度為排序復雜度O(n log n)


115 Calendar

       簡單模擬

       題意:給定一個月份M和天數N,求2001的那一天是星期幾。

       題解:可以預處理出每個月的天數,如果天數大于那個月對應的天數或者月數大于12肯定是不可行的,否則將給定月數M之前的 1M-1個月的天數相加再加上N得到SumSum mod 7就可以定位到星期幾了。


116 Index of super-prime

      素數篩選 + 記憶化搜索( 動態規劃)

      
       
題意:超級素數是素數下標也是素數的數,給定一個數N,問這個數能不能被一些超級素數組合出來,如果可以,需要滿足超級素數的個數最少,要求輸出一個方案。

       題解:先將滿足條件的 超級素數 篩選出來,3 5 11 17等等,然后對于輸入的N進行一次記憶化搜索(DP

       例如 N = 15

                     那么N的最優值一定是來自 15-3=12 15-5=10 15-11=4 這三個數,以此類推,遞歸出口是N = 0,每次計算完N就將它保存下來,下次就不用重復計算了。

      所以狀態轉移方程為:DP[i] =  min{ DP [ i - p ] ,  p 為超級素數 } + 1

       因為需要輸出組合的序列,所以每次搜索,需要保存當前最優值的前驅結點,最后一次性輸出即可。

 

117 Counting

      二分求余

       題意:給定N(N < 10001)個數Ai,求其中AiM K的倍數的數的個數。

       題解:對于ab mod c二分求解。 b == 0 則 ab mod c = 1 mod c; 否則 ab mod c = (A2)(B/2) * ( (B%2) ? A : 1 ) mod c;


118 Digital Root

      數學歸納法

       題意:定義f(n)  n的所有數字的和. 如果 f(n) 1位數字,那么它是n "digital root",否則n"digital root" f(n) "digital root"
      給定數組Ai­,求 A1*A2* … *AN  +  A1*A2*…*AN-1 + … +  A1*A2  + A1 "digital root"(N <= 1000)

       題解:定義d(n)n"digital root"利用數學歸納法可證明(N=1的情況網上遞推)

             1) d( A1*A2* … *AN ) = d( AN * d( A1*A2* … *AN-1 ) )

             2) d(A1 + A2 ) = d( d(A1) +d(A2) )

       利用這兩點,可以直接掃描A數組就可以計算出給定表達式的"digital root"了。

 

119 Magic Pairs

      枚舉 + 擴展歐幾里得

 

       題意:對于所有的(X,Y)在滿足 (A0 * X + B0 * Y) % N = 0 (1) 的情況下,使得(A * X + B * Y) % N = 0 (2) 也成立, 求這樣的 (A, B)

       題解:首先求出A0B0N三者的GCD,最后算出來的答案需要乘上這個GCD

       (A0 * X + B0 * Y) = K * N  (3)  (K為整數)

       (A * X + B * Y) = K' * N   (4) (K'為整數)

       X = (K*N - B0*Y) / A0

       代入(4),可得

       A*(K*N - B0*Y) + B*A0*Y = K'* N * A0 (5)

       化簡得 (6)

       (B*A0 - A*B0) * Y = (K'*A0 - K*A) * N (6)

       兩邊同時mod N

         (B*A0 - A*B0) * Y % N = 0 (7)

       由于Y為任意整數(Y有可能和N互素), 所以必須滿足  (B*A0 - A*B0) % N = 0  (8)

 

       可以化簡為  (B*A0 - A*B0) = K'' * N  (9)  (K''為整數)

       A0 * B + N * (-K'') = A*B0           (10)

       枚舉A的值,就可以把方程轉化成了 ax + by = c的形式,其中:

       a = A0, b = N, c = A*B0

       x = B, y = -K''

       利用擴展歐幾里得即可求得最小的滿足條件的x(B)的值了。

       最后的答案需要乘上 A0B0N 三者的GCD

 

posted on 2014-06-12 12:49 英雄哪里出來 閱讀(1430) 評論(0)  編輯 收藏 引用 所屬分類: Sgu 題解

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            亚洲大片免费看| 欧美日韩一区成人| 一区在线视频观看| 欧美a级片网| 欧美成人一区二区三区片免费| 在线成人免费观看| 最新日韩在线| 欧美精品一区二区三区很污很色的| 亚洲精品免费一区二区三区| 亚洲精品免费在线播放| 国产精品第2页| 久久久久成人精品| 欧美电影免费观看网站| 亚洲欧美日韩精品久久亚洲区| 亚洲一级高清| 久久国产精品网站| 亚洲精品一区在线| 亚洲在线成人精品| 亚洲国产高清在线| 9色porny自拍视频一区二区| 国产欧美精品在线播放| 欧美电影资源| 欧美图区在线视频| 欧美aaaaaaaa牛牛影院| 国产精品成人午夜| 欧美激情一区二区三区在线视频| 欧美特黄a级高清免费大片a级| 久久激情五月激情| 欧美噜噜久久久xxx| 久久精品一二三区| 欧美日韩在线播放三区四区| 久久久之久亚州精品露出| 欧美日韩成人一区| 你懂的网址国产 欧美| 欧美日韩在线免费| 欧美激情第五页| 国产偷国产偷亚洲高清97cao| 亚洲欧洲在线视频| 在线精品观看| 性18欧美另类| 午夜精品短视频| 欧美美女操人视频| 亚洲国产成人精品女人久久久 | 国产精品免费网站在线观看| 欧美激情亚洲自拍| 激情一区二区三区| 欧美一级专区| 欧美一二三区精品| 国产精品国产三级国产aⅴ9色| 亚洲国产成人tv| 亚洲国产成人不卡| 久久久久久网址| 久久综合色婷婷| 国内精品久久久久影院 日本资源| 99re6热在线精品视频播放速度| 最新国产成人在线观看| 久久蜜臀精品av| 男女av一区三区二区色多| 国产又爽又黄的激情精品视频| 亚洲网站在线| 亚洲欧美bt| 国产精品日韩精品欧美在线| 亚洲网站啪啪| 亚洲欧美一区二区在线观看| 欧美日韩中文| 一区二区三区 在线观看视| 中文久久精品| 欧美丝袜一区二区| 亚洲一区二区三| 久久国产手机看片| 国内揄拍国内精品久久| 欧美综合国产| 欧美电影在线免费观看网站| 亚洲国产欧美日韩| 欧美激情精品久久久久久久变态| 亚洲第一精品电影| 一区二区三区高清在线| 欧美日韩一级黄| 欧美一级视频免费在线观看| 国产日韩欧美在线观看| 欧美尤物一区| 亚洲高清视频在线观看| av成人免费在线观看| 国产精品国产三级国产aⅴ无密码 国产精品国产三级国产aⅴ入口 | 欧美 日韩 国产一区二区在线视频| 激情久久久久久久| 美女亚洲精品| 99这里有精品| 久久男女视频| 一本一道久久综合狠狠老精东影业| 欧美视频一区二区| 欧美亚洲视频在线观看| 欧美高清视频一区二区| 亚洲视频www| 黄色国产精品一区二区三区| 欧美激情第8页| 午夜精品婷婷| 亚洲国产综合视频在线观看| 午夜精品久久久久久99热软件 | 欧美视频精品一区| 欧美淫片网站| 一本色道88久久加勒比精品| 久久久亚洲高清| 99视频精品全部免费在线| 国产日韩欧美二区| 欧美精品一区二区三区很污很色的 | 亚洲欧美日韩精品一区二区| 另类成人小视频在线| 亚洲一区二区黄色| 亚洲国产成人av| 国产情人节一区| 欧美理论电影在线播放| 久久久国产亚洲精品| 制服丝袜激情欧洲亚洲| 亚洲夫妻自拍| 久久免费视频这里只有精品| 一区二区三区国产盗摄| 亚洲狠狠丁香婷婷综合久久久| 国产欧美日韩视频在线观看| 欧美激情综合五月色丁香小说| 香港久久久电影| 99热在线精品观看| 亚洲国产成人一区| 蜜月aⅴ免费一区二区三区| 欧美一区二区三区久久精品茉莉花| 亚洲三级免费电影| 亚洲高清123| 狠狠色综合色综合网络| 国产日韩精品一区二区| 国产精品国产三级国产aⅴ无密码| 欧美精品一区二区三区在线播放 | 国产日韩精品一区观看 | 久久在线免费| 久久精品国产精品亚洲综合| 午夜精品久久久| 亚洲制服av| 久久国产精品亚洲va麻豆| 亚洲在线一区二区三区| 一本色道久久88精品综合| 亚洲精品国产品国语在线app| 亚洲电影第三页| 亚洲国产va精品久久久不卡综合| 在线观看精品一区| 亚洲成色www久久网站| 伊大人香蕉综合8在线视| 国产综合久久久久久| 激情欧美一区二区| 在线精品国产欧美| 亚洲国产精品激情在线观看| 91久久亚洲| av不卡免费看| 亚洲永久在线观看| 欧美一区二区三区在线观看| 久久av二区| 蜜桃伊人久久| 亚洲三级毛片| 亚洲一区二区在线视频 | 一区二区日韩| 久久成人久久爱| 免费久久精品视频| 欧美日本精品在线| 国产精品网红福利| 韩国精品在线观看| 亚洲精品一区二区三区av| 一区二区三区三区在线| 欧美一区激情| 欧美高清不卡在线| 亚洲精品乱码| 午夜久久影院| 欧美黄色免费| 国产乱码精品1区2区3区| 影音先锋久久精品| 宅男噜噜噜66国产日韩在线观看| 欧美亚洲一级| 亚洲国产日韩欧美在线图片| 亚洲视频一区二区免费在线观看| 欧美在线不卡| 欧美日韩999| 激情五月婷婷综合| 一区二区三区偷拍| 久久亚洲精品一区| 日韩一本二本av| 久久亚洲视频| 国产日韩亚洲欧美| 夜夜爽99久久国产综合精品女不卡| 欧美亚洲免费电影| 亚洲国产三级| 久久久欧美精品| 国产精品久久久久免费a∨大胸| 亚洲国产欧美另类丝袜| 久久久99精品免费观看不卡| 亚洲伦理网站| 免费久久99精品国产自| 国产一区二区av| 亚洲综合精品一区二区| 91久久久亚洲精品| 久久免费高清| 黑人巨大精品欧美黑白配亚洲| 亚洲夜间福利| 亚洲免费观看在线视频|