• <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 糯米 閱讀(626) 評論(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 | 糯米
            亚洲欧美一区二区三区久久| 亚洲欧洲日产国码无码久久99| 久久亚洲精品成人AV| 国内精品久久久久| 国产午夜福利精品久久| 久久久久久国产精品免费免费| 久久综合亚洲色HEZYO社区| 久久精品国产免费观看| 久久免费高清视频| 2019久久久高清456| 99久久精品费精品国产 | 久久精品国产亚洲7777| 国产成人精品综合久久久久| 94久久国产乱子伦精品免费| 久久综合亚洲色HEZYO社区| 久久久久免费精品国产| 久久国产欧美日韩精品| 国产成人综合久久精品尤物| 97精品伊人久久大香线蕉| 99久久无码一区人妻| 人妻少妇久久中文字幕一区二区 | 亚洲中文字幕无码一久久区| 久久中文娱乐网| 亚洲国产精品高清久久久| 很黄很污的网站久久mimi色| 国内精品久久久久久野外| 久久精品国产亚洲精品2020 | 色妞色综合久久夜夜| 一级做a爰片久久毛片看看| 狠狠久久综合伊人不卡| 香蕉久久一区二区不卡无毒影院| 午夜天堂av天堂久久久| 久久久久久久波多野结衣高潮| 久久久久一本毛久久久| 国产综合免费精品久久久| 久久99精品久久久久久齐齐| 久久精品国内一区二区三区| 97久久超碰国产精品旧版| 久久精品成人免费网站| 99久久99久久精品国产| 久久亚洲国产成人影院网站|