• <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 糯米 閱讀(636) 評論(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 | 糯米
            久久久国产打桩机| 伊色综合久久之综合久久| 精品国产乱码久久久久软件| 成人国内精品久久久久影院VR| 久久久久久毛片免费播放| 天天爽天天狠久久久综合麻豆| 久久国产免费直播| 漂亮人妻被黑人久久精品| 99久久99久久精品国产片果冻| 久久天天躁夜夜躁狠狠躁2022| 久久天天躁狠狠躁夜夜躁2014| 亚洲综合伊人久久大杳蕉| 日产精品久久久久久久| 久久国产精品-国产精品| 麻豆精品久久精品色综合| 99久久精品国产一区二区| 久久久久18| 久久香蕉国产线看观看精品yw| 97r久久精品国产99国产精| 91久久精品国产免费直播| 久久久久亚洲爆乳少妇无| 久久午夜免费视频| 热re99久久精品国99热| 久久久精品一区二区三区| 亚洲欧美另类日本久久国产真实乱对白 | 精品国产乱码久久久久久呢| 亚洲国产另类久久久精品小说| 99国产欧美精品久久久蜜芽| 日韩精品久久久久久| 热综合一本伊人久久精品 | 久久国产乱子伦免费精品| 青青草原1769久久免费播放| 人妻无码精品久久亚瑟影视| 久久久一本精品99久久精品88| 久久精品亚洲欧美日韩久久| 嫩草伊人久久精品少妇AV| 欧美性猛交xxxx免费看久久久 | 亚洲第一极品精品无码久久| 国産精品久久久久久久| 久久亚洲春色中文字幕久久久| 午夜精品久久久久成人|