uva 11008

uva 11008

狀態壓縮+記憶化搜索

#include<stdio.h>
#include<string.h>
#define max 20
#define maxn 2<<18
#define INF 1<<30
int min(int a,int b)
{
    return a>b?b:a;
}
int g[max][max],f[maxn],limit,x[max],y[max];
int m,n;
int dp(int tree)
{

    if(f[tree]!=-1)
        return f[tree];
    int i,j,k,l,t;
    t=0;
    for(i=0;i<n;i++)
    {
        if((1<<i)&tree)
            t++;
    }
    if(t<=limit)
    {
        f[tree]=0;
        return 0;
    }
    if(t==1)
    {
        f[tree]=1;
        return 1;
    }
    else
    {
        f[tree]=INF;
        int ans;
        for(i=0;i<n;i++)
        {
            if( (1<<i)&tree )
                for(j=0;j<n;j++)
                {
					if(i==j)
						continue;
                    if( (1<< j)&tree && g[i][j] )
                    ans=dp(tree&(~g[i][j]))+1,
                    f[tree]=min(f[tree],ans);
                }
        }
        return f[tree];
    }
}
int main()
{
    int i,j,k,l,t;
    scanf("%d",&l);
    for(t=1;t<=l;t++)
    {
        scanf("%d %d",&n,&m);
        limit=n-m;
        for(i=0;i<n;i++)
        {
            scanf("%d %d",&x[i],&y[i]);
        }
        memset(g,0,sizeof(g));
        memset(f,-1,sizeof(f));
        for(i=0;i<n;i++)       //任意兩棵樹間連線中的樹的狀態
        {
            for(j=0;j<n;j++)
            {
                if(i==j)
                    continue;
                for(k=0;k<n;k++)
                {
                    if((y[i]-y[j])*(x[k]-x[j])==(y[k]-y[j])*(x[i]-x[j]))
                        g[i][j]=g[i][j]|(1<<k);
                }
            }
        }
        printf("Case #%d:\n%d\n",t,dp((1<<n)-1));
        if(t!=l)
            printf("\n");
    }
    return 0;
}