zoukankan      html  css  js  c++  java
  • BZOJ 1072 [SCOI2007]安排perm 如压力DP

    意甲冠军:联系

    方法:状压DP?

    题解:这题事实上没啥好写的。不算非常难,推一推就能搞出来。

    首先看到这个问题,对于被d整除这个条件,非常easy就想到是取余数为0,所以想到可能状态中刚開始含有取余数。

    先说我的第一个想法。f[i][j]表示选取i个数且此时的mod为j,这种思想是第一下蹦出来的,当时想的就是在线来搞终于的答案。只是转瞬即发现,这TM不就是暴力吗魂淡!并没有什么卵用,于是開始想这个状态可不能够做什么优化。

    显然第二维的j并不须要太大的优化。暂且先将其搁置一边。来考虑第一维的i怎么优化。

    把滚动栏拉到最下方发现这个i的范围最大是10,也就是说hash的空间为2^10-1,能够从状态的角度往下走,假设这个i代表hash状态,能不能表示他的转移方程呢?

    首先依照一般的想法,对于每一个i,可能从它的某一为1的位数变为零的状态转移过来,假设依照这个方向搞的话或许也能做出来,只是,我们更熟悉的操作是加上而并非从原状态减去某个位,所以思路又转向到逆向更新。

    也就是说我们能够初步写出来dp方程

    f[(i|(1<<(k1)))][(j10+a[k])%d]+=f[i][j]

    当中更新的时候要保证i&(1<<(k-1))==0

    然后我们做的仅仅须要枚举i,枚举j,枚举k

    嗯哼就是这样。然而这中途当然有非常多自己主动剪掉的部分。所以复杂度不是非常高。题目时间范围可承受。

    可是并没有完,这时候你会发现如此的代码连第一个小例子都过不了,也就是说。假设有反复的数字的话,我们会算非常多反复的情况。而怎么去重呢?仅仅须要将答案除掉反复的数的个数的排列就好了。也就是阶乘!

    后记:事实上写这个题解的目的就是将自己解一道题的全过程写出来。

    代码:

    #include <cstdio>
    #include <cstring>
    #include <iostream>
    #include <algorithm>
    #define N 15
    #define M 1100
    using namespace std;
    int tot,mod,l;
    int a[N];
    int cnt[N];
    char s[N];
    int f[1<<10][M];
    int A[]={0,1,2,6,24,120,720,5040,40320,362880,3628800};
    void dfs()
    {
        for(int i=0;i<(1<<l);i++)
        {
            for(int k=0;k<mod;k++)
            {
                if(f[i][k])
                {
                    for(int j=1;j<=l;j++)
                    {
                        if((i&(1<<(j-1)))==0)
                        {
                            f[(i|(1<<(j-1)))][(k*10+a[j])%mod]+=f[i][k];
                        }
                    }
                }
            }
        }
    }
    int main()
    {
        scanf("%d",&tot);
        while(tot--)
        {
            scanf("%s",s);
            l=strlen(s);
            memset(cnt,0,sizeof(cnt));
            for(int i=0;i<l;i++)
            {
                a[i+1]=s[i]-'0';
                cnt[a[i+1]]++;
            }
            scanf("%d",&mod);
            memset(f,0,sizeof(f));
            f[0][0]=1;
            dfs();
            int ans=f[(1<<l)-1][0];
            for(int i=0;i<=9;i++)
            {
                if(cnt[i])ans/=A[cnt[i]];
            }
            printf("%d
    ",ans);
        }
    }

    版权声明:本文博客原创文章。博客,未经同意,不得转载。

  • 相关阅读:
    linux ramdisk
    linux系统灵活运用灯[android课程3]
    linux编译注解
    Linux内核3.0移植并基于Initramfs根文件系统启动
    linux0.11文件分析
    2.2linux内核移植简介
    sudo: ./sd_fusing.sh:找不到命令
    Tiny4412之C语言实现流水灯,Tiny4412裸机程序[3]
    Exynos 4412的启动过程分析[2]
    POJ 3009 Curling 2.0 {深度优先搜索}
  • 原文地址:https://www.cnblogs.com/hrhguanli/p/4754869.html
Copyright © 2011-2022 走看看