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

oyjpArt ACM/ICPC算法程序設計空間

// I am new in programming, welcome to my blog
I am oyjpart(alpc12, 四城)
posts - 224, comments - 694, trackbacks - 0, articles - 6

基本參數搜索

Posted on 2008-06-03 15:45 oyjpart 閱讀(3140) 評論(14)  編輯 收藏 引用 所屬分類: ACM/ICPC或其他比賽

上次百度之星第三題竟然不會做,很是慚愧啊,腦袋生銹了。

后來從HUST上面找了道類似的題目,AC了。


The perfect hamilton path

Time Limit: 5 Sec  Memory Limit: 128 MB
Submissions: 72  Solved: 16

Description

There are N(2 <= N <= 13) cities and M bidirectional roads among the cities. There exist at most one road between any pair of the cities. Along every road, there are G pretty girls and B pretty boys(1 <= G,B <= 1000).
You want to visit every city exactly once, and you can start from any city you want to. The degree of satisfaction is the ratio of the number of the pretty girls to the number of the pretty boys. You want to know the highest degree of satisfation.

Input

There are multiply test cases.
First line: two integers N, M;
The following M lines: every line with four integers i, j, G, B, response that there is a road between i and j with G and B.

Output

The highest degree of the satisfation, rounded to the third place after the decimal point.

Sample Input

3 3
1 2 5 3
2 3 7 4
3 1 13 11

Sample Output

1.714

HINT

Source

dupeng


題目的意思是找到一個sigma(G)/sigma(B)最大的hamilton回路。
典型的參數搜索。二分或者迭代答案就可以了。

Solution:

#include <stdio.h>
#include 
<queue>
#include 
<cmath>
using namespace std;

const double EPS = 1e-4;
const int N = 15;
const int M = N * N;

#define Max(a, b) (a
>b?a:b)

inline 
int dblcmp(double a, double b) {
    
if(fabs(a-b) < EPS) return 0;
    
return a < b ? -1 : 1;
}

struct Node 
{
    
int x, mask;
    
double s;
    Node() {}
    Node(
int mm, int xx, double ss) {
        x 
= xx;
        mask 
= mm;
        s 
= ss;
    }
};

int n, m;

double adj[N][N];
int X[M], Y[M], G[M], B[M];

double dp[1<<N][N];

