C++博客-心如止水-随笔分类-题目分类：搜索

POJ 2531 Network Saboteur

lee1r — Sat, 06 Aug 2011 03:00:00 GMT

2985ms AC。
以下是我的代码：

/*
* Author:  lee1r
* Created Time:  2011/8/6 10:45:18
* File Name: poj2531.cpp
*/
#include<iostream>
#include<sstream>
#include<fstream>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<ctime>
#define L(x) ((x)<<1)
#define R(x) (((x)<<1)+1)
#define Half(x) ((x)>>1)
#define lowbit(x) ((x)&(-(x)))
using namespace std;
const int kInf(0x7f7f7f7f);
const double kEps(1e-8);
typedef long long int64;
typedef unsigned long long uint64;

const int kMaxn(27);

int n,g[kMaxn][kMaxn];
bool d[kMaxn];
int64 ans;

void dfs(int depth)
{
    if(depth>n)
    {
        int now(0);
        for(int i=2;i<=n;i++)
            for(int j=1;j<i;j++)
                if(d[i]==!d[j])
                    now+=g[i][j];
        if(ans<now)
            ans=now;
        return;
    }
    d[depth]=true;
    dfs(depth+1);
    d[depth]=false;
    dfs(depth+1);
}

int main()
{
    //freopen("data.in","r",stdin);

    while(scanf("%d",&n)==1)
    {
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                scanf("%d",&g[i][j]);

        ans=0;
        memset(d,false,sizeof(d));
        dfs(1);

        cout<<ans<<endl;
    }

    return 0;
}

lee1r 2011-08-06 11:00 发表评论

POJ 1321 棋盘问题

lee1r — Sat, 06 Aug 2011 01:45:00 GMT

以下是我的代码：

/*
* Author:  lee1r
* Created Time:  2011/8/6 9:35:45
* File Name: poj1321.cpp
*/
#include<iostream>
#include<sstream>
#include<fstream>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<ctime>
#define L(x) ((x)<<1)
#define R(x) (((x)<<1)+1)
#define Half(x) ((x)>>1)
#define lowbit(x) ((x)&(-(x)))
using namespace std;
const int kInf(0x7f7f7f7f);
const double kEps(1e-8);
typedef long long int64;
typedef unsigned long long uint64;

const int kMaxn(10);

int n,k;
int64 ans;
bool g[kMaxn][kMaxn],row[kMaxn],col[kMaxn];

void dfs(int depth,int now)
{
    if(now==k)
    {
        ans++;
        return;
    }
    if(depth>n)
        return;
    dfs(depth+1,now);
    for(int i=1;i<=n;i++)
        if(!row[i] && !col[i] && g[depth][i])
        {
            row[i]=col[i]=true;
            dfs(depth+1,now+1);
            row[i]=col[i]=false;
        }
}

int main()
{
    //freopen("data.in","r",stdin);

    while(scanf("%d%d",&n,&k)==2 && (n!=-1 || k!=-1))
    {
        for(int i=1;i<=n;i++)
        {
            getchar();
            for(int j=1;j<=n;j++)
            {
                char ch;
                scanf("%c",&ch);
                if(ch=='#')
                    g[i][j]=true;
                else
                    g[i][j]=false;
            }
        }

        ans=0;
        memset(row,false,sizeof(row));
        memset(col,false,sizeof(col));
        dfs(1,0);

        cout<<ans<<endl;
    }

    return 0;
}

lee1r 2011-08-06 09:45 发表评论

POJ 3009 Curling 2.0

lee1r — Sat, 06 Aug 2011 01:34:00 GMT

简单的迭代加深搜索求最优解。注意判断可以继续行走的条件。
以下是我的代码：

/*
* Author:  lee1r
* Created Time:  2011/8/6 8:51:53
* File Name: poj3009.cpp
*/
#include<iostream>
#include<sstream>
#include<fstream>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<ctime>
#define L(x) ((x)<<1)
#define R(x) (((x)<<1)+1)
#define Half(x) ((x)>>1)
#define lowbit(x) ((x)&(-(x)))
using namespace std;
const int kInf(0x7f7f7f7f);
const double kEps(1e-8);
typedef long long int64;
typedef unsigned long long uint64;

const int kMaxn(27);
const int dx[]={1,-1,0,0},dy[]={0,0,-1,1};

int n,m,s1,s2,g1,g2,r[kMaxn][kMaxn];
bool success;

int Check(int x,int y,int dir,int &xx,int &yy)
{
    int t(0);
    do
    {
        if(r[x][y]==3)
        {
            xx=x;
            yy=y;
            return 3;
        }
        else if(t<=1 && r[x][y]==1)
            return 0;
        else if(t>1 && r[x][y]==1)
        {
            xx=x;
            yy=y;
            return 1;
        }
        x+=dx[dir];
        y+=dy[dir];
        t++;
    }while(x>=1 && x<=n && y>=1 && y<=m);
    return 0;
}

void dfs(int depth,int x,int y,int maxdepth)
{
    if(depth>maxdepth)
    {
        if(g1==x && g2==y)
            success=true;
        return;
    }
    for(int i=0;i<4;i++)
    {
        int newx,newy,result(Check(x,y,i,newx,newy));
        if(result==1)
        {
            r[newx][newy]=0;
            dfs(depth+1,newx-dx[i],newy-dy[i],maxdepth);
            r[newx][newy]=1;
        }
        else if(result==3)
            dfs(depth+1,newx,newy,maxdepth);
    }
}

int main()
{
    //freopen("data.in","r",stdin);

    while(scanf("%d%d",&m,&n)==2 && (n || m))
    {
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
            {
                scanf("%d",&r[i][j]);
                if(r[i][j]==2)
                {
                    s1=i;s2=j;
                }
                else if(r[i][j]==3)
                {
                    g1=i;g2=j;
                }
            }

        int ans;
        success=false;
        for(ans=1;ans<=10;ans++)
        {
            dfs(1,s1,s2,ans);
            if(success)
                break;
        }

        if(ans>10)
            printf("%d\n",-1);
        else
            printf("%d\n",ans);
    }

    return 0;
}

lee1r 2011-08-06 09:34 发表评论

POJ 2488 A Knight's Journey

lee1r — Sat, 06 Aug 2011 00:25:00 GMT

经典搜索问题。
以下是我的代码：

/*
* Author:  lee1r
* Created Time:  2011/8/6 7:47:26
* File Name: poj2488.cpp
*/
#include<iostream>
#include<sstream>
#include<fstream>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<ctime>
#define L(x) ((x)<<1)
#define R(x) (((x)<<1)+1)
#define Half(x) ((x)>>1)
#define lowbit(x) ((x)&(-(x)))
using namespace std;
const int kInf(0x7f7f7f7f);
const double kEps(1e-8);
typedef long long int64;
typedef unsigned long long uint64;

const int kMaxn(27);
const int dx[]={-2,-2,-1,1,2,2,1,-1},dy[]={-1,1,2,2,1,-1,-2,-2};

int n,m;
string ans,d;
bool used[kMaxn][kMaxn];

string itoa(int x)
{
    string re,t;
    while(x>0)
    {
        t+=(x%10+'0');
        x/=10;
    }
    for(int i=t.size()-1;i>=0;i--)
        re+=t[i];
    return re;
}

void erase(string &t)
{
    int pos;
    for(pos=t.size()-1;pos>=0;pos--)
        if(isalpha(t[pos]))
            break;
    t.erase(pos,t.size()-pos);
}

