poj 2288 Islands and Bridges——状压dp(哈密尔顿回路)

题目：http://poj.org/problem?id=2288

不知为什么记忆化搜索就是WA得不得了！

#include<iostream>
#include<cstdio>
#include<cstring>
#define ll long long
using namespace std;
const int N=,Lm=(<<)+;
int q,n,m,lm,v[N];
ll ans,prn,dp[Lm][N][N],f[Lm][N][N];
bool b[N][N];
ll dfs(int S,int x,int px)
{
    if(dp[S][x][px]!=-)return dp[S][x][px];
    if(S==lm){f[S][x][px]=;return dp[S][x][px]=;}
    for(int i=;i<=n;i++)
        if(b[x][i]&&(S&(<<(i-)))==)
        {
            int D=(S|(<<(i-)));
            ll w=dfs(D,i,x)+v[i]+v[x]*v[i]+(b[px][i]?v[px]*v[x]*v[i]:);
            if(w>dp[S][x][px])
                dp[S][x][px]=w,f[S][x][px]=f[D][i][x];
            else if(w==dp[S][x][px])
                f[S][x][px]+=f[D][i][x];
        }
    return dp[S][x][px];
}
int main()
{
    scanf("%d",&q);
    while(q--)
    {
        memset(b,,sizeof b);
        scanf("%d%d",&n,&m);lm=(<<n)-;
        for(int i=;i<=n;i++)scanf("%d",&v[i]);
        int x,y;
        for(int i=;i<=m;i++)
        {
            scanf("%d%d",&x,&y);b[x][y]=;b[y][x]=;
        }
        if(n==){printf("%d 1\n",v[]);continue;}//
        memset(dp,-,sizeof dp);memset(f,,sizeof f);
        ans=-;prn=;
        for(int i=;i<=n;i++)
            for(int j=;j<=n;j++)
                if(b[i][j])
                {
                    int D=((<<(i-))|(<<(j-)));//注意括号！
                    ll w=dfs(D,i,j)+v[i]*v[j]+v[i]+v[j];
                    if(w>ans)ans=w,prn=f[D][i][j];
                    else if(w==ans)prn+=f[D][i][j];
                }
        printf("%I64d %I64d\n",prn?ans:,prn>>);//
    }
    return ;
}

改成刷表就可A了。注意初始化一遍，然后开始dp，而不是一边初始化一边dp很多次。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int N=,Lm=(<<)+;
int T,n,m,c[N],lm;
ll dp[Lm][N][N],f[Lm][N][N],ans,prn;
bool b[N][N];
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        memset(b,,sizeof b);
        scanf("%d%d",&n,&m);lm=(<<n)-;
        for(int i=;i<=n;i++)scanf("%d",&c[i]);
        int x,y;
        while(m--)
        {
            scanf("%d%d",&x,&y);
            b[x][y]=b[y][x]=;
        }
        if(n==){printf("%d 1\n",c[]);continue;}
        memset(dp,-,sizeof dp);memset(f,,sizeof f);
        for(int u=;u<=n;u++) for(int v=;v<=n;v++) if(b[u][v])
        {
            int s=((<<(u-))|(<<(v-)));
            dp[s][u][v]=c[u]+c[v]+c[u]*c[v];
            f[s][u][v]=;
        }
        for(int s=;s<=lm;s++)
            for(int i=;i<=n;i++) if(s&(<<(i-)))
                for(int j=;j<=n;j++) if(dp[s][i][j]!=-)
                    for(int k=;k<=n;k++) if(b[i][k]&&!(s&(<<(k-))))
                    {
                        int d=(s|(<<(k-)));
                        ll w=dp[s][i][j]+c[k]+c[i]*c[k]+(b[j][k]?c[i]*c[j]*c[k]:);
                        if(dp[d][k][i]<w)
                        {
                            dp[d][k][i]=w;
                            f[d][k][i]=f[s][i][j];
                        }
                        else if(dp[d][k][i]==w)
                            f[d][k][i]+=f[s][i][j];
                    }
        ans=-;prn=;//因为有"prn?"
        for(int i=;i<=n;i++) for(int j=;j<=n;j++) if(dp[lm][i][j]!=-)
            if(dp[lm][i][j]>ans)ans=dp[lm][i][j],prn=f[lm][i][j];
            else if(dp[lm][i][j]==ans)prn+=f[lm][i][j];
        printf("%I64d %I64d\n",prn?ans:,prn>>);
    }
    return ;
}