[CQOI2018]解锁屏幕

Posted 2021-01-31 asuldb

状压这个东西好像没有什么能优化的高级东西，像什么斜率优化，单调队列在状压的优化上都很少见

而最常见的状压优化就是预处理优化了，

这道题就预处理一下所有点对之间连线上的点，之后压成状态就能做到(O(2^n*n^2))

这道题的状态就非常简单了，就是一个小学生状压(dp[i][S])状态为(S)时最后一个点是(i)的方案数

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#define min(a,b) ((a)<(b)?(a):(b))
#define max(a,b) ((a)>(b)?(a):(b))
#define re register
#define lowbit(x) ((x)&(-x))
#define eps 1e-6
const int mod=100000007;
int n;
int dp[21][1048577];
int can[21][21];
inline int read()
{
    char c=getchar();
    int x=0,r=1;
    while(c<‘0‘||c>‘9‘) 
    {
        if(c==‘-‘) r=-1;
        c=getchar();
    }
    while(c>=‘0‘&&c<=‘9‘)
        x=(x<<3)+(x<<1)+c-48,c=getchar();
    return x*r;
}
inline int check(double x,double y)
{
    if(x+eps>y&&x-eps<y) return 1;
    return 0;
}
int X[21],Y[21];
inline int same(int a,int b,int c)
{
    if(X[c]<min(X[a],X[b])||X[c]>max(X[a],X[b])||Y[c]>max(Y[a],Y[b])||Y[c]<min(Y[a],Y[b])) return 0;
    if(X[a]==X[c]) return X[b]==X[c];
    if(X[b]==X[c]) return X[a]==X[c];
    return (Y[a]-Y[c])*(X[b]-X[c])==(Y[b]-Y[c])*(X[a]-X[c]);
}
inline int cnt(int x)
{
    int tot=0;
    while(x) tot++,x-=lowbit(x);
    return tot;
}
int main()
{
    n=read();
    for(re int i=1;i<=n;i++)
        X[i]=read(),Y[i]=read();
    int N=(1<<n)-1;
    for(re int i=1;i<=n;i++)
        for(re int j=i+1;j<=n;j++)
        {
            can[i][j]=can[j][i]=N;
            for(re int k=1;k<=n;k++)
                if(k!=i&&k!=j&&same(i,j,k)) can[i][j]^=(1<<(k-1)),can[j][i]^=(1<<(k-1));
        }
    for(re int i=1;i<=n;i++)
        dp[i][1<<(i-1)]++;
    for(re int i=1;i<N;i++)
        for(re int j=1;j<=n;j++)
            if(dp[j][i]&&(i&(1<<(j-1))))
            {
                for(re int k=1;k<=n;k++)
                if(!(i&(1<<(k-1)))&&((can[j][k]|i)==N)) 
                    dp[k][i|(1<<(k-1))]=(dp[k][i|(1<<(k-1))]+dp[j][i])%mod;
            }
    int ans=0;
    for(re int i=1;i<=N;i++)
    if(cnt(i)>=4) 
        for(re int j=1;j<=n;j++) 
            ans=(ans+dp[j][i])%mod;
    std::cout<<ans;
    return 0; 
}