void dfs(int depth,int x,int y)
{
    if(ans<=d)
        return;
    if(depth>n*m)
    {
        ans=d;
        return;
    }
    for(int i=0;i<8;i++)
    {
        int newx(x+dx[i]),newy(y+dy[i]);
        if(newx<1 || newx>n || newy<1 || newy>m || used[newx][newy])
            continue;
        d+=(char)(newx+'A'-1);
        d+=itoa(newy);
        used[newx][newy]=true;
        dfs(depth+1,newx,newy);
        erase(d);
        used[newx][newy]=false;
    }
}

int main()
{
    //freopen("data.in","r",stdin);

    int T;
    scanf("%d",&T);
    for(int case_num=1;case_num<=T;case_num++)
    {
        scanf("%d%d",&m,&n);
        ans="z";
        for(int i=1;i<=min(Half(n)+1,n);i++)
            for(int j=1;j<=min(Half(m)+1,m);j++)
            {
                d+=(char)(i+'A'-1);
                d+=itoa(j);
                used[i][j]=true;
                dfs(2,i,j);
                erase(d);
                used[i][j]=false;
            }

        printf("Scenario #%d:\n",case_num);
        if(ans<"z")
            cout<<ans<<endl;
        else
            cout<<"impossible"<<endl;
        cout<<endl;
    }

    return 0;
}

lee1r 2011-08-06 08:25 发表评论

POJ 2965 The Pilots Brothers' refrigerator

lee1r — Sun, 31 Jul 2011 01:39:00 GMT

与POJ 1753类似，但是需要加点优化，这里主要使用了位运算的方法。

#include<cstdio>
#include<cstring>
#define g(x,y) (((x-1)<<2)+y)
using namespace std;

int a,d;
int ansn,ans;

void f(int i,int &x,int &y)
{
    x=i/4+1;
    y=i%4;
    if(i%4==0)
    {
        x--;
        y=4;
    }
}

void handle(int &r,int x,int y)
{
    for(int i=1;i<=4;i++)
    {
        r^=(1<<g(x,i));
        r^=(1<<g(i,y));
    }
    r^=(1<<g(x,y));
}

void dfs(int last,int have)
{
    if(have>=ansn)
        return;

    int r(a);
    for(int i=1;i<=16;i++)
        if(d&(1<<i))
        {
            int x,y;
            f(i,x,y);
            handle(r,x,y);
        }
    if(r==0)
    {
        ansn=have;
        ans=d;
    }

    for(int i=last+1;i<=16;i++)
    {
        d=(d|(1<<i));
        dfs(i,have+1);
        d=(~((~d)|(1<<i)));
    }
}

int main()
{
    /*
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    //*/

    a=0;
    for(int i=1;i<=4;i++)
    {
        for(int j=1;j<=4;j++)
        {
            char t;
            scanf("%c",&t);
            if(t=='+')
                a=(a|(1<<g(i,j)));
            else
                a=(~((~a)|(1<<g(i,j))));
        }
        getchar();
    }

    ansn=0x7f7f7f7f;
    d=0;
    dfs(0,0);

    printf("%d\n",ansn);
    for(int i=1;i<=16;i++)
        if(ans&(1<<i))
        {
            int x,y;
            f(i,x,y);
            printf("%d %d\n",x,y);
        }

    return 0;
}

lee1r 2011-07-31 09:39 发表评论

POJ 1753 Flip Game

lee1r — Sun, 31 Jul 2011 01:36:00 GMT

分析得知最多只可能有16次操作，而且每次操作都不相同。因此枚举每个操作有没有执行即可。
以下是我的代码：

#include<cstdio>
#include<cstring>
using namespace std;
const int dx[]={-1,1,0,0},dy[]={0,0,1,-1};

char a[5][5];
bool d[17];
int ans;

bool Black(char r[5][5])
{
    for(int i=1;i<=4;i++)
        for(int j=1;j<=4;j++)
            if(r[i][j]=='w')
                return false;
    return true;
}
bool White(char r[5][5])
{
    for(int i=1;i<=4;i++)
        for(int j=1;j<=4;j++)
            if(r[i][j]=='b')
                return false;
    return true;
}

void f(int i,int &x,int &y)
{
    x=i/4+1;
    y=i%4;
    if(i%4==0)
    {
        x--;
        y=4;
    }
}

void handle(char r[5][5],int x,int y)
{
    for(int i=0;i<4;i++)
    {
        int newx(x+dx[i]),newy(y+dy[i]);
        if(newx<1 || newx>4 || newy<1 || newy>4)
            continue;
        if(r[newx][newy]=='b')
            r[newx][newy]='w';
        else
            r[newx][newy]='b';
    }
    if(r[x][y]=='b')
        r[x][y]='w';
    else
        r[x][y]='b';
}

void dfs(int last,int have)
{
    char r[5][5];
    memcpy(r,a,25*sizeof(char));
    for(int i=1;i<=16;i++)
        if(d[i])
        {
            int x,y;
            f(i,x,y);
            handle(r,x,y);
        }
    if(Black(r) || White(r))
    {
        if(ans>have)
            ans=have;
    }

    for(int i=last+1;i<=16;i++)
    {
        d[i]=true;
        dfs(i,have+1);
        d[i]=false;
    }
}

int main()
{
    /*
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    //*/

    for(int i=1;i<=4;i++)
    {
        for(int j=1;j<=4;j++)
            scanf("%c",&a[i][j]);
        getchar();
    }

    ans=0x7f7f7f7f;
    memset(d,false,sizeof(d));
    dfs(0,0);

    if(ans<0x7f7f7f7f)
        printf("%d\n",ans);
    else
        printf("Impossible\n");

    return 0;
}

lee1r 2011-07-31 09:36 发表评论

NOIP 2009 TG 4 靶型数独

lee1r — Wed, 20 Jul 2011 12:02:00 GMT

#include<cstdio>
using namespace std;
const int score[10][10]={{0,0,0,0,0,0,0,0,0,0},
                         {0,6,6,6,6,6,6,6,6,6},
                         {0,6,7,7,7,7,7,7,7,6},
                         {0,6,7,8,8,8,8,8,7,6},
                         {0,6,7,8,9,9,9,8,7,6},
                         {0,6,7,8,9,10,9,8,7,6},
                         {0,6,7,8,9,9,9,8,7,6},
                         {0,6,7,8,8,8,8,8,7,6},
                         {0,6,7,7,7,7,7,7,7,6},
                         {0,6,6,6,6,6,6,6,6,6}};

int r[10][10],ans(-1),tot(0);
bool squ[4][4][10],row[10][10],col[10][10];

void f(int i,int &x,int &y)
{
    x=i/9+1;
    y=i%9;
    if(i%9==0)
    {
        x--;
        y=9;
    }
}
void g(int x,int y,int &xx,int &yy)
{
    xx=x/3;
    if(x%3)
        xx++;
    yy=y/3;
    if(y%3)
        yy++;
}

void dfs(const int &depth,const int &now)
{
    tot++;
    if(tot>=13000000)
        return;

    if(depth<=0)
    {
        if(ans<now)
            ans=now;
        return;
    }

    int x,y,xx,yy;
    f(depth,x,y);
    g(x,y,xx,yy);
    if(r[x][y])
        dfs(depth-1,now+r[x][y]*score[x][y]);
    else
    {
        for(int i=1;i<=9;i++)
            if(!squ[xx][yy][i] && !row[x][i] && !col[y][i])
            {
                squ[xx][yy][i]=row[x][i]=col[y][i]=true;
                dfs(depth-1,now+i*score[x][y]);
                squ[xx][yy][i]=row[x][i]=col[y][i]=false;
            }
    }
}

