HDU - 4332 - 多校四1002 - 矩阵+状压

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <cstdlib>
#include <algorithm>
using namespace std;

__int64 tri[300][300];
int ok[300];
int pos[300];
int sta;

struct node
{
    __int64 x[300][300];
}tmp[35], ans;

void state()//把256*256的矩阵变成128*128的，否则会超时。
{
    sta = 0;
    memset(ok, 0, sizeof(ok));
    int i, j, tmp;
    for(i=0; i<=255; i++)
    {
        tmp = 0;
        for(j=0; j<8; j++)
        if( ((1<<j)&i) ) tmp++;
        if( tmp%2==0 ) ok[i] = 1;
        else ok[i] = 0;
    }
    for(i=0; i<=255; i++)
    if(ok[i])
        pos[i] = sta++;
}

__int64 can(int a, int b)//状态a 转到状态b
{
    int i, now;
    if(a==255) swap(a, b);
    if(a==255) return 2;//状态255 转到状态255 有两种转移
    if(a==0) now=0;

    else //now为第一个1
    {
        for(i=0; i<8; i++)
        {
            int j ;
            if(i==0) j = 7;
            else j = i-1;
            if( ((1<<j)&a)==0 && ((1<<(i))&a) )
            {
                now = i;break;
            }
        }
    }
    int tt;
    for(i=now, tt=0; tt<8; tt++,i=(i+1)%8)
    {
        int tmp = (i+1)%8;
        if( ((1<<i)&a)==0 )
        {
            if( ((1<<i)&b)==0 ) return 0;
        }
        else
        {
            if( ((1<<tmp)&a) && ((1<<i)&b) && ((1<<tmp)&b) ) i++; //横着放
            else if( ((1<<i)&b)==0 ) continue; //竖着放
            else return 0;
        }
    }
    return 1;
}

void make_tri() //初始化最初的矩阵
{
    int i, j, k;
    for(i=0; i<=255; i++)
    for(j=0; j<=255; j++)
    if(ok[i] && ok[j])
    {
        tri[ pos[i] ][ pos[j] ] = can(i, j);
    }
    tri[sta-1][sta-1] = 2;
}

void init()
{
    int i, j;
    for(i=0; i<sta; i++)
    for(j=0; j<sta; j++)
    {
        ans.x[i][j] = (i==j);
    }
}

void mulans(int a)
{
    node c;
    int i, j, k;
    for(i=0; i<sta; i++)
    for(j=0; j<sta; j++)
    {
        c.x[i][j] = 0;
        for(k=0; k<sta; k++)
        {
            c.x[i][j] += ans.x[i][k] * tmp[a].x[k][j];
            c.x[i][j] %= mod;
        }
    }
    for(i=0; i<sta; i++)
    for(j=0; j<sta; j++)
        ans.x[i][j] = c.x[i][j];
}

void multmp(int a, int b)
{
    node c;
    int i, j, k;
    for(i=0; i<sta; i++)
    for(j=0; j<sta; j++)
    {
        tmp[a].x[i][j] = 0;
        for(k=0; k<sta; k++)
        {
            tmp[a].x[i][j] += tmp[b].x[i][k] * tmp[b].x[k][j];
            tmp[a].x[i][j] %= mod;
        }
    }
}

void make_tmp()
{
    int i, j;
    for(i=0; i<sta; i++)
    for(j=0; j<sta; j++)
    {
        tmp[0].x[i][j] = tri[i][j];
    }
    for(i=1; i<33; i++)
    {
        multmp(i, i-1);
    }
}

int main()
{
    int t, tt=0, n;
    state();
    make_tri();
    make_tmp();
    scanf("%d", &t);
    while(t--)
    {
        scanf("%d", &n);
        init();
        for(int i=0; (1<<i)<=n; i++)
        {
            if((1<<i)&n) mulans(i);
        }
        printf("Case %d: %I64d\n", ++tt, ans.x[sta-1][sta-1]);
    }
	return 0;
}