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

            糯米

            TI DaVinci, gstreamer, ffmpeg
            隨筆 - 167, 文章 - 0, 評論 - 47, 引用 - 0
            數據加載中……

            POJ 2132 Cow Math 二分

            思路:
            首先每條路徑的值都可以分解一下質因數,就可以表示為多個質數的冪相乘的形式,
             比如 2^6 * 3^8 * 17^22 * 23^1。

            三個數字a, b, c求最大公約數,分解完質因數后:
            如果a擁有2^8,b擁有2^10,c擁有2^4。那最大公約數必然擁有2^4,取最小的一個。
            對于每個質數 2, 3, 5, 7。。都是這個道理。

            如果是求最小公倍數,在剛剛的例子里,就是取最大的一個了。

            在點之間行走的過程,可以這樣來看。在點1的時候GCF的值是所有質數的最大次冪的乘積。
            GCF的值必定是越走越小。
            每經過一條路徑,CGF各個質因數的冪都必須小于等于路徑的對應的值。
            就好比路徑就只能容納這么大的流量。然后到達點2的時候,看看哪條路徑的流量最大。
            看起來像最大流問題,但不是最大流問題。

            我們沒辦法遍歷一次圖,就求出哪條路徑的流量最大。
            但由于路徑的權值最大才2000,質因數的冪最大也只有11(2^11 = 2048),大不了每個冪都試一次。
            用二分法就可以了。

            對于每一個質數,求到達點2 的時候的最大的冪。
            最后再乘起來,就是答案了。
            可見這種方法還是很巧妙的,效率也很高,0ms AC。

            注意:
            不需要高精度。但需要用__int64來保存答案。

            #include <stdio.h>

            #define MAX_W 2048
            #define MAX_N 32 

            int N, visit[MAX_N], map[MAX_N][MAX_N], tm;
            int prime[MAX_W], prime_cnt, max_cnt[MAX_W];

            int dfs(int idx, int val, int cnt)
            {
                
            int i, j, k;

                
            if (idx == 2)
                    
            return 1;

                visit[idx] 
            = tm;
                
            for (i = 1; i <= N; i++{
                    
            if (visit[i] == tm)
                        
            continue;
                    j 
            = map[idx][i];
                    
            for (k = 0; j && !(j % val); k++)
                        j 
            /= val;
                    
            if (k < cnt)
                        
            continue;
                    
            if (dfs(i, val, cnt))
                        
            return 1;
                }


                
            return 0;
            }


            __inline 
            int calc(int val, int r)
            {
                
            int l, m;

                l 
            = 0;
                
            while (l <= r) {
                    m 
            = (l + r) / 2;
                    tm
            ++;
                    
            if (dfs(1, val, m))
                        l 
            = m + 1;
                    
            else
                        r 
            = m - 1;
                }


                
            return r;
            }


            int main()
            {
                
            int i, j, val, p, cnt;
                __int64 r;

                freopen(
            "e:\\test\\in.txt""r", stdin);

                prime[prime_cnt
            ++= 2;
                
            for (i = 3; i < MAX_W; i++{
                    
            for (j = 0; j < prime_cnt && (i % prime[j]); j++);
                    
            if (j == prime_cnt)
                        prime[prime_cnt
            ++= i;
                }

                
                scanf(
            "%d"&N);
                
            for (i = 1; i <= N; i++)
                    
            for (j = 1; j <= N; j++)
                        scanf(
            "%d"&map[i][j]);
                
                
            for (i = 2; i <= N; i++{
                    val 
            = map[1][i];
                    
            for (j = 0; j < prime_cnt && val >= 1; j++{
                        p 
            = prime[j];
                        
            for (cnt = 0!(val % p); cnt++)
                            val 
            /= p;
                        
            if (cnt > max_cnt[j])
                            max_cnt[j] 
            = cnt;
                    }

                }

                
                
            for (i = 0; i < prime_cnt; i++{
                    
            if (!max_cnt[i])
                        
            continue;
                    max_cnt[i] 
            = calc(prime[i], max_cnt[i]);
                }


                r 
            = 1;
                
            for (i = 0; i < prime_cnt; i++{
                    
            if (!max_cnt[i])
                        
            continue;
                    
            for (cnt = 0; cnt < max_cnt[i]; cnt++)
                        r 
            *= prime[i];
                }

                printf(
            "%I64d\n", r);

                
            return 0;
            }

            posted on 2010-03-14 14:37 糯米 閱讀(623) 評論(1)  編輯 收藏 引用 所屬分類: POJ

            評論

            # re: POJ 2132 Cow Math 二分[未登錄]  回復  更多評論   

            POJ 2132 Cow Math 二分

            思路:
            首先每條路徑的值都可以分解一下質因數,就可以表示為多個質數的冪相乘的形式,
            比如 2^6 * 3^8 * 17^22 * 23^1。

            三個數字a, b, c求最大公約數,分解完質因數后:
            如果a擁有2^8,b擁有2^10,c擁有2^4。那最大公約數必然擁有2^4,取最小的一個。
            對于每個質數 2, 3, 5, 7。。都是這個道理。

            如果是求最小公倍數,在剛剛的例子里,就是取最大的一個了。

            在點之間行走的過程,可以這樣來看。在點1的時候GCF的值是所有質數的最大次冪的乘積。
            GCF的值必定是越走越小。
            每經過一條路徑,CGF各個質因數的冪都必須小于等于路徑的對應的值。
            就好比路徑就只能容納這么大的流量。然后到達點2的時候,看看哪條路徑的流量最大。
            看起來像最大流問題,但不是最大流問題。

            我們沒辦法遍歷一次圖,就求出哪條路徑的流量最大。
            但由于路徑的權值最大才2000,質因數的冪最大也只有11(2^11 = 2048),大不了每個冪都試一次。
            用二分法就可以了。

            對于每一個質數,求到達點2 的時候的最大的冪。
            最后再乘起來,就是答案了。
            可見這種方法還是很巧妙的,效率也很高,0ms AC。

            注意:
            不需要高精度。但需要用__int64來保存答案。


            #include <stdio.h>

            #define MAX_W 2048
            #define MAX_N 32

            int N, visit[MAX_N], map[MAX_N][MAX_N], tm;
            int prime[MAX_W], prime_cnt, max_cnt[MAX_W];

            int dfs(int idx, int val, int cnt)
            {
            int i, j, k;

            if (idx == 2)
            return 1;

            visit[idx] = tm;
            for (i = 1; i <= N; i++) {
            if (visit[i] == tm)
            continue;
            j = map[idx][i];
            for (k = 0; j && !(j % val); k++)
            j /= val;
            if (k < cnt)
            continue;
            if (dfs(i, val, cnt))
            return 1;
            }

            return 0;
            }

            __inline int calc(int val, int r)
            {
            int l, m;

            l = 0;
            while (l <= r) {
            m = (l + r) / 2;
            tm++;
            if (dfs(1, val, m))
            l = m + 1;
            else
            r = m - 1;
            }

            return r;
            }

            int main()
            {
            int i, j, val, p, cnt;
            __int64 r;

            freopen("e:\\test\\in.txt", "r", stdin);

            prime[prime_cnt++] = 2;
            for (i = 3; i < MAX_W; i++) {
            for (j = 0; j < prime_cnt && (i % prime[j]); j++);
            if (j == prime_cnt)
            prime[prime_cnt++] = i;
            }

            scanf("%d", &N);
            for (i = 1; i <= N; i++)
            for (j = 1; j <= N; j++)
            scanf("%d", &map[i][j]);

            for (i = 2; i <= N; i++) {
            val = map[1][i];
            for (j = 0; j < prime_cnt && val >= 1; j++) {
            p = prime[j];
            for (cnt = 0; !(val % p); cnt++)
            val /= p;
            if (cnt > max_cnt[j])
            max_cnt[j] = cnt;
            }
            }

            for (i = 0; i < prime_cnt; i++) {
            if (!max_cnt[i])
            continue;
            max_cnt[i] = calc(prime[i], max_cnt[i]);
            }

            r = 1;
            for (i = 0; i < prime_cnt; i++) {
            if (!max_cnt[i])
            continue;
            for (cnt = 0; cnt < max_cnt[i]; cnt++)
            r *= prime[i];
            }
            printf("%I64d\n", r);

            return 0;
            }
            2014-08-07 14:52 | 糯米
            色综合久久久久网| 97精品伊人久久大香线蕉app| AAA级久久久精品无码片| 精品国产乱码久久久久久呢| 麻豆av久久av盛宴av| 亚洲日韩欧美一区久久久久我| 亚洲国产成人精品女人久久久 | 色综合久久久久网| 久久综合狠狠综合久久激情 | 国产精品一久久香蕉国产线看| 久久九九精品99国产精品| 久久99久久99小草精品免视看| 久久综合中文字幕| 久久精品国产亚洲5555| 久久久久久久综合日本| 久久亚洲AV无码精品色午夜麻豆| 久久这里只有精品首页| av无码久久久久久不卡网站| 97久久精品人人澡人人爽| 亚洲国产精品一区二区三区久久| 亚洲午夜久久久久久久久久| 久久99精品久久久久久动态图| 办公室久久精品| 久久久国产精华液| AA级片免费看视频久久| 国产亚洲精品久久久久秋霞| 91精品国产综合久久婷婷| 久久人妻少妇嫩草AV无码蜜桃| 午夜天堂精品久久久久| 久久精品一区二区三区中文字幕| 日韩精品久久久久久免费| 久久精品国产亚洲Aⅴ香蕉| 午夜天堂精品久久久久| 久久人人超碰精品CAOPOREN| 99久久超碰中文字幕伊人| 久久亚洲精品国产精品婷婷| 欧美日韩中文字幕久久伊人| 无码精品久久久天天影视| 天堂无码久久综合东京热| 94久久国产乱子伦精品免费 | 日韩精品久久久肉伦网站|