24/7/2012 ICPC培训 第九天

又来吐吐糟啦。。。

今天还算顺利，总共刷出了5题。

前两题是并查集，就是昨天没搞定的。后面三题都是最短路径了。

 

简单说说这几题吧。

关于并查集的两题其实都是卡在了路径压缩的处理上。很坑爹的说。其实也就是在路径压缩时

利用数组记录下一些信息。在需要的时候在合理的计算出来。

其实，挺难理解的。。。

我也只是一知半解了。。。

 

第一题（HDU1829）

主要是判断两个人是否在同一连通分量的同一层次上。那个计算层次的两个公式，较难理解。

代码：

#include<iostream>#include<cstring> using namespace std;const int maxn=2001;int father[maxn],layer[maxn],numBug,actionBug;bool flag;void init(){    memset(layer,0,sizeof(layer));         for(int i=1;i<=numBug;i++)    {        father[i]=i;    }}int find(int x){    if(x!=father[x])    {        int tmp=father[x];         father[x]=find(father[x]);         layer[x]=(layer[x]+layer[tmp])%2;    }        return father[x]; }void combine(int x,int y){    int fx=find(x);    int fy=find(y);        if(fx==fy)    {        if(layer[x]==layer[y])  flag=true;          }    else    {        father[fx]=fy;        layer[fx]=(layer[x]+layer[y]+1)%2;     }}void deal(){    flag=false;        int x,y;    for(int i=1;i<=actionBug;i++)    {        scanf("%d%d",&x,&y);        if(!flag) combine(x,y);    }}        int main(){    int cas,t=1;        scanf("%d",&cas);    while(cas--)    {        scanf("%d%d",&numBug,&actionBug);                init();                deal();                printf("Scenario #%d:\n",t++);        if(flag) printf("Suspicious bugs found!\n\n");        else printf("No suspicious bugs found!\n\n");    }        system("pause");     return 0;}

 

第二题（HDU3635）

这题的难点在于求ball的转移次数。

暴力解法是每次Q A时都暴力找到和A同祖先的ball，然后把转移次数加一，O(n^2)的算法，TLE了。

另外的解法就是在路径压缩时利用路径压缩的特点来解，难点也就是这个了。

如果你是我的祖先，你的转移次数加上我的转移次数才是我的转移次数，否则你早已不是我的祖先

了。只有直接转移祖先路径才不会被压缩，否则压缩。

代码：

#include<iostream>#include<cstring>using namespace std;const int maxn=10001;int father[maxn];int transport[maxn];int numBall[maxn]; int N,K,t=1;int find(int x){    if(x!=father[x])    {        int temp=father[x];        father[x]=find(father[x]);        transport[x]+=transport[temp];    }        return father[x];        //if(x==father[x])    //{    //    return x;    //}        //return father[x]=find(father[x]);}void combine(int x,int y){    int fx=find(x);    int fy=find(y);        if(fx==fy) return;        father[fx]=fy;    transport[fx]++;    numBall[fy]+=numBall[fx];}void init(){    memset(transport,0,sizeof(transport));        for(int i=1;i<=N;i++)    {        father[i]=i;        numBall[i]=1;    }} void inputAndQuary(){    int i,a,b;    char ch;        printf("Case %d:\n",t++);         for(i=0;i<K;i++)    {        getchar();            scanf("%c",&ch);                if(ch=='T')        {            scanf("%d%d",&a,&b);                        combine(a,b);        }        else        {            scanf("%d",&a);                        int tmp=find(a);            printf("%d %d %d\n",tmp,numBall[tmp],transport[a]);        }    }}                                  int main(){    int cas;        scanf("%d",&cas);    while(cas--)    {        scanf("%d%d",&N,&K);                init();                   inputAndQuary();    }        system("pause");    return 0;}

 

后面三题全是最短路径问题。（HDU2066、2544、1874）

还都是用SPFA算法来解的。代码都非常相似，我就放在一次讲啦。

不过，问题时我用模板过了2066、2544，却卡在了1874，数次WA呀。。。

这就让我纠结了。