int main()
{
    //*
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    //*/

    for(int i=1;i<=9;i++)
        for(int j=1;j<=9;j++)
        {
            scanf("%d",&r[i][j]);
            int xx,yy;
            g(i,j,xx,yy);
            squ[xx][yy][r[i][j]]=row[i][r[i][j]]=col[j][r[i][j]]=true;
        }

    dfs(81,0);
    printf("%d\n",ans);

    return 0;
}

lee1r 2011-07-20 20:02 发表评论

POJ 1270 Following Orders

lee1r — Fri, 27 May 2011 10:40:00 GMT

题目大意：给出变量列表和约束关系，计算出所有的可能的大小关系。
比较标准的做法应该构造一个图（这样效率比较高）。不过数据太弱了，直接全排列所有的可能，然后判断是否符合就AC了。
以下是我的代码：

#include<iostream>
#include<string>
#include<algorithm>
#include<cstdio>
#include<cctype>
using namespace std;
const int kMaxn(107);

int main()
{
    /*
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    //*/

    int T(0);
    string line1,line2;
    while(getline(cin,line1) && getline(cin,line2))
    {
        T++;
        if(T!=1)
            cout<<endl;

        string var;
        for(int i=0;i<line1.size();i++)
            if(isalpha(line1[i]))
                var+=line1[i];
        sort(var.begin(),var.end());

        int cnt(0);
        char ch1,ch2,x[kMaxn],y[kMaxn];
        for(int i=0,j=0;i<line2.size();i++)
        {
            if(isalpha(line2[i]))
            {
                j++;
                if(j&1)
                    ch1=line2[i];
                else
                {
                    ch2=line2[i];
                    cnt++;
                    x[cnt]=ch1;
                    y[cnt]=ch2;
                }
            }
        }

        do
        {
            bool success(true);
            for(int i=1;i<=cnt;i++)
                if(find(var.begin(),var.end(),y[i])-find(var.begin(),var.end(),x[i])<0)
                {
                    success=false;
                    break;
                }
            if(success)
                cout<<var<<endl;
        }while(next_permutation(var.begin(),var.end()));
    }

    return 0;
}

lee1r 2011-05-27 18:40 发表评论

UVa 993 Product of digits

lee1r — Tue, 24 May 2011 03:04:00 GMT

题目大意：给出一个数字，找到一个最小的数字，使得这个数字各位数的乘积等于给定的数字。
网上许多人的做法我不晓得如何证明做法是正确的。
由于一个数字的因子个数不会很多，采用DFS即可。
以下是我的代码：

#include<vector>
#include<cstdio>
using namespace std;

int len;
vector<int> ans,now;

bool update()
{
    if(len>now.size())
        return true;
    for(int i=0;i<len;i++)
        if(ans[i]>now[i])
            return true;
        else if(ans[i]<now[i])
            return false;
    return false;
}

void dfs(int depth,int x,int last)
{
    if(x==1)
    {
        if(update())
        {
            ans=now;
            len=depth-1;
        }
        return;
    }
    if(depth>len)
        return;
    for(int i=last;i<=9;i++)
        if(x%i==0)
        {
            now.push_back(i);
            dfs(depth+1,x/i,i);
            now.pop_back();
        }
}

int main()
{
    #ifndef ONLINE_JUDGE
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    #endif

    int T;
    scanf("%d",&T);
    while(T--)
    {
        int n;
        scanf("%d",&n);

        if(n==0 || n==1)
        {
            printf("%d\n",n);
            continue;
        }
        len=0x7f7f7f7f;
        now.clear();
        dfs(1,n,2);

        if(len==0x7f7f7f7f)
            printf("%d\n",-1);
        else
        {
            for(int i=0;i<ans.size();i++)
                printf("%d",ans[i]);
            printf("\n");
        }
    }

    return 0;
}

lee1r 2011-05-24 11:04 发表评论

经典迭代加深搜索——埃及分数

lee1r — Mon, 09 May 2011 14:59:00 GMT

给出一个分数，比如19/45，把它写成若干个形如1/Ri的分数的和的形式，比如19/45=1/5+1/6+1/18，要求分母不能重复使用并且使用的分数的个数最少。（如果有多组个数相同的解，最后的分数的分母越小越好，这对于题目来说是次要的。）
1、分母从小到大搜索
为了避免重复搜索
2、使用迭代加深搜索
求“步骤数最少”这类问题，基本上有两种似乎：广搜、迭代加深搜索。对于这道题来说，如果广搜将永远得不到结果，分母可以无限大！但是迭代加深搜索就比较好，虽然做了许多重复工作，但状态空间至少被限制住了。如果当前正在枚举的分母，使得接下来的选择即使每次都选择最大，达到最大深度的时候也不可能达到目标分数，那么当前正在枚举的分母及比它还大的分母，都不需要枚举了。这样可以给分母确定一个上界。另外，已经得到的结果加上当前枚举的分母对应的分数，要小于等于目标分数，这样给分母确定了一个下界（可以在O(1)的复杂度内确定这个下界）。

在下面的代码中，因为确定上下界都使用了浮点运算，最终还是需要把当前结果和目标结果比较。
浮点运算～真让人纠结的东西！

以下是我的代码：

#include<algorithm>
#include<cstdio>
#include<cmath>
using namespace std;

int gcd(int a,int b)
{
    for(int t=a%b;t;a=b,b=t,t=a%b); return b;
}
struct Type
{
    Type():a_(0),b_(1) {}
    Type(int a,int b):a_(a),b_(b) {}

    int a_,b_;
};
Type operator+(const Type &a,const Type &b)
{
    Type re(a.a_*b.b_+a.b_*b.a_,a.b_*b.b_);
    int t(gcd(re.a_,re.b_));
    re.a_/=t;
    re.b_/=t;
    return re;
}
bool operator==(const Type &a,const Type &b)
{
    return (a.a_*b.b_==b.a_*a.b_);
}

Type target;
int ans,r[10000],t[10000];
bool found;

void dfs(const int &depth,const int &last,const Type &now)
{
    if(depth>ans)
    {
        if(now==target)
        {
            if(found)
            {
                bool cover(false);
                for(int i=ans;i>=1;i--)
                    if(r[i]>t[i])
                    {
                        cover=true;
                        break;
                    }
                    else if(r[i]<t[i])
                        break;
                if(cover)
                    for(int i=1;i<=ans;i++)
                        r[i]=t[i];
            }
            else
            {
                for(int i=1;i<=ans;i++)
                    r[i]=t[i];
                found=true;
            }
        }
        return;
    }
    for(int i=max(last+1,(int)ceil((double)(target.b_*now.b_)/(target.a_*now.b_-target.b_*now.a_))); ;i++)
    {
        Type news(now+Type(1,i));
        if(depth+(int)ceil((double)(target.a_*news.b_-target.b_*news.a_)*(i+1)/(target.b_*news.b_))>ans)
            break;
        t[depth]=i;
        dfs(depth+1,i,news);
    }
}

int main()
{
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);

    while(scanf("%d%d",&target.a_,&target.b_)==2)
    {
        found=false;
        for(ans=1; ;ans++)
        {
            dfs(1,0,Type());
            if(found)
                break;
        }

        for(int i=1;i<=ans;i++)
        {
            if(i!=1)
                printf(" ");
            printf("%d",r[i]);
        }
        printf("\n");
    }

    return 0;
}

lee1r 2011-05-09 22:59 发表评论

