• <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>

            superman

            聚精會神搞建設 一心一意謀發展
            posts - 190, comments - 17, trackbacks - 0, articles - 0
               :: 首頁 :: 新隨筆 :: 聯系 :: 聚合  :: 管理

            Section 1.4 - Arithmetic Progressions

            Posted on 2009-03-21 18:54 superman 閱讀(125) 評論(0)  編輯 收藏 引用 所屬分類: USACO
             1 #include <iostream>
             2 
             3 using namespace std;
             4 
             5 struct Record
             6 {
             7     int a, b;
             8 
             9     bool operator < (const Record & r) const
            10     {
            11         if (b == r.b)
            12             return a < r.a;
            13         return b < r.b;
            14     }
            15 }   rec[10000];
            16 
            17 int recCnt;
            18 
            19 int n, m;
            20 int bisqCnt, bisq[125000 + 1];
            21 bool isBisq[125000 + 1];
            22 
            23 int main()
            24 {
            25     freopen("ariprog.in""r", stdin);
            26     freopen("ariprog.out""w", stdout);
            27 
            28     cin >> n >> m;
            29 
            30     for (int p = 0; p <= m; p++)
            31     for (int q = 0; q <= m; q++)
            32         isBisq[p * p + q * q] = true;
            33 
            34     for (int i = 0; i <= 2 * m * m; i++)
            35         if (isBisq[i])
            36             bisq[bisqCnt++= i;
            37 
            38     for (int i = 0; i <= bisqCnt - n; i++)
            39     {
            40         for (int j = i + 1; j <= bisqCnt - n + 1; j++)
            41         {
            42             int a = bisq[i];
            43             int b = bisq[j] - bisq[i];
            44 
            45             int k;
            46             for (k = 2; k < n; k++)
            47             {
            48                 if (a + b * k > 2 * m * m)
            49                     break;
            50                 if (isBisq[a + b * k] == false)
            51                     break;
            52             }
            53             if (k == n)
            54             {
            55                 rec[recCnt].a = a;
            56                 rec[recCnt].b = b;
            57                 recCnt++;
            58             }
            59         }
            60     }
            61 
            62     sort(rec, rec + recCnt);
            63 
            64     if (recCnt == 0)
            65         cout << "NONE" << endl;
            66 
            67     for (int i = 0; i < recCnt; i++)
            68         cout << rec[i].a << ' ' << rec[i].b << endl;
            69 
            70     return 0;
            71 }
            72 
            久久天天躁狠狠躁夜夜不卡 | 久久综合九色综合欧美就去吻 | 国产精品亚洲综合久久| 伊人热热久久原色播放www | 99久久国产综合精品五月天喷水 | 久久精品无码专区免费东京热| 新狼窝色AV性久久久久久| 久久精品国产亚洲欧美| 欧美午夜A∨大片久久| 狠狠久久亚洲欧美专区| 久久婷婷色香五月综合激情 | 国产精品成人精品久久久| 国产精品久久新婚兰兰| 久久高潮一级毛片免费| 精品免费久久久久久久| 久久精品卫校国产小美女| 99久久精品无码一区二区毛片 | 久久亚洲精品国产精品| 久久AAAA片一区二区| 精品久久久无码21p发布| 欧美va久久久噜噜噜久久| 三级片免费观看久久| 97香蕉久久夜色精品国产| 97久久国产综合精品女不卡| 久久国产乱子精品免费女| 色婷婷久久综合中文久久蜜桃av| 婷婷综合久久狠狠色99h| 国内精品人妻无码久久久影院| 久久99热这里只有精品66| 久久亚洲日韩看片无码| 色综合久久久久无码专区| 久久水蜜桃亚洲av无码精品麻豆| 久久久黄色大片| 91久久精品91久久性色| 久久亚洲国产精品五月天婷| 久久精品日日躁夜夜躁欧美| 久久这里只精品国产99热| 久久精品国产日本波多野结衣| 久久久久久狠狠丁香| 久久精品中文无码资源站| 国产精品一区二区久久精品涩爱|