神马情况呀。前面都过了。

找这种错误是最难的，最纠结的。因为我认为木有问题呀。从逻辑到数据都是经过检验的呀。

花了很长时间，一直忍着没找度娘。功夫不负有心人呀。最后发现模板中有一个函数的代码

有问题，就这样也过了2066、2544。。。这数据没的说了。。。

代码：

#include<iostream>#include<vector>#include<queue>#include<cstring>using namespace std;const int maxn=1001;const int INF=0x7fffffff;struct node{    int to;    int cost;};vector<node> myV[maxn];  //利用临界表存储图 int numRoad,numFrom,numTo; //路径数、可以开始的地点数、想去的地点数 int minPath[maxn],canFrom[maxn],wantTo[maxn]; //最短路、可以出发的城市、想去的城市 bool inQ[maxn];  //是否入队 /*  有bug的函数 bool judgeExistAndSmall(int a,int b,int time){    bool flag=false;        for(int i=0;i<myV[a].size();i++)    {         if(myV[a][i].to==b)        {            if(myV[a][i].cost>time)            {                myV[a][i].cost=time;                flag=true;                break;            }        }    }        if(flag)    {        for(int j=0;j<myV[b].size();j++)        {            if(myV[b][j].to==a)            {                myV[b][j].cost=time;                break;            }        }    }        if(flag) return true;    return false;}*/bool judgeExistAndSmall(int a,int b,int time){    bool flagExist=false,flagSmall=false;  //路径是否存在、存在的路径是否更小         for(int i=0;i<myV[a].size();i++)    {         if(myV[a][i].to==b)        {            flagExist=true;             if(myV[a][i].cost>time)            {                myV[a][i].cost=time;                flagSmall=true;                break;            }        }    }        if(flagSmall)  //双向图，只有存在的路径更小才需要处理     {        for(int j=0;j<myV[b].size();j++)        {            if(myV[b][j].to==a)            {                myV[b][j].cost=time;                break;            }        }    }        if(flagExist) return true;   //只要存在路径就返回true     return false;}     void inputItial(){    int i,a,b,time;        for(i=0;i<maxn;i++)  //清空再使用     {        myV[i].clear();    }        for(i=0;i<numRoad;i++)    {        scanf("%d%d%d",&a,&b,&time);                if(!judgeExistAndSmall(a,b,time))  //无向图、判断是否存在了路、存在的路是大是小         {            node tmp;            tmp.cost=time;                        tmp.to=b;            myV[a].push_back(tmp);                        tmp.to=a;            myV[b].push_back(tmp);                  }    }        for(i=0;i<numFrom;i++)    {        scanf("%d",&canFrom[i]);    }        for(i=0;i<numTo;i++)    {        scanf("%d",&wantTo[i]);    }}        void SPFA()   //最短路径快速算法 Shortest Path Faster Algorithm{    int minCost=INF;        for(int i=0;i<numFrom;i++)    {        memset(inQ,false,sizeof(inQ));        for(int j=0;j<maxn;j++)        {            minPath[j]=INF;        }        minPath[canFrom[i]]=0;                queue<int> myQ;        myQ.push(canFrom[i]);        inQ[canFrom[i]]=true;                int now,to,cost;        while(!myQ.empty())        {            now=myQ.front();            myQ.pop();                        for(int k=0;k<myV[now].size();k++)            {                to=myV[now][k].to;                cost=myV[now][k].cost+minPath[now];                if(minPath[to]>cost)                {                    minPath[to]=cost;                                        if(!inQ[to])                    {                        inQ[to]=true;                        myQ.push(to);                    }                }            }                        inQ[now]=false;        }                for(int kk=0;kk<numTo;kk++)        {            minCost=min(minCost,minPath[wantTo[kk]]);        }    }        printf("%d\n",minCost);}     int main(){    while(scanf("%d%d%d",&numRoad,&numFrom,&numTo)==3)    {        inputItial();                SPFA();    }     system("pause");       return 0;}

 

最后，到这么晚才写实在不好意思。

看电视的说。。。