UVa 10422 Knights in FEN

lee1r — Thu, 21 Apr 2011 14:04:00 GMT

题目大意是：在5×5的棋盘上有12个白子和12个黑子，另外有一个空格，白子或黑子可以按照骑士的行走方式走到空格位置。给定初始状态和目标状态，判断是否能够在10步（包含10步）以内达到目标状态。
典型的BFS，但是迭代加深的DFS效率同样很高，而且编程复杂度很低。在这里有一个优化：比如当前所做的搜索最大深度设定为maxdepth，比较当前状态和目标状态有多少处不同，结果假设为diff，则如果当前深度depth+diff/2>maxdepth，则停止当前搜索。比如有2处不同，那么最好情况就是2/2=1次转移到目标状态。
以下是我的代码：

#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;
const char kTarget[5][5]={ {'1','1','1','1','1'},
                           {'0','1','1','1','1'},
                           {'0','0',' ','1','1'},
                           {'0','0','0','0','1'},
                           {'0','0','0','0','0'} };
const int dx[]={-2,-2,2,2,-1,1,-1,1},dy[]={-1,1,-1,1,-2,-2,2,2};

char r[5][5];
int ans;
bool success;

int Diff()
{
    int re(0);
    for(int i=0;i<5;i++)
        for(int j=0;j<5;j++)
            re+=(r[i][j]!=kTarget[i][j]);
    return (re/2);
}

void dfs(int depth,int x,int y)
{
    if(depth>=ans)
    {
        if(memcmp(r,kTarget,25*sizeof(char))==0)
            success=true;
        return;
    }
    //*
    if(depth+Diff()>ans)
        return;
    //*/
    for(int i=0;i<8 && !success;i++)
    {
        int newx(x+dx[i]),newy(y+dy[i]);
        if(newx<0 || newx>= 5 || newy<0 || newy>=5)
            continue;
        swap(r[x][y],r[newx][newy]);
        dfs(depth+1,newx,newy);
        swap(r[x][y],r[newx][newy]);
    }
}

int main()
{
    #ifndef ONLINE_JUDGE
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    #endif

    int T;
    scanf("%d",&T);
    getchar();
    while(T--)
    {
        int x,y;
        for(int i=0;i<5;i++)
        {
            for(int j=0;j<5;j++)
            {
                scanf("%c",&r[i][j]);
                if(r[i][j]==' ')
                {
                    x=i;
                    y=j;
                }
            }
            getchar();
        }

        success=false;
        for(ans=0;ans<=10;ans++)
        {
            dfs(0,x,y);
            if(success)
                break;
        }

        if(success)
            printf("Solvable in %d move(s).\n",ans);
        else
            printf("Unsolvable in less than 11 move(s).\n");
    }

    return 0;
}

lee1r 2011-04-21 22:04 发表评论

UVa 301 Transportation

lee1r — Thu, 21 Apr 2011 01:44:00 GMT

生成子集问题，每个订单选或不选的问题。
一开始心想22个元素的子集规模并不大，420万左右，于是使用了先生成子集，再检测判断的方法，结果TLE了；
然后采用边生成子集，边检测是否符合条件的搜索，然后AC之～
以下是我的代码：

#include<algorithm>
#include<cstdio>
using namespace std;
const int kMaxk(27);

int n,m,k,ans,b[kMaxk],e[kMaxk],p[kMaxk],cnt[kMaxk];
bool used[kMaxk];

void dfs(int s,int now)
{
    ans=max(ans,now);
    for(int i=s+1;i<=k;i++)
    {
        used[i]=true;
        bool success(true);
        for(int j=b[i];j<e[i];j++)
        {
            cnt[j]+=p[i];
            if(cnt[j]>n)
                success=false;
        }
        if(success)
            dfs(i,now+(e[i]-b[i])*p[i]);
        used[i]=false;
        for(int j=b[i];j<e[i];j++)
            cnt[j]-=p[i];
    }
}

int main()
{
    #ifndef ONLINE_JUDGE
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    #endif

    while(scanf("%d%d%d",&n,&m,&k)==3 && (n || m || k))
    {
        for(int i=1;i<=k;i++)
            scanf("%d%d%d",&b[i],&e[i],&p[i]);

        ans=0;
        dfs(0,0);

        printf("%d\n",ans);
    }

    return 0;
}

lee1r 2011-04-21 09:44 发表评论

UVa 167 The Sultan's Successors

lee1r — Wed, 20 Apr 2011 12:22:00 GMT

八皇后问题，只不过每一个点上加了一个权值，求一个权值最大的解。
以下是我的代码：

#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;

int ans,r[8][8],d[8];
bool used1[8],used2[20],used3[20];

void dfs(int depth)
{
    if(depth>=8)
    {
        int t(0);
        for(int i=0;i<8;i++)
            t+=r[i][d[i]];
        ans=max(ans,t);
        return;
    }
    for(int i=0;i<8;i++)
        if(!used1[i] && !used2[depth+i] && !used3[8+depth-i-1])
        {
            d[depth]=i;
            used1[i]=used2[depth+i]=used3[8+depth-i-1]=true;
            dfs(depth+1);
            used1[i]=used2[depth+i]=used3[8+depth-i-1]=false;
        }
}

int main()
{
    #ifndef ONLINE_JUDGE
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    #endif

    int T;
    scanf("%d",&T);
    while(T--)
    {
        for(int i=0;i<8;i++)
            for(int j=0;j<8;j++)
                scanf("%d",&r[i][j]);

        ans=0;
        memset(d,0,8*sizeof(int));
        memset(used1,false,8*sizeof(bool));
        memset(used2,false,20*sizeof(bool));
        memset(used3,false,20*sizeof(bool));
        dfs(0);

        printf("%5d\n",ans);
    }

    return 0;
}

lee1r 2011-04-20 20:22 发表评论

UVa 331 Mapping the Swaps

lee1r — Wed, 20 Apr 2011 05:27:00 GMT

求将给定序列排成升序的最小相邻交换的方案总数。
用BFS编程复杂度较高，需要队列、判重……由于数据规模很小，使用迭代加深的DFS也可。
以下是我的代码：

#include<algorithm>
#include<cstdio>
using namespace std;

int n,ans,r[7];
bool success;

bool OK()
{
    for(int i=1;i<=n-1;i++)
        if(r[i]>r[i+1])
            return false;
    return true;
}

void dfs(int depth,int maxdepth)
{
    if(depth>maxdepth)
    {
        if(OK())
        {
            success=true;
            if(depth>1)
                ans++;
        }
        return;
    }
    for(int i=1;i<=n-1;i++)
    {
        swap(r[i],r[i+1]);
        dfs(depth+1,maxdepth);
        swap(r[i],r[i+1]);
    }
}

int main()
{
    #ifndef ONLINE_JUDGE
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    #endif

    int T(0);
    while(scanf("%d",&n)==1 && n)
    {
        T++;
        for(int i=1;i<=n;i++)
            scanf("%d",&r[i]);

        ans=0;
        success=false;
        for(int i=0; ;i++)
        {
            dfs(1,i);
            if(success)
                break;
        }

        printf("There are %d swap maps for input data set %d.\n",ans,T);
    }

    return 0;
}

lee1r 2011-04-20 13:27 发表评论

UVa 539 The Settlers of Catan

lee1r — Wed, 20 Apr 2011 02:39:00 GMT

果搜即可。
从一个起点开始出发，尝试每条没有走过的边，直到不能走为止。
注意到和顶点数相比，边数要少得多，所以使用邻接表存储。
以下是我的代码：

