【JZOJ 3383】【NOIP2013模拟】太鼓达人（欧拉图？）

问题描述
七夕祭上，Vani牵着cl的手，在明亮的灯光和欢乐的气氛中愉快地穿行。这时，在前面忽然出现了一台太鼓达人机台，而在机台前坐着的是刚刚被精英队伍成员XLk、Poet_shy和lydrainbowcat拯救出来的的applepi。看到两人对太鼓达人产生了兴趣，applepi果断闪人，于是cl拿起鼓棒准备挑战。然而即使是在普通难度下，cl的路人本性也充分地暴露了出来。一曲终了，不但没有过关，就连鼓都不灵了。Vani十分过意不去，决定帮助工作人员修鼓。

鼓的主要元件是M个围成一圈的传感器。每个传感器都有开和关两种工作状态，分别用1和0表示。显然，从不同的位置出发沿顺时针方向连续检查K个传感器可以得到M个长度为K的01串。Vani知道这M个01串应该是互不相同的。而且鼓的设计很精密，M会取到可能的最大值。现在Vani已经了解到了K的值，他希望你求出M的值，并给出字典序最小的传感器排布方案。
输入
一个整数K。
输出
一个整数M和一个二进制串，由一个空格分隔。表示可能的最大的M，以及字典序最小的排布方案，字符0表示关，1表示开。你输出的串的第一个字和最后一个字是相邻的。
样例输入
3
样例输出
8 00010111

样例解释：

得到的8个01串分别是000、001、010、101、011、111、110和100。注意前后是相邻的。长度为3的二进制串总共只有8种，所以M = 8一定是可能的最大值。
算法讨论
正解貌似是什么欧拉图？嘛先贴着好了。

很显然，第一问的答案就是 2^n。
第二问，构造一个有 2^(n-1)个节点的图，对应 2^(n-1)个 n-1 位二进制数。从代表数 k 的节
点(0<=k<2^(n-1))向代表数(k<<1)&(1<<(n-1))的节点，和代表数((k<<1)+1)&(1<<(n-1))的节点
分别连一条边。可以发现这样的图中，所有点的入度和出度都是 2，因此这是一个欧拉图。
因此我们从 0 号点开始 dfs 寻找一个欧拉回路，回溯的时候记录到栈中，最后倒序输出即可。
因为要求字典序最小，dfs 的时候要注意搜索顺序，先走 0 边，再走 1 边。这个算法寻找欧
拉回路，每个点、每条边只访问一次，是 O(V+E)级别的。
时间复杂度 O(2^n)，预计得分 100 分。

第一问就是个组合。第二问我就只是用了个暴力，显然位于高位的1越少越好，那么我们就从前往后枚举，能放0就放0，然后把那些二进制数转成十进制标记一下就好。

#include <cstdio>
#define MAX_N 2500
using namespace std;
int a[MAX_N],k,n;
bool f[MAX_N],flag;

int mul(int x)
{
    int s=1;
    for (int i=k;i>=k-x+1;i--)
        s*=i;
    return s;
}

int mul1(int x)
{
    int s=1;
    for (int i=2;i<=x;i++)
        s*=i;
    return s;
}

bool check(int x)
{
    int p=1,s=0;
    for (int i=x;i>=x-k+1;i--)
    {
        s+=a[i]*p;
        p*=2;
    }
    if (!f[s])
    {
        f[s]=1;
        return 1;
    }
    return 0;
}

bool check1()
{
    int p=1,s=0;
    for (int i=k-1;i>=1;i--)
    {
        s+=a[i]*p;
        p*=2;
    }
    s+=a[n]*p*2;
    if (!f[s])
        return 1;
    return 0;
}

void dfs(int dep)
{
    if (dep>n)
    {
        if (check1())
            flag=1;
        return;
    }
    a[dep]=0;
    if (check(dep))
        dfs(dep+1);
    if (!flag)
    {
        a[dep]=1;
        if (check(dep))
            dfs(dep+1);
    }
}

int main()
{
    scanf("%d",&k);
    n=k+2; a[1]=0; f[0]=1;
    for (int i=k-1;i>=2;i--)
    {
        n+=mul(i) / mul1(i);
        a[i]=0;
    }
    dfs(k+1);
    printf("%d ",n);
    for (int i=1;i<=n;i++)
        printf("%d",a[i]);
}