double go(double ans) {
    
int i, j;
    
for(i = 0; i < n; ++i) {
        adj[i][i] 
= 0;
        
for(j = i+1; j < n; ++j) {
            adj[i][j] 
= adj[j][i] = -10e300;
        }
    }
    
for(i = 0; i < m; ++i) {
        adj[X[i]
-1][Y[i]-1= G[i]-ans * B[i];
        adj[Y[i]
-1][X[i]-1= adj[X[i]-1][Y[i]-1];
    }

    
for(i = 0; i < (1<<n); ++i) {
        
for(j = 0; j < n; ++j)
            dp[i][j] 
= -10e100;
    }
    queue
<Node> Q;
    
for(i = 0; i < n; ++i) {
        Q.push(Node(
1<<i, i, 0.0));
        dp[
1<<i][i] = 0;
    }
    
while(Q.size()) {
        
int f = Q.front().mask, x = Q.front().x;
        
double s = Q.front().s;
        
double& d = dp[f][x];
        Q.pop();
        
if(s < d) continue;
        
for(i = 0; i < n; ++i) if((f&(1<<i)) == 0) {
            
if(dp[f|1<<i][i] < s + adj[x][i]) {
                dp[f
|1<<i][i] = s + adj[x][i];
                Q.push(Node(f
|1<<i, i, s + adj[x][i]));
            }
        }
    }

    
double max = -10e100;
    
for(i = 0; i < n; ++i) {
        max 
= Max(max, dp[(1<<n)-1][i]);
    }
    
return max;
}

int main()
{
    
// freopen("t.in", "r", stdin);

    
int i;
    
double ans;
    
while(scanf("%d %d"&n, &m) != EOF) {
        
double min = 2000, max = 0;
        
for(i = 0; i < m; ++i) {
            scanf(
"%d %d %d %d"&X[i], &Y[i], &G[i], &B[i]);
            
if(B[i] < min) min = B[i];
            
if(G[i] > max) max = G[i];
        }
        
double lo = 0, hi = max/min;
        
int ok = 0;
        
for(i = 0; ; ++i) {
            
double mid = lo + (hi-lo)/2;
            
if(dblcmp((ans=go(mid)), 0.0> 0) {
                lo 
= mid;
            } 
else if(dblcmp(ans, 0.0== 0) {
                printf(
"%.3lf\n", mid);
                ok 
= 1;
                
break;
            } 
else {
                hi 
= mid;
            }
        }

        
if(!ok) { int a = 0; a = 1/a; }
    }

    
return 0;
}

 


Feedback

# re: 基本參數搜索  回復  更多評論   

2008-06-04 13:43 by w
你好,這個程序我看不懂……能講一下思路嗎?

# re: 基本參數搜索  回復  更多評論   

2008-06-04 14:56 by oyjpart
你可以參考《算法藝術與信息學競賽》303-304頁
3.地震--最有比率生成樹 一節的解答
和這個非常類似

就是2分枚舉那個答案,然后將除的表達式的權 轉化成+-*表達式的權,再這個基礎上求目標函數。 如果目標函數 != 0,則枚舉的答案應該向使目標函數更接近0的方向取值,

go函數實際求的就是最大權的hamilton回路。用的是基本的壓縮狀態廣搜。

# re: 基本參數搜索  回復  更多評論   

2008-06-04 15:02 by Surfing
我的解法

#include <stdio.h>

#define N 13

typedef struct _T_AdjNode
{
int nBoys;
int nGirls;
double dRatio;
}TAdjNode;

TAdjNode g_AdjNode[N][N];
int g_Path[2][N];
int g_PathIndex[2] = {0};
double g_dRatio[2] = {0.0};
int nCities, nRoads;

int FindNextNode(int nPathIndex, int nLine)
{
double dRatio = 0;
int nNode = 0;
int i = 0;
int j = 0;
bool bExist = false;

for (j = 0; j < nCities; j++)
{
for (i = 0; i < g_PathIndex[nPathIndex]; i++)
{
if (j == g_Path[nPathIndex][i])
{
bExist = true;
break;
}
}
if (bExist)
{
bExist = false;
continue;
}
if (g_AdjNode[nLine][j].dRatio > dRatio)
{
dRatio = g_AdjNode[nLine][j].dRatio;
nNode = j;
}
}

return nNode;
}

int FindPath(int nPathIndex, int nNode)
{
int nNextNode = 0;
static int nBoys = 0, nGirls = 0;

g_Path[nPathIndex][g_PathIndex[nPathIndex]] = nNode;
g_PathIndex[nPathIndex]++;
if (g_PathIndex[nPathIndex] >= nCities)
{
g_dRatio[nPathIndex] = (double)nGirls / nBoys;
return 0;
}

nNextNode = FindNextNode(nPathIndex, nNode);
nBoys += g_AdjNode[nNode][nNextNode].nBoys;
nGirls += g_AdjNode[nNode][nNextNode].nGirls;
FindPath(nPathIndex, nNextNode);

return 0;
}

int main()
{
int i,j,nGirls,nBoys;
char q = '0';
int nPathIndex = 0;

nCities = nRoads = 0;
i = j = nGirls = nBoys = 0;

printf("Input the number of cities and roads:\n");
scanf("%d %d", &nCities, &nRoads);

if (nCities < 1 || nRoads < 1)
{
return 1;
}

do
{
printf("Input the road index and the number of girls and boys sequentially : "
"from to girls boys\n");
scanf("%d %d %d %d", &i, &j, &nGirls, &nBoys);
getchar();

g_AdjNode[i - 1][j - 1].nBoys = nBoys;
g_AdjNode[i - 1][j - 1].nGirls = nGirls;
g_AdjNode[i - 1][j - 1].dRatio = (double)nGirls / nBoys;
g_AdjNode[j - 1][i - 1].nBoys = nBoys;
g_AdjNode[j - 1][i - 1].nGirls = nGirls;
g_AdjNode[j - 1][i - 1].dRatio = g_AdjNode[i - 1][j - 1].dRatio;

printf("Input finished?(y/n)");
scanf("%c", &q);
getchar();
} while ('y' != q);

//process here
nPathIndex = 0;
for (i = 0; i < nCities; i++)
{
FindPath(nPathIndex, 0);
nPathIndex = g_dRatio[0] <= g_dRatio[1] ? 0 : 1;
g_PathIndex[nPathIndex] = 0;
}

//output the result
nPathIndex = g_dRatio[0] >= g_dRatio[1] ? 0 : 1;
printf("The max ratio is %.3lf\n", g_dRatio[nPathIndex]);\
printf("The best path : \n");
for (i = 0; i < nCities; i++)
{
printf("%d\t", g_Path[nPathIndex][i]);
}
printf("\n");

return 0;
}

# re: 基本參數搜索  回復  更多評論   

2008-06-04 15:10 by Surfing
一點小問題,更正一下

if (g_PathIndex[nPathIndex] >= nCities)
{
g_dRatio[nPathIndex] = (double)nGirls / nBoys;
nGirls = nBoys = 0;
return 0;
}

# re: 基本參數搜索  回復  更多評論   

2008-06-04 17:13 by oyjpart
@Surfing
嘿嘿,謝謝分享

# re: 基本參數搜索  回復  更多評論   

2008-06-05 22:27 by w
多謝,受教了

# re: 基本參數搜索  回復  更多評論   

2008-06-05 23:07 by oyjpart
不謝

# re: 基本參數搜索  回復  更多評論   

2008-06-09 23:54 by richardxx
我做了百度那題,但比賽完才想起我貼的那個模版有點問題,最后果然只有4.5分,和沒做沒區別~~

# re: 基本參數搜索  回復  更多評論   

2008-06-10 12:03 by oyjpart
@richardxx
呵呵 進復賽了就可以了不 看我們這種初賽就被水掉的菜菜。。

# re: 基本參數搜索  回復  更多評論   

2008-06-10 20:01 by 小Young
跟著大牛漲經驗值!

# re: 基本參數搜索  回復  更多評論   

2008-06-10 20:34 by oyjpart
汗。。。
您謙虛了。。。

# re: 基本參數搜索  回復  更多評論   

2008-06-11 19:12 by 小Young
請問這題你用隊列有什么用途啊?
這題不用隊列也可以啊.

# re: 基本參數搜索  回復  更多評論   

2008-06-11 22:19 by oyjpart
@ 小Young
就是廣搜用的隊列
不用隊列你的意思是深搜么?

# re: 基本參數搜索  回復  更多評論   

2008-07-26 06:09 by lengbufang
看看!!!
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            久久爱另类一区二区小说| 亚洲欧美视频在线| 欧美α欧美αv大片| 看欧美日韩国产| 日韩视频在线观看| 一区二区三区日韩精品视频| 国产免费成人在线视频| 老鸭窝亚洲一区二区三区| 玖玖综合伊人| 午夜精品福利一区二区蜜股av| 亚洲一区久久| 亚洲福利一区| 亚洲网站视频| 亚洲国产经典视频| 一区二区欧美在线| 在线高清一区| 亚洲性线免费观看视频成熟| 黄色亚洲免费| 一区二区三区 在线观看视频| 国产伪娘ts一区| 亚洲激情一区| 国产一区二区三区在线观看免费视频| 老司机成人在线视频| 欧美三日本三级少妇三99| 久久米奇亚洲| 欧美视频在线免费| 欧美成人亚洲成人| 国产精品羞羞答答xxdd| 亚洲国产成人不卡| 国产午夜亚洲精品理论片色戒| 亚洲大胆人体在线| 国产在线视频欧美| 在线视频精品| 亚洲美女精品久久| 久久精品国产精品亚洲| 午夜精品免费在线| 欧美日韩国内| 亚洲黄色一区| 伊人久久综合97精品| 午夜一区二区三视频在线观看| 99国产精品久久久久老师| 久久久久久一区二区| 欧美在线一二三四区| 国产精品久久久久久久久搜平片| 亚洲第一级黄色片| 在线播放国产一区中文字幕剧情欧美| 亚洲一区二区3| 亚洲一区二区黄| 欧美精品九九| 亚洲精品乱码久久久久久按摩观| 亚洲国产欧美国产综合一区| 久久国产99| 久久久久久欧美| 国产一区二区三区直播精品电影 | 亚洲精品在线一区二区| 亚洲电影天堂av| 久久久久高清| 免费看的黄色欧美网站| 狠狠色2019综合网| 久久国产黑丝| 美女精品一区| 亚洲国产视频一区二区| 裸体丰满少妇做受久久99精品 | 亚洲精品国产视频| 欧美91精品| 亚洲激情av| 亚洲性感美女99在线| 欧美日韩一区在线观看视频| 99视频热这里只有精品免费| 亚洲视频在线播放| 国产精品久久久久77777| 亚洲午夜av电影| 欧美在线免费视屏| 黄色成人在线| 麻豆精品网站| 99国产精品视频免费观看| 亚洲先锋成人| 国内精品美女av在线播放| 久久婷婷人人澡人人喊人人爽 | 久久精品三级| 亚洲欧洲精品成人久久奇米网| 欧美成人激情在线| 99精品免费网| 久久综合国产精品| 亚洲精品在线观| 国产精品一区2区| 久久久久久婷| 亚洲美女区一区| 久久激情五月丁香伊人| 亚洲激情精品| 国产精品久久久久免费a∨大胸| 欧美亚洲在线视频| 亚洲国产一区二区视频| 香蕉av福利精品导航| 亚洲国产精品一区二区www| 欧美区日韩区| 久久精品一本| 亚洲天堂av高清| 欧美freesex交免费视频| 亚洲深夜激情| 在线观看欧美精品| 国产精品伦理| 欧美日本韩国在线| 久久久精品五月天| 99视频热这里只有精品免费| 麻豆精品91| 欧美一区二区三区四区高清| 亚洲精选国产| 在线观看视频一区二区| 国产精品一区二区视频 | 欧美影片第一页| 99精品视频一区| 欧美高清视频在线 | 午夜亚洲福利| 一区二区三区欧美日韩| 曰本成人黄色| 国产中文一区| 国产精品综合视频| 欧美日韩激情小视频| 免费一级欧美在线大片| 欧美在线视频观看| 亚洲欧美99| 中文精品视频一区二区在线观看| 欧美成人一二三| 久久综合伊人77777麻豆| 性做久久久久久久免费看| 一区二区成人精品 | 国产精品日韩在线播放| 欧美日韩在线三区| 欧美日韩国产区| 欧美精品一区三区| 欧美啪啪成人vr| 欧美激情精品久久久久久变态| 麻豆精品在线播放| 免费在线观看成人av| 美女脱光内衣内裤视频久久网站| 亚洲欧美日韩中文播放| 亚洲欧美日韩成人| 午夜电影亚洲| 久久狠狠婷婷| 美女福利精品视频| 欧美1区2区3区| 欧美激情一区二区三区高清视频| 欧美xxx成人| 欧美日韩久久| 欧美三级在线| 国产女人水真多18毛片18精品视频| 国产精品毛片va一区二区三区 | 久久亚洲精品中文字幕冲田杏梨| 久久精品91| 欧美sm视频| 欧美日韩一级黄| 国产精品一区二区久久久| 国产欧美一区二区三区在线老狼 | 9久草视频在线视频精品| 亚洲丝袜av一区| 久久精品2019中文字幕| 美日韩精品视频| 亚洲毛片在线观看.| 亚洲永久免费精品| 久久精品国产成人| 欧美精品久久99久久在免费线| 欧美精品午夜| 国产欧美一级| 亚洲日本国产| 欧美一区成人| 欧美国产精品中文字幕| 一本色道**综合亚洲精品蜜桃冫| 亚洲在线成人| 欧美jizz19hd性欧美| 国产精品福利影院| 一区二区在线视频播放| 在线综合亚洲欧美在线视频| 亚洲欧美激情一区二区| 欧美88av| 亚洲一区影音先锋| 鲁大师成人一区二区三区| 国产精品国产三级国产a| 黄色成人在线观看| 亚洲一区二区少妇| 欧美3dxxxxhd| 亚洲免费视频成人| 欧美福利网址| 国模精品一区二区三区色天香| 妖精视频成人观看www| 久久―日本道色综合久久| 亚洲精品免费一二三区| 久久精品夜色噜噜亚洲aⅴ| 欧美性淫爽ww久久久久无| 又紧又大又爽精品一区二区| 午夜精品久久久久久久| 亚洲欧洲一区二区三区久久| 性欧美超级视频| 欧美性理论片在线观看片免费| 在线免费观看一区二区三区| 欧美一区免费| 亚洲专区在线| 国产精品av一区二区| 亚洲人成毛片在线播放| 美国十次成人|