#include<list>
#include<cstdio>
#include<cstring>
using namespace std;
const int kMaxn(27);

int n,m,ans;
list<int> g[kMaxn];
bool used[kMaxn][kMaxn];

void dfs(int node,int length)
{
    if(ans<length)
        ans=length;
    for(list<int>::iterator i=g[node].begin();i!=g[node].end();i++)
        if(!used[node][*i])
        {
            used[node][*i]=used[*i][node]=true;
            dfs(*i,length+1);
            used[node][*i]=used[*i][node]=false;
        }
}

int main()
{
    #ifndef ONLINE_JUDGE
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    #endif

    while(scanf("%d%d",&n,&m)==2 && (n || m))
    {
        for(int i=0;i<n;i++)
            g[i].clear();
        for(int i=0;i<m;i++)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            g[u].push_back(v);
            g[v].push_back(u);
        }

        ans=0;
        for(int i=0;i<n;i++)
        {
            memset(used,false,kMaxn*kMaxn*sizeof(bool));
            dfs(i,0);
        }

        printf("%d\n",ans);
    }

    return 0;
}

lee1r 2011-04-20 10:39 发表评论

UVa 140 Bandwidth

lee1r — Tue, 19 Apr 2011 13:00:00 GMT

生成全排列，然后从中选择最优解即可。
输入输出被我弄得有点麻烦了……
以下是我的代码：

#include<iostream>
#include<string>
#include<algorithm>
#include<cctype>
using namespace std;
const int kMaxn(10);
const int kInf(930315);

int n,ans,r[kMaxn],v[kMaxn];
bool map[kMaxn][kMaxn],used[kMaxn];

void dfs(int depth)
{
    if(depth>n)
    {
        int now(0);
        for(int i=1;i<=n;i++)
            for(int j=n;j>i;j--)
                if(map[r[i]][r[j]])
                {
                    now=max(now,j-i);
                    break;
                }

        if(ans>now)
        {
            ans=now;
            for(int i=1;i<=n;i++)
                v[i]=r[i];
        }
        return;
    }
    for(int i=1;i<=n;i++)
        if(!used[i])
        {
            r[depth]=i;
            used[i]=true;
            dfs(depth+1);
            r[depth]=0;
            used[i]=false;
        }
}

int main()
{
    string graph;
    while(cin>>graph && graph!="#")
    {
        bool t[27]={false};
        for(int i=0;i<graph.size();i++)
            if(isalpha(graph[i]))
                t[graph[i]-'A']=true;
        n=0;
        for(int i=0;i<26;i++)
            if(t[i])
                n++;
        for(int i=0;i<graph.size();i++)
            if(isalpha(graph[i]))
            {
                int cnt(0);
                for(int j=0;j<=graph[i]-'A';j++)
                    if(t[j])
                        cnt++;
                graph[i]=cnt+'0';
            }

        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                map[i][j]=false;

        for(int i=0;i<graph.size();i++)
        {
            int k(graph[i]-'0');
            i+=2;
            while(i<graph.size() && graph[i]!=';')
            {
                map[k][graph[i]-'0']=true;
                map[graph[i]-'0'][k]=true;
                i++;
            }
        }

        ans=kInf;
        fill(used,used+kMaxn,false);
        dfs(1);

        for(int i=1;i<=n;i++)
        {
            int j,cnt(0);
            for(j=0;j<26;j++)
                if(t[j])
                {
                    cnt++;
                    if(cnt>=v[i])
                        break;
                }

            cout<<(char)(j+'A')<<" ";
        }
        cout<<"-> "<<ans<<endl;
    }

    return 0;
}

lee1r 2011-04-19 21:00 发表评论

UVa 11218 KTV

lee1r — Tue, 19 Apr 2011 07:31:00 GMT

采用很暴力的办法即可。
以下是我的代码：

#include<algorithm>
#include<cstdio>
using namespace std;
const int kMaxn(87);

struct Type
{
    int a_,b_,c_,s_;
};

int main()
{
    /*
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    //*/

    int T(0),n;
    while(scanf("%d",&n)==1 && n)
    {
        Type r[kMaxn];
        for(int i=1;i<=n;i++)
            scanf("%d%d%d%d",&r[i].a_,&r[i].b_,&r[i].c_,&r[i].s_);

        int ans(-1);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                for(int k=1;k<=n;k++)
                {
                    bool used[10]={false};
                    used[r[i].a_]=true;
                    used[r[i].b_]=true;
                    used[r[i].c_]=true;
                    used[r[j].a_]=true;
                    used[r[j].b_]=true;
                    used[r[j].c_]=true;
                    used[r[k].a_]=true;
                    used[r[k].b_]=true;
                    used[r[k].c_]=true;

                    bool success(true);
                    for(int l=1;l<=9;l++)
                        if(!used[l])
                        {
                            success=false;
                            break;
                        }

                    if(!success)
                        continue;

                    ans=max(ans,r[i].s_+r[j].s_+r[k].s_);
                }

        T++;
        printf("Case %d: %d\n",T,ans);
    }

    return 0;
}

lee1r 2011-04-19 15:31 发表评论

UVa 639 Don't Get Rooked

lee1r — Mon, 18 Apr 2011 07:47:00 GMT

数据规模很小，直接搜索即可，每次试着在可以放置的位置上放置。
以下是我的代码：

#include<iostream>
#include<cstdio>
using namespace std;
const int kMaxn(7);

int n,ans,map[kMaxn][kMaxn];

void getxy(int depth,int &x,int &y)
{
    x=depth/n;
    y=depth%n;
    if(y)
        x++;
    if(!y)
        y=n;
}

bool check(int x,int y)
{
    for(int i=y-1;i>=1 && map[x][i]!=2;i--)
        if(map[x][i]==1)
            return false;
    for(int i=y+1;i<=n && map[x][i]!=2;i++)
        if(map[x][i]==1)
            return false;
    for(int i=x-1;i>=1 && map[i][y]!=2;i--)
        if(map[i][y]==1)
            return false;
    for(int i=x+1;i<=n && map[i][y]!=2;i++)
        if(map[i][y]==1)
            return false;
    return true;
}

void dfs(int depth,int now)
{
    if(depth>n*n)
    {
        if(ans<now)
            ans=now;
        return;
    }

    int x,y;
    getxy(depth,x,y);

    dfs(depth+1,now);
    if(!map[x][y] && check(x,y))
    {
        map[x][y]=1;
        dfs(depth+1,now+1);
        map[x][y]=0;
    }

}

int main()
{
    /*
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    //*/

    while(cin>>n && n)
    {
        cin.get();
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
            {
                char t;
                cin>>t;
                if(t=='.')
                    map[i][j]=0;
                else
                    map[i][j]=2;
            }

        ans=0;
        dfs(1,0);

        cout<<ans<<endl;
    }

    return 0;
}

lee1r 2011-04-18 15:47 发表评论

UVa 11205 The broken pedometer

lee1r — Fri, 15 Apr 2011 09:29:00 GMT

枚举含p个元素的集合的子集，求出该子集和给定的若干个集合的交集，如果得到的这些交集都不相同，则该子集符合条件。题目要求求出符合条件的子集中最少含有多少个元素。
以下是我的代码：

#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;
const int kMaxp(17);
const int kMaxn(107);

bool cmp(const bitset<kMaxp> &a,const bitset<kMaxp> &b)
{
    return (a.to_ulong()<b.to_ulong());
}

int p,n,ans;
bitset<kMaxp> r[kMaxn];
bitset<kMaxp> used,t[kMaxn];

