POJ 3925 - 상태 DP.비트 연산

2554 단어 OO연상
북대의 를 연구하여 이 문제를 발견한...DFS 열거점에 대해 설명합니다..그리고 최소 생성 트리로 답을 찾아서...마침 얼마 전에 상태 DP 질문을 했는데...그냥 상태 DP 물로 넘어가서...
한 종류의 점 개수가 n=15 정도인 문제에 대해...상태 DP를 민감하게 연상해야 하는데...n비트 2진수로 비교적 좋은 모든 점의 상태에서...이 문제는 바로 이렇다..x(0<=x<=2^n)로 현재 나무에 어떤 점이 있는지 표시...
여기 두 자리 연산이 들어갔는데...
하나는 십진수 x가 이진수 아래의 k위가 1인지 아닌지를 판단하는 것이다.x&(1<또 하나는 십진수 x가 이진수 아래에서 몇 자리가 1인지 판단하는 것이다.x=x & (x-1) x==0까지 실행하는 횟수는 x중 1의 개수...이유..왜냐하면 항상 맨 끝에 있는 1을 빼요.
사실 쓰고 나서...본질적으로 매거점을 최소화해서 트리를 만드는 게 저랑 같은 상태인 DP랑 똑같아요...왜냐하면 제가 업데이트 과정에서 프라임의 과정이었기 때문에...
Program:
#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<math.h>
#include<map>
#include<queue>
#include<stack>
#define ll long long
#define oo 1000000000
#define pi acos(-1)
using namespace std;    
int n,m,p[20],arc[20][20],g,ans,dp[40000]; 
double minimal;
int main()
{  
     int i,x,y,k,w,z,v;
     while (~scanf("%d%d",&n,&m))
     {
           if (!n && !m) break;
           for (i=1;i<=n;i++) scanf("%d",&p[i]);
           for (y=1;y<=n;y++)
              for (x=1;x<=n;x++) 
                 scanf("%d",&arc[y][x]);
           g=(1<<n)-1;
           memset(dp,-1,sizeof(dp));
           dp[0]=0;
           minimal=1e+20;
           for (x=1;x<=g;x++)
           {
                v=0;
                for (k=1;k<=n;k++)
                   if (x & (1<<k-1)) //    k    1 
                   {
                          v+=p[k];
                          w=x-(1<<k-1);
                          for (y=1;y<=n;y++)
                            if (w & (1<<y-1))
                              if (dp[x]==-1 || dp[x]>dp[w]+arc[y][k])
                                  dp[x]=dp[w]+arc[y][k];
                   }
                k=0;
                z=x;
                while (z)
                {
                      k++;
                      z=z & (z-1);
                } //             .. 
                if (k==m && minimal>(dp[x]*1.0/v))
                {
                        minimal=dp[x]*1.0/v;
                        ans=x;
                }
           }
           for (x=1;x<=n;x++)
              if (ans & (1<<x-1)) printf("%d ",x);
           printf("
"); } return 0; }

좋은 웹페이지 즐겨찾기