• <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 
            人妻久久久一区二区三区| 久久福利青草精品资源站| 欧美粉嫩小泬久久久久久久 | 亚洲AV伊人久久青青草原| 精品久久久一二三区| 久久99国产综合精品女同| Xx性欧美肥妇精品久久久久久| 色婷婷噜噜久久国产精品12p| 亚洲精品乱码久久久久久蜜桃不卡| 久久99精品国产| 日韩中文久久| 好属妞这里只有精品久久| 久久九色综合九色99伊人| 性欧美丰满熟妇XXXX性久久久| 国产午夜精品久久久久九九电影| 久久精品极品盛宴观看| 久久久国产精品福利免费| 精品久久久久久国产| 久久国产免费| 国产69精品久久久久99| 中文字幕无码免费久久| 亚洲国产综合久久天堂| 99精品久久久久久久婷婷| 潮喷大喷水系列无码久久精品| 伊人色综合九久久天天蜜桃| 欧美精品一本久久男人的天堂| 久久精品中文字幕一区| 久久久久久国产a免费观看黄色大片| 精品久久一区二区三区| 人妻无码久久一区二区三区免费 | 久久99精品国产99久久| 国产亚洲精品久久久久秋霞 | 午夜不卡久久精品无码免费| 日韩久久久久中文字幕人妻| 99热热久久这里只有精品68| 久久精品男人影院| 久久精品国产福利国产秒| 久久久综合九色合综国产| 99热成人精品免费久久| 久久国产免费直播| 亚洲综合久久夜AV |