bool OK()
{
    for(int i=1;i<n;i++)
        if(t[i]==t[i+1])
            return false;
    return true;
}

void dfs(int depth)
{
    if(depth>p)
    {
        for(int i=1;i<=n;i++)
            t[i]=r[i];
        for(int i=1;i<=n;i++)
            t[i]=t[i]&used;
        sort(t+1,t+n+1,cmp);
        if(OK())
            ans=min(ans,(int)used.count());
        return;
    }
    dfs(depth+1);
    used[depth]=1;
    dfs(depth+1);
    used[depth]=0;
}

int main()
{
    /*
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    //*/

    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&p,&n);
        for(int i=1;i<=n;i++)
            r[i].reset();
        for(int i=1;i<=n;i++)
            for(int j=1;j<=p;j++)
            {
                int t;
                scanf("%d",&t);
                r[i][j]=t;
            }

        ans=315;
        used.reset();
        dfs(1);

        printf("%d\n",ans);
    }

    return 0;
}

lee1r 2011-04-15 17:29 发表评论

HDU 1181 变形课

lee1r — Mon, 07 Mar 2011 12:57:00 GMT

构图，BFS。
以下是我的代码：

#include<iostream>
#include<string>
#include<queue>
#include<cstdio>
#include<cstring>
using namespace std;
const string kEndSign ("0");
const int kBeginNum ('b'-'a'+1),kEndNum ('m'-'a'+1);

int main()
{
    /*
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    //*/

    string s;
    bool map[27][27]={false};
    while(cin>>s)
    {
        if(s==kEndSign)
        {
            bool visited[27]={false},success(false);
            queue<int> q;
            q.push(kBeginNum);
            visited[kBeginNum]=true;
            while(!q.empty())
            {
                int now=q.front();q.pop();
                if(now==kEndNum)
                {
                    success=true;
                    break;
                }
                for(int i=1;i<=26;i++)
                    if(!visited[i] && map[now][i])
                    {
                        q.push(i);
                        visited[i]=true;
                    }
            }
            if(success)
                cout<<"Yes."<<endl;
            else cout<<"No."<<endl;

            memset(map,false,27*27*sizeof(bool));
        }
        else
            map[s[0]-'a'+1][s[s.size()-1]-'a'+1]=true;
    }

return 0;
}

lee1r 2011-03-07 20:57 发表评论

HDU 1175 连连看

lee1r — Tue, 01 Mar 2011 03:29:00 GMT

主要在于判断转折了几次。
以下是我的代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<ctime>
using namespace std;
const int kMaxn = 1007;
const int dx[]={-1,0,1,0},dy[]={0,1,0,-1};

int n,m,map[kMaxn][kMaxn];
bool used[kMaxn][kMaxn];
int x1,y1,x2,y2;

bool dfs(int nowx,int nowy,int lastx,int lasty,int turn_times)
{
    if(turn_times>2)
        return false;
    if(nowx==x2 && nowy==y2)
        return true;
    used[nowx][nowy]=true;
    for(int i=0;i<4;i++)
    {
        int newx(nowx+dx[i]),newy(nowy+dy[i]);
        if(newx<1 || newx>n || newy<1 || newy>m)
            continue;
        if(used[newx][newy])
            continue;
        if(map[newx][newy] && (newx!=x2 || newy!=y2))
            continue;
        if((lastx==nowx && nowx==newx) || (lasty==nowy && nowy==newy))
        {
            if(dfs(newx,newy,nowx,nowy,turn_times))
                return true;
        }
        else
        {
            if(dfs(newx,newy,nowx,nowy,turn_times+1))
                return true;
        }
    }
    used[nowx][nowy]=false;
    return false;
}

int main()
{
    /*
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    //*/

    while(cin>>n>>m && (n || m))
    {
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
                cin>>map[i][j];
        int question_num;
        cin>>question_num;
        while(question_num--)
        {
            cin>>x1>>y1>>x2>>y2;
            if(!map[x1][y1] || !map[x2][y2] || map[x1][y1]!=map[x2][y2])
            {
                cout<<"NO"<<endl;
                continue;
            }
            memset(used,false,kMaxn*kMaxn*sizeof(bool));
            if(dfs(x1,y1,x1,y1,0))
                cout<<"YES"<<endl;
            else cout<<"NO"<<endl;
        }
    }
    /*
    cout<(clock())/CLOCKS_PER_SEC<    //*/

return 0;
}

lee1r 2011-03-01 11:29 发表评论

UVa 861 Little Bishops

lee1r — Fri, 12 Nov 2010 14:04:00 GMT

我写的搜索程序太慢了……没敢提交上去……写了个交表的程序。猥琐地飘过……
以下是我的代码：

#include<iostream>
using namespace std;

const long cheat[9][67]={{0},
                         {1,1},
                         {1,4,4},
                         {1,9,26,26,8},
                         {1,16,92,232,260,112,16},
                         {1,25,240,1124,2728,3368,1960,440,32},
                         {1,36,520,3896,16428,39680,53744,38368,12944,1600,64},
                         {1,49,994,10894,70792,282248,692320,1022320,867328,389312,81184,5792,128},
                         {1,64,1736,26192,242856,1444928,5599888,14082528,22522960,22057472,12448832,3672448,489536,20224,256}};

int main()
{
    long n,k;

    while(cin>>n>>k)
    {
       if(n==0&&k==0) break;

       cout<<cheat[n][k]<<endl;
    }
return 0;
}

lee1r 2010-11-12 22:04 发表评论

tyvj P1080 N皇后

lee1r — Fri, 22 Oct 2010 13:35:00 GMT

lie[i]==false表示第i列还没有任何棋子；dui1[i]==false表示向右倾斜的、从上到下的第i条对角线上还没有任何棋子；dui2[i]==false表示向左倾斜的、从上到下的第i条对角线上还没有任何棋子。
利用这些记录优化搜索。
n==13的数据我用了403ms。
以下是我的代码：

#include<iostream>
#include<string.h>
#define maxn 17
using namespace std;

long n,cnt,ans[maxn];
bool lie[maxn],dui1[2*maxn],dui2[2*maxn];

long f1(long x,long y)
{
    return (x-y+n);
}
long f2(long x,long y)
{
    return (x+y-1);
}

void dfs(long dep)
{
    if(dep>n)
    {
        cnt++;
        if(cnt>3) return;
        for(long i=1;i<=n;i++)
        {
            if(i!=1) cout<<" ";
            cout<<ans[i];
        }
        cout<<endl;
        return;
    }
    for(long i=1;i<=n;i++)
        if(!lie[i]&&!dui1[f1(dep,i)]&&!dui2[f2(dep,i)])
        {
            lie[i]=dui1[f1(dep,i)]=dui2[f2(dep,i)]=true;
            ans[dep]=i;
            dfs(dep+1);
            lie[i]=dui1[f1(dep,i)]=dui2[f2(dep,i)]=false;
        }
}

int main()
{
    cin>>n;

    cnt=0;
    memset(lie,false,sizeof(lie));
    memset(dui1,false,sizeof(dui1));
    memset(dui2,false,sizeof(dui2));
    dfs(1);
    cout<<cnt<<endl;

return 0;
}

lee1r 2010-10-22 21:35 发表评论

tyvj P1030 乳草的入侵

lee1r — Fri, 22 Oct 2010 03:29:00 GMT

Easy BFS。
Here is my code:

#include<iostream>
#define maxn 107
using namespace std;
const long xd[]={-1,-1,0,1,1,1,0,-1},yd[]={0,-1,-1,-1,0,1,1,1};

