P4876 近似排列计数

最新推荐文章于 2024-04-16 12:06:03 发布

HT008_123

最新推荐文章于 2024-04-16 12:06:03 发布

阅读量331

点赞数 1

分类专栏：题目分析 DP DFS 状态压缩

本文链接：https://blog.csdn.net/qq_35914587/article/details/78161777

版权

题目分析同时被 3 个专栏收录

671 篇文章 4 订阅

订阅专栏

91 篇文章 1 订阅

订阅专栏

DFS

42 篇文章 0 订阅

订阅专栏

本文探讨了一个计数问题的两种解决方案：一种是通过全排列的方法进行暴力判断；另一种则是采用状态压缩动态规划（状压DP）的方式解决。通过这两种方法的对比，可以帮助读者更好地理解如何针对特定的问题选择合适的算法。

摘要由CSDN通过智能技术生成

不用看啦，蒟蒻博主没AC，只打了50分的代码！
1：全排列方法暴力判断

#include <cstdio>
#include <iostream>
#include <cstring> 
using namespace std;
const int mod=1e9+7;
bool vis[21];
int a[21],ans;
void dfs(int x,int p,int k)
{
    if(x==p+1)
    {
        ans=(ans+1)%mod;
        return;
    }
    int l,r;
    l=max(1,x-k),r=min(p,x+k);
    if(a[x]) dfs(x+1,p,k);
    else
    {
      if(x>=3)
      if(!vis[x-k])//关键：假设这一位没有用这个数，后面的一定不可能用这个，直接返回就好！
      {
        vis[max(1,x-k)]=1,dfs(x+1,p,k),vis[max(1,x-k)]=0;
        return;
      }
      for(int i=l;i<=r;i++)
      {
        if(!vis[i]) vis[i]=1,dfs(x+1,p,k),vis[i]=0;
      }
    }
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        memset(a,0,sizeof(a));
        memset(vis,0,sizeof(vis));
        ans=0;
        int n,m,k;
        scanf("%d%d%d",&n,&m,&k);
        for(int i=1;i<=m;i++)
         {
            int x,y;
            scanf("%d%d",&x,&y);
            a[x]=y;
            vis[y]=1;
         }
        if(k==0) ans=1;
        else dfs(1,n,k);
        printf("%d\n",ans); 
    }

}

2：状压DP
用一个20位串表示1-n是否用过，枚举每个位置，看可以上面填的数是否用过转移即可

#include <cstdio>
#include <iostream>
#include <cstring> 
using namespace std;
const int mod=1e9+7;
const int maxm=(1<<20)+1;
int dp[maxm],a[21];
int bx[]={-1,0,1,-2,2};
void work1(int n,int t)
{
    memset(dp,0,sizeof(dp));
    dp[0]=1;
    for(int i=1;i<=n;i++)
     {
        if(a[i])//此处已有数
        {
            for(int k=(1<<n)-1;k>=0;k--)
             if(!(k&(1<<a[i]-1))) //判断有没有用过a[i]这个数
             dp[k|(1<<a[i]-1)]=(dp[k|(1<<a[i]-1)]+dp[k])%mod;//没有向加入a[i]的串中转移方案数
            continue;
        }
        int r;
        if(t==1) r=2;
        else r=4;
        for(int k=(1<<n)-1;k>=0;k--)
         for(int j=0;j<=r;j++) 
          {
            int x=i+bx[j];//同上
            if(x<1||x>n) continue;//不能越界！
            if(!(k&(1<<x-1)))  dp[k|(1<<x-1)]=(dp[k|(1<<x-1)]+dp[k])%mod;
          }
     }
     printf("%d\n",dp[(1<<n)-1]);
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
     int n,m,k;
     memset(a,0,sizeof(a));
     scanf("%d%d%d",&n,&m,&k);
     for(int i=1;i<=m;i++)
     {
        int x,y;
        scanf("%d%d",&x,&y);
        a[x]=y;
     }

     if(k==0) printf("1\n");
     else work1(n,k);
    }
}