BZOJ5299 [Cqoi2018]解锁屏幕 【状压dp】
2018-05-15 15:36
627 查看
题目链接
题解
就一个毒瘤卡常题。。写了那么久
设\(f[i][s]\)表示选了集合\(s\)中的点,最后一个是\(i\),进行转移
要先预处理出两点间的点,然后卡卡常就可以过了
#include<algorithm> #include<iostream> #include<cstring> #include<cstdio> #include<cmath> #include<queue> #include<map> #define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt) #define REP(i,n) for (register int i = 1; i <= (n); i++) #define mp(a,b) make_pair<int,int>(a,b) #define cls(s) memset(s,0,sizeof(s)) #define cp pair<int,int> #define LL long long int #define lbt(x) (x & -x) #define res register using namespace std; const int maxn = 23,maxm = (1 << 21),INF = 1000000000,P = 100000007; inline int read(){ int out = 0,flag = 1; char c = getchar(); while (c < 48 || c > 57){if (c == '-') flag = -1; c = getchar();} while (c >= 48 && c <= 57){out = (out << 3) + (out << 1) + c - 48; c = getchar();} return out * flag; } int X[maxn],Y[maxn]; inline void add(int& x,int y){ x += y; if (x >= P) x -= P; } int n,bin[maxn],f[maxn][maxm],vis[maxn][maxm],a[maxn][maxn]; inline bool onit(int k,int i,int j){ if (X[k] < min(X[i],X[j]) || X[k] > max(X[i],X[j]) || Y[k] < min(Y[i],Y[j]) || Y[k] > max(Y[i],Y[j])) return false; return (X[j] - X[i]) * (Y[k] - Y[i]) - (Y[j] - Y[i]) * (X[k] - X[i]) == 0; } inline bool is4(int s){ int cnt = 0; while (s) cnt++,s -= lbt(s); return cnt >= 4; } inline void init(){ for (res int i = 0; i < n; i++) for (res int j = 0; j < n; j++) if (j != i) for (res int k = 0; k < n; k++) if (k != i && k != j && onit(k,i,j)) a[i][j] |= bin[k]; } inline void solve(){ int ans = 0,maxv = (1 << n) - 1; for (int i = 0; i < n; i++) f[i][bin[i]] = 1; for (res int s = 0; s <= maxv; s++) for (res int u = 0; u < n; u++) if (f[u][s]){ for (int j = 0; j < n; j++) if (!(s & bin[j])){ if ((s & a[u][j]) != a[u][j]) continue; add(f[j][s | bin[j]],f[u][s]); } } for (res int s = 0; s <= maxv; s++) if (is4(s)) for (int i = 0; i < n; i++) add(ans,f[i][s]); printf("%d\n",ans); } int main(){ bin[0] = 1; for (res int i = 1; i < maxn; i++) bin[i] = bin[i - 1] << 1; n = read(); for (res int i = 0; i < n; i++) X[i] = read(),Y[i] = read(); init(); solve(); return 0; }
相关文章推荐
- bzoj5299: [Cqoi2018]解锁屏幕
- [CQOI 2018]解锁屏幕
- [CQOI2018] 解锁屏幕
- BZOJ_5301_[Cqoi2018]异或序列&&CF617E_莫队
- BZOJ_5296_[Cqoi2018]破解D-H协议_BSGS
- BZOJ 2669: [cqoi2012]局部极小值 状压dp 容斥原理
- bzoj 5297: [Cqoi2018]社交网络
- bzoj 5301: [Cqoi2018]异或序列
- bzoj 5298: [Cqoi2018]交错序列
- BZOJ5296 [CQOI2018] 破解D-H协议 【数学】【BSGS】
- BZOJ5300 [Cqoi2018]九连环 【数学】【FFT】
- BZOJ 3930: [CQOI2015]选数 递推
- BZOJ3930 【CQOI2015】选数
- bzoj5195 [Usaco2018 Feb]Directory Traversal(树+递推)
- Android唤醒、解锁屏幕代码实例
- bzoj 4814: [Cqoi2017]小Q的草稿【计算几何】
- BZOJ3503: [Cqoi2014]和谐矩阵
- [BZOJ5286][HNOI2018]转盘(线段树)
- BZOJ 2668: [cqoi2012]交换棋子
- bzoj 5306: [Haoi2018]染色