typedef struct
{
    long counter,front,rear,x[maxn*maxn],y[maxn*maxn];
}queue;
void Clear(queue &q)
{
    q.counter=0;q.front=0;q.rear=-1;
}
bool Empty(queue &q)
{
    return (q.counter==0);
}
void Push(queue &q,long rx,long ry)
{
    q.counter++;q.rear++;
    if(q.rear>=maxn*maxn) q.rear=0;
    q.x[q.rear]=rx;q.y[q.rear]=ry;
}
void Pop(queue &q,long &rx,long &ry)
{
    q.counter--;rx=q.x[q.front];ry=q.y[q.front];
    q.front++;
    if(q.front>=maxn*maxn) q.front=0;
}

long n,m,mx,my,ans,g[maxn][maxn],d[maxn][maxn];
queue q;

int main()
{
    /*
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    //*/
    cin>>n>>m>>mx>>my;
    for(long i=1;i<=m;i++)
    {
        char ch;
        for(long j=1;j<=n;j++)
        {
            cin>>ch;
            if(ch=='.') g[i][j]=1;
            else g[i][j]=0;
        }
    }
    long t=my;my=mx;mx=m-t+1;
    //  Input
    /*
    for(long i=1;i<=m;i++)
    {
        for(long j=1;j<=n;j++)
            cout<        cout<    }
    cout<    //*/

    Clear(q);Push(q,mx,my);
    d[mx][my]=0;g[mx][my]=0;
    while(!Empty(q))
    {
        long xx,yy,nx,ny;Pop(q,xx,yy);
        for(long k=0;k<8;k++)
        {
            nx=xx+xd[k];
            ny=yy+yd[k];
            if(nx>=1&&nx<=m&&ny>=1&&ny<=n&&g[nx][ny])
            {
                Push(q,nx,ny);
                d[nx][ny]=d[xx][yy]+1;g[nx][ny]=0;
            }
        }
    }

    ans=0;
    for(long i=1;i<=m;i++)
        for(long j=1;j<=n;j++)
            if(ans<d[i][j])
                ans=d[i][j];
    cout<<ans<<endl;

return 0;
}

lee1r 2010-10-22 11:29 发表评论

rqnoj 34 紧急救援

lee1r — Fri, 15 Oct 2010 14:43:00 GMT

BFS。找到终点之后可以立即跳出循环，提高效率，但是因为数据规模不大，我没有这么做，见谅~~
以下是我的代码：

#include<iostream>
#define maxn 1001
#define INF 20000007
using namespace std;
const long xd[]={-1,1,0,0},yd[]={0,0,1,-1};

typedef struct
{
    long counter,front,rear,x[maxn*maxn],y[maxn*maxn];
}queue;
void Clear(queue &q)
{
    q.counter=0;q.front=0;q.rear=-1;
}
bool Empty(queue &q)
{
    return (q.counter==0);
}
void Push(queue &q,long rx,long ry)
{
    q.counter++;q.rear++;
    q.x[q.rear]=rx;q.y[q.rear]=ry;
}
void Pop(queue &q,long &rx,long &ry)
{
    q.counter--;rx=q.x[q.front];ry=q.y[q.front];
    q.front++;
}

long n,bx,by,ex,ey,d[maxn][maxn];
bool g[maxn][maxn];
queue q;

int main()
{
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);

    cin>>n;
    for(long i=1;i<=n;i++)
    {
        for(long j=1;j<=n;j++)
        {
            char r;
            cin>>r;
            g[i][j]=(r=='0'?false:true);
            d[i][j]=INF;
        }
        getchar();
    }
    cin>>bx>>by;
    cin>>ex>>ey;
    //  Input

    Clear(q);
    d[bx][by]=0;
    Push(q,bx,by);
    while(!Empty(q))
    {
        long xx,yy;Pop(q,xx,yy);
        for(long k=0;k<4;k++)
        {
            long rx=xx+xd[k],ry=yy+yd[k];
            if(rx>=1&&rx<=n&&ry>=1&&ry<=n&&!g[rx][ry]&&d[rx][ry]>=INF)
            {
                d[rx][ry]=d[xx][yy]+1;
                Push(q,rx,ry);
            }
        }
    }

    cout<<d[ex][ey]<<endl;

return 0;
}

lee1r 2010-10-15 22:43 发表评论

rqnoj 70 八数码难题

lee1r — Wed, 22 Sep 2010 02:26:00 GMT

Bfs + Hash is ok.
I found that the functions named mem... are veru useful.
Here is my code:

#include<iostream>
#include<string.h>
#define maxn 1000000
using namespace std;
const int xd[]={-1,1,0,0},yd[]={0,0,-1,1};
int st[maxn][3][3],front=1,rear=1,aim[3][3]={{1,2,3},{8,0,4},{7,6,5}},dist[maxn];
int head[maxn],next[maxn];

long hash(long x)
{
    long re=0;
    for(long i=0;i<3;i++)
        for(long j=0;j<3;j++)
            re=(re*10+st[x][i][j])%maxn;
    return re;
}
void insert(long x)
{
    long h=hash(x),p=head[h];
    while(p)
    {
        if(memcmp(st[p],st[x],sizeof(st[p]))==0) return;
        p=next[p];
    }
    next[x]=head[h];head[h]=x;
    rear++;dist[rear]=dist[front]+1;
}
long bfs()
{
    while(front<=rear)
    {
        if(memcmp(st[front],aim,sizeof(aim))==0) return front;
        //  Found answer , return the position of the final state
        long x,y;
        bool find_zero=false;
        for(x=0;x<3&&!find_zero;x++)
            for(y=0;y<3&&!find_zero;y++)
                if(st[front][x][y]==0)
                    find_zero=true;
        x--;y--;
        //  Find the position of figure zero
        for(long k=0;k<4;k++)
        {
            long newx=x+xd[k],newy=y+yd[k];
            if(newx>=0&&newx<3&&newy>=0&&newy<3)
            {
                memcpy(st[rear+1],st[front],sizeof(st[front]));
                st[rear+1][x][y]=st[rear+1][newx][newy];
                st[rear+1][newx][newy]=0;
                //  Produce new states
                insert(rear+1);
                //  Try to insert
            }
        }
        front++;
    }
    return -1;
}
int main()
{
    string s;
    cin>>s;
    for(long i=0;i<3;i++)
        for(long j=0;j<3;j++)
            st[1][i][j]=s.at(3*i+j)-'0';
    dist[1]=0;
    //  Input
    long ans=bfs();
    //  Bfs
    cout<<dist[ans]<<endl;
    //  Output
return 0;
}

lee1r 2010-09-22 10:26 发表评论

Ural 1005 Stone pile

lee1r — Sun, 05 Sep 2010 14:27:00 GMT

刚看完题目想Greedy，试着举反例，没找到。但是看数据规模像是搜索，还是特简单的那种DFS，发现效率还挺高！提交，AC。
后来看到有DP的做法。原问题等价于从n个物品中选取若干个，其重量不超过sum/2，且重量达到最大。只要令权值等于重量即可。
以下是我的代码，DFS程序：

#include<iostream>
#include<stdlib.h>
#define maxn 27
#define INF 20000007
using namespace std;
long n,sum,ans,w[maxn];
void dfs(long dep,long now)
{
    if(dep>n)
    {
        if(ans>abs(now-(sum-now))) ans=abs(now-(sum-now));
        return;
    }
    dfs(dep+1,now);
    dfs(dep+1,now+w[dep]);
}
int main()
{
    cin>>n;
    sum=0;
    for(long i=1;i<=n;i++)
    {
        cin>>w[i];
        sum+=w[i];
    }
    //  Input
    ans=INF;
    dfs(1,0);
    cout<<ans<<endl;
return 0;
}

