青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

糯米

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 糯米 閱讀(643) 評論(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 | 糯米
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            国产精品入口福利| 亚洲欧洲精品一区二区精品久久久 | 亚洲综合色丁香婷婷六月图片| 欧美激情一区二区三区在线视频| 久久久久久久久久久久久久一区| 在线观看成人小视频| 久久色在线观看| 免费观看成人| 亚洲图片欧美一区| 亚洲在线中文字幕| 激情久久久久久久| 亚洲激情成人网| 欧美日韩免费观看一区| 亚洲欧美在线磁力| 久久久久成人精品免费播放动漫| 亚洲激情图片小说视频| 亚洲精品视频一区二区三区| 国产精品久久毛片a| 久久久av毛片精品| 欧美人与性动交a欧美精品| 亚洲综合色视频| 狼人社综合社区| 午夜精彩国产免费不卡不顿大片| 性久久久久久久久| 亚洲巨乳在线| 欧美在线视频观看免费网站| 亚洲级视频在线观看免费1级| 一本一本久久a久久精品综合妖精| 国产视频久久久久久久| 亚洲国产精品一区二区尤物区| 国产精品地址| 欧美成人资源网| 国产精品免费网站| 亚洲国产一区二区a毛片| 国产精品婷婷午夜在线观看| 欧美大片免费久久精品三p| 国产精品啊啊啊| 欧美激情亚洲精品| 国产一区在线视频| 亚洲美女中文字幕| 亚洲第一页在线| 午夜综合激情| 性欧美暴力猛交另类hd| 欧美伦理a级免费电影| 美女脱光内衣内裤视频久久影院| 国产精品成人一区| 亚洲理论在线观看| 亚洲国产美女精品久久久久∴| 午夜久久久久久久久久一区二区| 一区二区三区视频在线看| 久久久久女教师免费一区| 欧美一级欧美一级在线播放| 欧美女主播在线| 亚洲高清毛片| 亚洲国内高清视频| 久久久久一区二区三区四区| 久久激情婷婷| 国产精品一区二区在线| 99精品国产在热久久| 亚洲美洲欧洲综合国产一区| 免费在线播放第一区高清av| 六月天综合网| 亚洲成色精品| 麻豆精品一区二区av白丝在线| 久久亚洲综合| 一区在线电影| 久久一区激情| 欧美福利视频一区| 亚洲人体大胆视频| 欧美精品精品一区| 亚洲精品影院在线观看| 一区二区高清在线| 欧美午夜剧场| 国产精品99久久久久久久久久久久 | 欧美四级电影网站| 一区二区三区欧美日韩| 亚洲图中文字幕| 国产精品一级| 午夜精品短视频| 久久久久久久综合| 91久久国产综合久久蜜月精品| 嫩草伊人久久精品少妇av杨幂| 欧美黑人国产人伦爽爽爽| 日韩网站在线观看| 国产精品久久久久久久一区探花| 亚洲综合视频一区| 久久香蕉国产线看观看av| 亚洲国产日韩精品| 欧美日本韩国在线| 中文日韩电影网站| 久久亚洲美女| 99热在线精品观看| 国产精品日韩一区二区| 欧美一区二区三区精品电影| 免费永久网站黄欧美| 亚洲视频在线观看免费| 国产视频一区在线| 免费欧美在线| 亚洲综合视频一区| 亚洲高清在线观看一区| 亚洲网站在线看| 国外成人在线| 欧美午夜一区| 久久男人资源视频| 一区二区电影免费观看| 久久夜色精品国产| 亚洲调教视频在线观看| 黄页网站一区| 国产精品美女久久久免费| 久久野战av| 亚洲综合清纯丝袜自拍| 亚洲高清视频的网址| 久久久国产精品一区二区三区| 亚洲精品久久久久久一区二区| 国产精品一区二区你懂得| 欧美成人免费一级人片100| 亚洲欧美日韩国产成人精品影院| 欧美激情视频一区二区三区不卡| 性一交一乱一区二区洋洋av| 亚洲国产裸拍裸体视频在线观看乱了中文 | 国产亚洲aⅴaaaaaa毛片| 欧美日韩国产精品| 久久久久国产免费免费| 亚洲一区二区三区四区五区午夜 | 亚洲精品欧洲精品| 久久久久久夜精品精品免费| 亚洲香蕉在线观看| 亚洲人成人77777线观看| 国产精品日日摸夜夜添夜夜av| 欧美精品www| 蜜桃精品久久久久久久免费影院| 欧美在线综合| 亚洲免费视频网站| 一区二区三区精品在线| 亚洲国产精品久久久久婷婷884| 老**午夜毛片一区二区三区| 亚洲欧美在线免费观看| 亚洲一区二区三区高清| 99riav久久精品riav| 亚洲精品欧美激情| 亚洲国产一区二区三区在线播| 极品日韩av| 亚洲大片av| 亚洲日本成人| 亚洲精品黄色| 99国产精品久久久久久久成人热| 亚洲国产欧美一区二区三区丁香婷| 一区二区在线观看av| 国精品一区二区三区| 国内久久婷婷综合| 国内精品久久久久伊人av| 狠狠色伊人亚洲综合网站色| 国产日韩亚洲欧美综合| 国内久久婷婷综合| 极品尤物一区二区三区| 在线观看欧美| 亚洲日本中文字幕| 夜夜嗨一区二区| 亚洲一级黄色| 久久国产精品黑丝| 久久免费视频在线| 亚洲成人在线免费| 国产精品qvod| 国产视频欧美视频| 亚洲第一黄色| av不卡在线观看| 欧美一区二区三区喷汁尤物| 久久综合导航| 亚洲日本理论电影| 中国女人久久久| 欧美中文字幕在线观看| 免费久久99精品国产自在现线| 欧美另类亚洲| 国产亚洲成精品久久| 亚洲福利在线看| 亚洲午夜电影网| 久久久另类综合| 亚洲国产精品久久人人爱蜜臀| 一本一本久久a久久精品综合麻豆| 午夜国产精品视频| 欧美大尺度在线观看| 国产精品毛片大码女人| 亚洲国产精品欧美一二99| 亚洲一区二区在| 久久综合综合久久综合| 亚洲精品网址在线观看| 欧美一区二区三区在| 欧美日韩精品三区| 韩国v欧美v日本v亚洲v| 亚洲视频 欧洲视频| 欧美aaa级| 亚洲一区免费看| 欧美激情免费在线| 国产真实乱偷精品视频免| 一区二区免费在线播放| 久久综合网络一区二区| 一本久久精品一区二区| 欧美成人精品激情在线观看| 国产三级精品三级| 亚洲天堂久久|