hdu 3909 数独扩展

思路：做法与9*9的一样。只不过是变量。

#include<set>
#include<map>
#include<cmath>
#include<queue>
#include<cstdio>
#include<vector>
#include<string>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define pb push_back
#define mp make_pair
#define Maxn 4100
#define Maxm 200010
#define Y 1100
#define LL __int645rssss
#define Abs(x) ((x)>0?(x):(-x))
#define lson(x) (x<<1)
#define rson(x) (x<<1|1)
#define inf 100000
#define lowbit(x) (x&(-x))
#define clr(x,y) memset(x,y,sizeof(x))
#define Mod 1000000007
using namespace std;
int H[Maxn],R[Maxn*Y],L[Maxn*Y],D[Maxn*Y],U[Maxn*Y],Q[Maxn*Y],X[Maxn*Y],C[Maxn*Y],S[Maxn*Y];
int cnt,n,g[][],id,N,num;
void init(int n)
{
    int i;
    for(i=;i<=n;i++){
        R[i]=i+;
        L[i+]=i;
        U[i]=D[i]=i;
        S[i]=;
    }
    clr(H,-);
    R[n]=;
    id=n+;
}
void ins(int r,int c)
{
    U[id]=c;
    D[id]=D[c];
    U[D[c]]=id;
    D[c]=id;
    S[c]++;
    if(H[r]<)
        H[r]=R[id]=L[id]=id;
    else {
        R[id]=R[H[r]];
        L[id]=H[r];
        L[R[H[r]]]=id;
        R[H[r]]=id;
    }
    C[id]=c;
    X[id++]=r;
}
void Remove(int c)
{
    R[L[c]]=R[c];
    L[R[c]]=L[c];
    for(int i=D[c];i!=c;i=D[i]){
        for(int j=R[i];j!=i;j=R[j]){
            D[U[j]]=D[j];
            U[D[j]]=U[j];
            S[C[j]]--;
        }
    }
}
void Resume(int c)
{
    R[L[c]]=c;
    L[R[c]]=c;
    for(int i=D[c];i!=c;i=D[i]){
        for(int j=R[i];j!=i;j=R[j]){
            U[D[j]]=j;
            D[U[j]]=j;
            S[C[j]]++;
        }
    }
}
bool dfs(int step,int f)
{
    int i,k,c,j,temp;
    temp=inf;
    if(R[]==) {
        num=step;
        if(cnt||f){
        cnt++;
        return true;
        }
        cnt++;
        return false;
    }
    for(i=R[];i;i=R[i]) if(S[i]<temp){
        temp=S[i];
        c=i;
    }
    Remove(c);
    for(i=D[c];i!=c;i=D[i]){
        Q[step]=X[i];
        for(j=R[i];j!=i;j=R[j]){
            Remove(C[j]);
        }
        if(dfs(step+,f))
            return true;
        for(j=L[i];j!=i;j=L[j]){
            Resume(C[j]);
        }
    }
    Resume(c);
    return false;
}
void build()
{
    int i,j,k;
    int r,c;
    init(N*N*);
    for(i=;i<=N;i++){
        for(j=;j<=N;j++){
            if(g[i][j]){
                r=(i-)*N*N+(j-)*N+g[i][j];
                c=(i-)*N+g[i][j];
                ins(r,c);
                c=N*N+(j-)*N+g[i][j];
                ins(r,c);
                c=*N*N+(i-)*N+j;
                ins(r,c);
                c=*N*N+((i-)/n*n+(j+n-)/n-)*N+g[i][j];
                ins(r,c);
            }
            else{
                for(k=;k<=N;k++){
                r=(i-)*N*N+(j-)*N+k;
                c=(i-)*N+k;
                ins(r,c);
                c=N*N+(j-)*N+k;
                ins(r,c);
                c=*N*N+(i-)*N+j;
                ins(r,c);
                c=*N*N+((i-)/n*n+(j+n-)/n-)*N+k;
                ins(r,c);
                }
            }
        }
    }
}
void solve()
{
    int i,j,k,r,c;
    cnt=;
    build();
    dfs(,);
    if(!cnt) {
        printf("No Solution\n");
        return;
    }
    if(cnt==){
        printf("Multiple Solutions\n");
        return;
    }
    for(i=;i<=N;i++){
        for(j=;j<=N;j++){
            if(!g[i][j]) continue;;
            c=g[i][j];
            g[i][j]=;
            cnt=;
            build();
            dfs(,);
            if(cnt<)
                break;
            g[i][j]=c;
        }
        if(j<=N) break;
    }
    if(i<=N){
        printf("Not Minimal\n");
        return;
    }
    cnt=;
    build();
    dfs(,);
    for(i=;i<num;i++){
        r=(Q[i]-)/(N*N)+;
        c=(Q[i]-(r-)*(N*N)-)/N+;
        k=(Q[i]-)%N+;
        g[r][c]=k;
    }
    for(i=;i<=N;i++){
        for(j=;j<=N;j++){
            if(g[i][j]<)
                printf("%d",g[i][j]);
            else{
                printf("%c",'A'+g[i][j]-);
            }
        }
        printf("\n");
    }
}
int main()
{
    int i,j,k;
    char str[];
    while(scanf("%d",&n)!=EOF){
        clr(g,);
        N=n*n;
        for(i=;i<=N;i++){
            scanf("%s",str);
            for(j=;j<N;j++) if(str[j]!='.'){
                if(str[j]<=''&&str[j]>='')
                g[i][j+]=str[j]-'';
                else
                g[i][j+]=str[j]-'A'+;
            }
        }
        solve();
    }
    return ;
}