lee1r 2010-09-05 22:27 发表评论

POJ 2245 Lotto

lee1r — Thu, 08 Jul 2010 13:30:00 GMT

Summer camp还有上机测试，很长时间没有coding了，准备一下，做了几道题目。
以下是我的代码：

#include<iostream>
using namespace std;
long n,a[17],r[7];
void dfs(long dep,long s)
{
    if(dep>6)
    {
        for(long i=1;i<=6;i++)
        {
            if(i>1) cout<<" ";
            cout<<r[i];
        }
        cout<<endl;
        return;
    }
    for(long i=s+1;i<=n;i++)
    {
        r[dep]=a[i];
        dfs(dep+1,i);
    }
}
int main()
{
    /*
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
    //*/
    bool first=true;
    while(cin>>n&&n!=0)
    {
        if(first) first=false;
        else cout<<endl;
        for(long i=1;i<=n;i++)
            cin>>a[i];
        dfs(1,0);
    }
return 0;
}

lee1r 2010-07-08 21:30 发表评论

vijos P1016 北京2008的挂钟

lee1r — Wed, 07 Apr 2010 13:47:00 GMT

我的第一次……
第一次使用九维数组……
简单地BFS即可。
以下是我的代码：

#include<stdio.h>
#include<string.h>
#define max_size 300000
const long change[  ][9]={{1,1,0,1,1,0,0,0,0},
                          {1,1,1,0,0,0,0,0,0},
                          {0,1,1,0,1,1,0,0,0},
                          {1,0,0,1,0,0,1,0,0},
                          {0,1,0,1,1,1,0,1,0},
                          {0,0,1,0,0,1,0,0,1},
                          {0,0,0,1,1,0,1,1,0},
                          {0,0,0,0,0,0,1,1,1},
                          {0,0,0,0,1,1,0,1,1}};
struct
{
    long count,front,rear,item[max_size][9];
}q;
void Clear()
{
    q.count=0;q.front=0;q.rear=-1;
}
void Push(long *x)
{
    q.count++;
    q.rear++;
    for(long i=0;i<9;i++)
      q.item[q.rear][i]=x[i];
}
void Pop(long *x)
{
    q.count--;
    for(long i=0;i<9;i++)
      x[i]=q.item[q.front][i];
    q.front++;
}
bool Empty()
{
    return (q.count==0);
}

int d[4][4][4][4][4][4][4][4][4];
int f[4][4][4][4][4][4][4][4][4];
long r[17];

void init()
{
    for(long i=0;i<9;i++) scanf("%ld",&r[i]);
    memset(d,-1,sizeof(d));
    memset(f,-1,sizeof(f));
}
bool ok(long *x)
{
    for(long i=0;i<9;i++)
      if(x[i]) return false;
    return true;
}
void bfs()
{
    long t[17],s[17];
    Clear();
    d[r[0]][r[1]][r[2]][r[3]][r[4]][r[5]][r[6]][r[7]][r[8]]=0;
    Push(r);
    while(!Empty())
    {
       Pop(t);
       if(ok(t)) break;

       for(long i=0;i<9;i++)
       {
          for(long j=0;j<9;j++) s[j]=(t[j]+change[i][j])%4;
          if(d[s[0]][s[1]][s[2]][s[3]][s[4]][s[5]][s[6]][s[7]][s[8]]==-1)
          {
             d[s[0]][s[1]][s[2]][s[3]][s[4]][s[5]][s[6]][s[7]][s[8]]=d[t[0]][t[1]][t[2]][t[3]][t[4]][t[5]][t[6]][t[7]][t[8]]+1;
             f[s[0]][s[1]][s[2]][s[3]][s[4]][s[5]][s[6]][s[7]][s[8]]=i;
             Push(s);
          }
       }
    }
}
void output()
{
    long t[17],ans[30],tot,tmp;
    for(long i=0;i<9;i++) t[i]=0;
    tot=0;
    while(f[t[0]][t[1]][t[2]][t[3]][t[4]][t[5]][t[6]][t[7]][t[8]]!=-1)
    {
       tot++;
       tmp=f[t[0]][t[1]][t[2]][t[3]][t[4]][t[5]][t[6]][t[7]][t[8]];
       ans[tot]=tmp+1;
       for(long i=0;i<9;i++)
         if(change[tmp][i])
           t[i]=(t[i]+3)%4;
    }
    bool first=true;
    for(long i=tot;i>=1;i--)
    {
       if(first) first=false;
       else printf(" ");
       printf("%ld",ans[i]);
    }
    printf("\n");
}
int main()
{
    init();
    bfs();
    output();
return 0;
}

lee1r 2010-04-07 21:47 发表评论

UVa 208 Firetruck

lee1r — Sat, 20 Mar 2010 04:03:00 GMT

题目大意：给出一个无边权无向图，起始结点为1，终止结点为t，求出从1到t的全部路径。
先求出和t属于同一个连通分量的结点，然后DFS。
以下是我的代码：

#include<stdio.h>
#include<string.h>
#define maxn 1007
#define max(a,b) (a>b?a:b)
long cas,aim,n,m,ans,a[maxn],way[maxn];
bool g[maxn][maxn],used[maxn];
void swap(long &a,long &b)
{
    long t=a;a=b;b=t;
}
void travel(long node)
{
    used[node]=true;
    m++;
    a[m]=node;
    for(long i=1;i<=n;i++)
      if(!used[i]&&g[node][i])
        travel(i);
}
void Qsort(long *a,long begin,long end)
{
    long i=begin,j=end,mid=a[(begin+end)/2];
    do{
         while(a[i]<mid) i++;
         while(a[j]>mid) j--;
         if(i<=j)
         {
            swap(a[i],a[j]);
            i++;j--;
         }
    }while(i<=j);
    if(begin<j) Qsort(a,begin,j);
    if(i<end)   Qsort(a,i,end);
}
void dfs(long node,long dep)
{
    if(node==aim)
    {
       bool first=true;
       for(long i=1;i<=dep;i++)
       {
          if(first) first=false;
          else printf(" ");
          printf("%ld",way[i]);
       }
       putchar('\n');
       ans++;
       return;
    }
    for(long i=1;i<=m;i++)
      if(!used[a[i]]&&g[node][a[i]])
      {
         way[dep+1]=a[i];
         used[a[i]]=true;
         dfs(a[i],dep+1);
         used[a[i]]=false;
      }
}
int main()
{
    long s,t;
    cas=0;
    while(scanf("%ld",&aim)==1)
    {
       cas++;
       n=ans=0;
       memset(g,false,sizeof(g));
       while(scanf("%ld%ld",&s,&t)==2)
       {
          if(s==0&&t==0) break;
          n=max(n,s);n=max(n,t);
          g[s][t]=g[t][s]=true;
       }
       memset(used,false,sizeof(used));
       m=0;
       travel(aim);
       Qsort(a,1,m);
       memset(way,0,sizeof(way));
       memset(used,false,sizeof(used));
       used[1]=true;
       way[1]=1;
       printf("CASE %ld:\n",cas);
       dfs(1,1);
       printf("There are %ld routes from the firestation to streetcorner %ld.\n",ans,aim);
    }
return 0;
}

lee1r 2010-03-20 12:03 发表评论