- 1、本文档共13页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
Dancing Links, FFT 模板
Dancing Links
#include iostream
#include cstring
#include cstdio
using namespace std;
typedef long long ll;
templateint N,int M
struct DLX {
int n, m, size;
int U[M], D[M], L[M], R[M], Row[M], Col[M];
int H[M], S[M];
int ansd, ans[N];
int Limit;
void init(int _n,int _m,int lim) {
n = _n; m = _m; Limit = lim;
for(int i = 0; i = m; ++i) {
S[i] = 0;
U[i] = D[i] = i;
L[i] = i - 1;
R[i] = i + 1;
}
R[m] = 0; L[0] = m;
size = m;
for(int i = 1; i = n; ++i)
H[i] = -1;
ansd = 0;
}
void add(int r,int c) {
++S[Col[++size]=c];
Row[size]=r;
D[size]=D[c];
U[D[c]]=size;
U[size]=c;
D[c]=size;
if(H[r]0)H[r]=L[size]=R[size]=size;
else {
R[size]=R[H[r]];
L[R[H[r]]]=size;
R[H[r]]=size;
L[size]=H[r];
}
}
void remove(int c) {
L[R[c]]=L[c],R[L[c]]=R[c];
for(int i = D[c]; i!=c; i=D[i]) {
for(int j = R[i]; j!=i; j=R[j])
U[D[j]]=U[j],D[U[j]]=D[j],--S[Col[j]];
}
}
void resume(int c) {
for(int i = U[c]; i!=c; i=U[i]) {
for(int j = R[i]; j!=i; j=R[j])
U[D[j]]=D[U[j]]=j,++S[Col[j]];
}
L[R[c]]=R[L[c]]=c;
}
void select(int r) {
for(int i = R[H[r]]; i != H[r]; i = R[i])
remove (Col[i]);
remove(Col[H[r]]);
ans[ansd++] = r;
}
bool main() {
return Dance(ansd);
}
bool Dance(int d) {
if(d Limit) return false;
if(R[0] == 0) { ansd = d; return true; }
int c = R[0];
for(in
文档评论(0)