C++博客-Uriel's Corner-文章分类-POJ

POJ 3450 Corporate Identity---KMP or strstr or Suffix Array

Uriel — Tue, 10 May 2011 07:24:00 GMT

摘要: 最近一直在准备省赛,基本都是复习算法和组队赛,基本没切什么题. 昨天想省赛前POJ刷到1200先,就跟光光对拍了下POJ过的题. 发现了这道字符串,顺便复习下KMP和后缀数组啥的. 这题跟POJ 3080基本完全一样,当年写过个解题报告 http://www.cppblog.com/Urie... 阅读全文

Uriel 2011-05-10 15:24 发表评论

POJ 3613 Cow Relays---Floyd+矩阵相乘

Uriel — Mon, 02 May 2011 09:39:00 GMT

这是在DY大牛第六期专题里列出的题目之一.

话说搞ACM也2年多了, 竟然用矩阵乘法+Floyd的题目都还没有做过... 上周的组队赛遇到一道题才听说... (没有拜读matrix67的神文... 检讨下... )

/**//*
POJ 3613 Cow Relays

-------Classify: Floyd+矩阵乘法
----Description: 求从S到T恰好经过K条边(可重复走)的最短路
-----------------K<=1e6, m<100(边数)
---Sample Input:

2 6 6 4----------//K m s t
11 4 6-----------//length x<->y
4 4 8
8 4 9
6 6 8
2 6 9
3 8 9

--Sample Output:

10

-----Time Limit: 1000Ms
---------Source: USACO 2007 November Gold
-------Solution: 

    By Solution Report

    01邻接矩阵A的K次方C=A^K,C[i][j]表示i点到j点正好经过K条边的路径数

          对应于这道题,对邻接图进行K次floyd之后,C[i][j]就是点i到j正好经过K条边的最短路


---------Status: AC C++ 532Ms
-----------Date: 2011.05.02
------Reference: http://hi.baidu.com/aconly/blog/item/23fb73acc1d874004b36d634.html
-----------Code: 
*/

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#define N 1010
#define INF 0x3f3f3f3f
typedef int M[N][N];

int k, m, s, t, v;
int ans[N][N], adj[N][N], vis[N], vtx[N], dis[N][N], tp[N][N];

void floyd(M c, M a, M b) {
    int i, j, k;
    for (k = 0; k < v; ++k)
        for (i = 0; i < v; ++i)
            for (j = 0; j < v; ++j)
                if (c[vtx[i]][vtx[j]] > a[vtx[i]][vtx[k]] + b[vtx[k]][vtx[j]])
                    c[vtx[i]][vtx[j]] = a[vtx[i]][vtx[k]] + b[vtx[k]][vtx[j]];
}

void copy(M a, M b) {
    int i, j;
    for (i = 0; i < v; ++i)
        for (j = 0; j < v; ++j) {
            a[vtx[i]][vtx[j]] = b[vtx[i]][vtx[j]];
            b[vtx[i]][vtx[j]] = INF;
        }
}

void BS(int k) {
    while (k) {
        if (k & 1) {
            floyd(dis, ans, adj);
            copy(ans, dis);
        }
        floyd(tp, adj, adj);
        copy(adj, tp);
        k >>= 1;
    }
}

int main() {
    int i, j, x, y, w;
    scanf("%d %d %d %d", &k, &m, &s, &t);
    for (i = 0; i <= 1001; ++i) {
        for (j = 0; j <= 1001; ++j) {
            adj[i][j] = INF;
            tp[i][j] = INF;
            dis[i][j] = INF;
            ans[i][j] = INF;
        }
        ans[i][i] = 0;
    }
    v = 0;
    memset(vis, 0, sizeof(vis));
    for (i = 0; i < m; ++i) {
        scanf("%d %d %d", &w, &x, &y);
        if (!vis[x]) {
            vis[x] = 1;
            vtx[v++] = x;
        }
        if (!vis[y]) {
            vis[y] = 1;
            vtx[v++] = y;
        }
        if (adj[x][y] > w)
            adj[x][y] = adj[y][x] = w;
    }
    BS(k);
    printf("%d\n", ans[s][t]);
    return 0;
}

Uriel 2011-05-02 17:39 发表评论

POJ 1625 Censored!---AC 自动机+DP+高精度

Uriel — Sat, 30 Oct 2010 15:40:00 GMT

摘要: 第二道AC自动机，对DP状态转移有一点感觉了。。。这题基本和2778一样，不过要高精度，所以就想用JAVA试试，正好JAVA很久不碰了，跟没学过一样，Regional之前练习练习，结果做得我痛不欲生。。00：00-03：10，0... 阅读全文

Uriel 2010-10-30 23:40 发表评论

POJ 3250 Bad Hair Day---栈的应用

Uriel — Thu, 21 Oct 2010 07:08:00 GMT

        跟别人Compare得到的水题，虽然是想拿来轻松一下，不过还是有些收获的~
        看完题本来想5min切掉的，结果搞了半小时。。= =。。思路有点混乱。。。

        stk[i]表示牛i右侧不比i低的第一头牛，累加更新。。(单步之后才搞清楚。。)

        加了读入优化之后32Ms。。还行。。

//Problem: 3250  User: Uriel
//Memory: 1004K  Time: 32MS
//Language: G++  Result: Accepted

#include<stdio.h>
#include<stdlib.h>

int stk[80010],a[800010];

int in(){
    char ch;
    int a=0;
    while((ch=getchar())==' ' || ch=='\n');
    a*=10;
    a+=ch-'0';
    while((ch=getchar())!=' ' && ch!='\n'){
        a*=10;
        a+=ch-'0';
    }
    return a;
}

int main(){
    int i,j,n,k=1;
    __int64 res=0,t;
    n=in();
    for(i=0;i<n;++i)a[i]=in();
    a[n]=0x3fffffff;
    for(i=n-1;i>=0;--i){
        t=0;
        stk[i]=i+1;
        while(a[i]>a[stk[i]] && stk[i]<n){
            t+=stk[stk[i]]-stk[i];
            stk[i]=stk[stk[i]];
        }
        res+=t;
    }
    printf("%I64d\n",res);
    return 0;
}

Uriel 2010-10-21 15:08 发表评论

POJ 3759 Simple Distributed computing system---最大流

Uriel — Wed, 13 Oct 2010 15:23:00 GMT

摘要: 网络流一拖再拖，发现再不临时抱佛脚一下Regional我们队要悲剧了。。= =。。于是刚看了几天最大流。。看某分类说这题是最大流，但是这题AC甚少，看到Discuss说很裸于是看了题，但是题目... 阅读全文

Uriel 2010-10-13 23:23 发表评论

POJ 1851 Map---计算几何

Uriel — Sun, 26 Sep 2010 18:23:00 GMT

摘要: ECUST ACHunter威武，今天哈尔滨现场赛一度rank6(算上旅游队)，后2h可惜了，福州继续加油~~ 继续刷计算几何。。不会... 阅读全文

Uriel 2010-09-27 02:23 发表评论

POJ 1380 Equipment Box && POJ 2172 Bricks---计算几何

Uriel — Mon, 20 Sep 2010 14:06:00 GMT

        POJ 1380是去年十月做的，今天发现POJ 2172与POJ 1380一样的。。拿原来代码直接贴了。。

        1380是问一个矩形能否塞进另一个矩形，2172其实一回事，虽然变成三维，但是只要判三次就行了。。

如图，假设小矩形正好卡进去的情况，求出边长Len需要多少，Len>a则塞不进

代码中：∠1=∠sita      ∠1+∠2=∠p

以下是2172的代码，1380判定过程跟这个一样，而且只要判一次就行

//Problem: 2172  User: Uriel
//Memory: 208K  Time: 32MS
//Language: C++  Result: Accepted
//2010.09.20

#include<math.h>
#include<stdio.h>
#include<stdlib.h>

struct point {
    double x,y,h;
}a,b;

void Swap(double & a,double & b){
    double t;
    t=a;
    a=b;
    b=t;
}

bool check(double c,double d,double a,double b){
    if(a*b<c*d)return false;
    if(a<b)Swap(a,b);
    if(c<d)Swap(c,d);
    if(b<d)return false;
    if(a>=c && b>=d)return true;
    double dis=sqrt(c*c+d*d);
    double p=asin(b/dis);
    double q=asin(d/dis);
    double sita=p-q;
    double Len=c*cos(sita)+d*sin(sita);
    if(Len<=a)return true;
    return false;
}

int main(){
    scanf("%lf %lf %lf %lf %lf",&a.x,&a.y,&a.h,&b.x,&b.y);
    if(check(a.x,a.y,b.x,b.y)||check(a.x,a.h,b.x,b.y)||check(a.y,a.h,b.x,b.y))puts("YES");
    else
        puts("NO");
    return 0;
}

Uriel 2010-09-20 22:06 发表评论

POJ 2284 That Nice Euler Circuit---计算几何

Uriel — Sat, 18 Sep 2010 17:45:00 GMT

摘要: 这几天做计算几何一直不顺，各种WA。。这题做了一天多。。WA到哭了。。这题思路很简单，用平面图的欧拉定理，V+F-E=2；其中V是顶点数，F是分割出的面数，E是边数。 &... 阅读全文

Uriel 2010-09-19 01:45 发表评论

POJ 1931 Biometrics---计算几何

Uriel — Fri, 17 Sep 2010 09:59:00 GMT

        5门专业必修+4门大类选修真纠结。。微机、计算机图像处理还行，数据库也还行，计算机组成原理算得头大，最纠结的还是编译原理，K老师的课漏听30min肯定就听不懂了，作业又一堆。。每天白天都用来上课+写作业了。。就晚上选修课结束回1教切会题。。

        看起来很水的一道计算几何卡了两天才过。。悲剧啊。。

        一开始没注意相似的点是对应的，不会说第一个多边形的第1个点对应第二个多边形的第3个点之类的，想了几个WS算法，WA。。今天注意到了这点，结果改来改去还是不对。。最后看了某解题报告，把判转向那里改了一下，总算过了。。

        代码比较挫。。很多无用的东西懒得擦掉了。。

#include<math.h>
#include<stdio.h>
#include<stdlib.h>
#include<algorithm>
using namespace std;
#define eps 1e-6
const double INF=1e20;

struct point{
    double x,y;
    point operator-(point &b){
        point c;
        c.x = x - b.x;
        c.y = y - b.y;
        return c;
    }
}p1[11],p2[11];

int n;
int dir1[11],dir2[11];
double len1[11],len2[11];

double dis(point a,point b){
    return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}

double Cos(point p0,point p1,point p2){
    return (dis(p0,p1)+dis(p1,p2)-dis(p0,p2))/(2*sqrt(dis(p0,p1))*sqrt(dis(p1,p2)));
}

int cross_product(point p1,point p2,point p0){
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}

double dot_product(point p1,point p2,point p0){
    return (p1.x-p0.x)*(p2.x-p0.x)+(p1.y-p0.y)*(p2.y-p0.y);
}

double cross(point a, point b){
    return a.x * b.y - a.y * b.x;
}

int get_dir(point p1,point p2,point p3){
    double t1=cross(p2-p1,p3-p2);
    if(fabs(t1)<eps)return 1;
    if(t1<0)return 2;
    else return 3;
}

int main(){
    int f1,f2,i,j,cnt;
    bool ok;
    double minn1,len,minn2,cons;
    while(scanf("%d",&n),n){
        for(i=0;i<n;i++)scanf("%lf %lf",&p1[i].x,&p1[i].y);
        for(i=0;i<n;i++)scanf("%lf %lf",&p2[i].x,&p2[i].y);
        cons=sqrt(dis(p1[0],p1[1]))/sqrt(dis(p2[0],p2[1]));
        for(i=0;i<n;i++){
            len1[i]=sqrt(dis(p1[i],p1[(i+1)%n]))/cons;
            len2[i]=sqrt(dis(p2[i],p2[(i+1)%n]));
        }
        ok=true;
        for(i=0;i<n;i++){
            if(fabs(len1[i]-len2[i])>eps){
                ok=false;
                break;
            }
        }
        if(ok){
            for(i=0;i<n;i++){
                if(fabs(fabs(Cos(p1[i],p1[(i+1)%n],p1[(i+2)%n]))-fabs(Cos(p2[i],p2[(i+1)%n],p2[(i+2)%n])))>eps){
                    ok=false;
                    break;
                }
            }
        }
        if(ok){
            for(i=0;i<n;i++){
                if(get_dir(p1[i],p1[(i+1)%n],p1[(i+2)%n])!=get_dir(p2[i],p2[(i+1)%n],p2[(i+2)%n])){
                    ok=false;
                    break;
                }
            }
        }
        if(ok)puts("similar");
        else
            puts("dissimilar");
    }
    return 0;
}

Uriel 2010-09-17 17:59 发表评论

POJ 3129 How I Wonder What You Are!----计算几何

Uriel — Tue, 07 Sep 2010 15:45:00 GMT

       找了一晚上计算几何题，看一道一道不会的。。终于找到一道大水。。。

       题意是有天空有N颗星星，给出它们的三维坐标。。在(0,0,0)处放着M个望远镜，给出它们的仰角(tx,ty,tz)以及望远镜圆锥形的顶角phi，求总共能看见多少星星。。

      方法就是暴力N*M，一个个的判，貌似是第一还是第二次用浙大模板的三维几何~~

//Problem: 3129  User: Uriel
//Memory: 228K  Time: 47MS
//Language: C++  Result: Accepted

#include<math.h>
#include<stdio.h>
#include<stdlib.h>

struct point{
    bool flag;
    double x,y,z;
}p[550];

struct tele{
    double x,y,z,tx,ty,tz,phi;
}tel[55];

int n,m;

double vlen(point p){
    return sqrt(p.x*p.x+p.y*p.y+p.z*p.z);
}

double dmult(point u,point v){
    return u.x*v.x+u.y*v.y+u.z*v.z;
}

point subt(point u,point v){
    point ret;
    ret.x=u.x-v.x;
    ret.y=u.y-v.y;
    ret.z=u.z-v.z;
    return ret;
}

double angle_cos(point u1,point u2,point v1,point v2){
    return dmult(subt(u1,u2),subt(v1,v2))/vlen(subt(u1,u2))/vlen(subt(v1,v2));
}

int main(){
    int i,j,cnt;
    while(scanf("%d",&n),n){
        for(i=0;i<n;i++){
            p[i].flag=false;
            scanf("%lf %lf %lf",&p[i].x,&p[i].y,&p[i].z);
        }
        scanf("%d",&m);
        for(i=0;i<m;i++){
            scanf("%lf %lf %lf %lf",&tel[i].tx,&tel[i].ty,&tel[i].tz,&tel[i].phi);
            tel[i].x=0.0;
            tel[i].y=0.0;
            tel[i].y=0.0;
        }
        cnt=0;
        for(i=0;i<m;i++){
            point tmp,ori;
            ori.x=0.0;
            ori.y=0.0;
            ori.z=0.0;
            tmp.x=tel[i].tx;
            tmp.y=tel[i].ty;
            tmp.z=tel[i].tz;
            for(j=0;j<n;j++){
                if(p[j].flag)continue;
                if(acos(angle_cos(ori,tmp,ori,p[j]))<tel[i].phi){
                    cnt++;
                    p[j].flag=true;
                }
            }
        }
        printf("%d\n",cnt);
    }
    return 0;
}

Uriel 2010-09-07 23:45 发表评论

POJ 1578 Instruens Fabulam---模拟

Uriel — Fri, 03 Sep 2010 14:37:00 GMT

摘要: 很久不切模拟，今天比赛的时候手有点生了。。。这题断断续续做了两天，被各种事情打乱。。纯模拟，方法没什么可说的，这题跟3699有点类似，比3699简单很多... 阅读全文

Uriel 2010-09-03 22:37 发表评论

POJ 1687 Buggy Sat---计算几何

Uriel — Sat, 28 Aug 2010 15:22:00 GMT

计算几何大水题，题目求哪个多边形是最外层的，只要算每个多边形的面积，找出最大的一个就行，因为题目给的坐标标号已经是排过序的，所以求面积只要叉积即可~

//Problem: 1687  User: Uriel
//Memory: 388K  Time: 0MS
//Language: G++  Result: Accepted
//Computational Geometry -> Cross Product
//2010.08.29

#include<math.h>
#include<stdio.h>
#include<stdlib.h>

struct point{
    int x,y;
}p[55];

int cross_product(point a,point b)
{
    return a.x*b.y-a.y*b.x;
}

int main(){
    int i,j,cse,n,x,y,tmp,id,pre,sum,maxs,ans;
    scanf("%d",&cse);
    while(cse--){
        scanf("%d",&n);
        for(i=0;i<n;i++){
            scanf("%d %d",&p[i].x,&p[i].y);
        }
        scanf("%d",&x);
        maxs=-1;ans=0;
        for(i=0;i<x;i++){
            scanf("%d",&y);
            sum=0;
            for(j=0;j<y;j++){
                scanf("%d",&id);
                id--;
                if(j==0){
                    pre=id;
                    tmp=id;
                }
                else{
                    sum+=cross_product(p[pre],p[id]);
                    pre=id;
                }
            }
            sum+=cross_product(p[id],p[tmp]);
            if(sum>maxs){
                maxs=sum;
                ans=i+1;
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}

Uriel 2010-08-28 23:22 发表评论

POJ 3481 Double Queue---SBT/STL

Uriel — Wed, 18 Aug 2010 14:52:00 GMT

        这题貌似最优解法是SBT，这个只看过，完全没有敲过，想到STL set，于是水过去了。。。
        对于STL，现在基本也只会priority queue, set, heap, map还基本没怎么写过，只能应付应付比赛的时候最基本的用法。。。
        这题用set，C++比G++快一倍，282Ms，有点慢，不过代码量很少~~

//Problem: 3481  User: Uriel
//Memory: 840K  Time: 282MS
//Language: C++  Result: Accepted
//Version: STL => set
//2010.08.18

#include<set>
#include<stdio.h>
#include<stdlib.h>
using namespace std;

struct node{
    int id,p;
}q;

bool operator<(const node &a,const node &b){
    return a.p>b.p;
}

int main(){
    int x;
    set<node> st;
    while(scanf("%d",&x),x){
        if(x==1){
            scanf("%d %d",&q.id,&q.p);
            st.insert(q);
        }
        else if(x==2){
            if(st.empty()){
                puts("0");
                continue;
            }
            printf("%d\n",st.begin()->id);
            st.erase(st.begin());
        }
        else if(x==3){
            if(st.empty()){
                puts("0");
                continue;
            }
            set<node>::iterator it;
            it=st.end();
            it--;
            printf("%d\n",it->id);
            st.erase(st.find(*it));
        }
    }
    return 0;
}

Uriel 2010-08-18 22:52 发表评论

POJ 1208 The Blocks Problem---模拟

Uriel — Thu, 12 Aug 2010 04:16:00 GMT

摘要: 集训结束了，最后一天想切道水题休息休息，结果被虐一晚上。。。看了Discuss才知道题目都理解错了。。晚上回宿舍改改，早上终于AC了。。 //Problem: 1208 User: U... 阅读全文

Uriel 2010-08-12 12:16 发表评论

POJ 1791 Paper Cutting---模拟

Uriel — Wed, 11 Aug 2010 06:36:00 GMT

        又是一道去年没切掉的模拟，去年大略看了一下，题目不是很懂就没做。
        今天又翻出来这道题。。果真大水。。= =

        题意是印刷切割名片，印的时候要A*B张一起印刷，每张名片大小为C*D，印刷纸张大小E*F，问印完之后至少切一刀能把名片都分开(不能把名片叠在一起一刀切)，明白意思之后就显然大水一道了~~

//Problem: 1791  User: Uriel
//Memory: 368K  Time: 0MS
//Language: G++  Result: Accepted
//Simulation
//2010.08.11
#include<stdio.h>
#include<stdlib.h>
#include<algorithm>
using namespace std;
#define INF 100000000
int A,B,C,D,E,F,minx,cut;

int main(){
    int hh,ww;
    while(scanf("%d %d %d %d %d %d",&C,&D,&A,&B,&E,&F),A|B|C|D|E|F){
        minx=INF;
        if(A>B)swap(A,B);
        if(C>D)swap(C,D);
        if(E>F)swap(E,F);
        //-----------------------------------------------case 1
        cut=0;
        hh=A*C;ww=B*D;
        if(hh>E || ww>F)goto case2;
        if(E>hh)cut++;
        if(F>ww)cut++;
            cut+=C*D-1;
        if(cut<minx)minx=cut;
        //-----------------------------------------------case 2
case2:    cut=0;
        hh=A*D;ww=B*C;
        if(hh>ww)swap(hh,ww);
        if(hh>E || ww>F)goto flag;
        if(E>hh)cut++;
        if(F>ww)cut++;
        cut+=C*D-1;
        if(cut<minx)minx=cut;
flag:    if(minx<INF)printf("The minimum number of cuts is %d.\n",minx);
        else
            puts("The paper is too small.");
    }
    return 0;
}

Uriel 2010-08-11 14:36 发表评论

POJ 1169 Packing Rectangles---模拟

Uriel — Tue, 10 Aug 2010 03:57:00 GMT

摘要: 比较麻烦的模拟，倒不是方法难，主要情况多，过程很简单，就是题目给的那几种情况枚举(开始以为随便放置，所以一直以为是大自然题，所以放了一年没做。。)，题目的case 4 & case 5 可以合并~ 每种case都要考虑横放... 阅读全文

Uriel 2010-08-10 11:57 发表评论

POJ 2942 Knights of the Round Table---双连通分量 & 判二分图

Uriel — Sun, 08 Aug 2010 13:02:00 GMT

一切均参考了DY大牛的Blog，http://www.answeror.com/?s=2942
丑陋的代码就不贴了。。。

Uriel 2010-08-08 21:02 发表评论

POJ 3713 Transferring Sylla---判3-连通

Uriel — Sun, 08 Aug 2010 12:54:00 GMT

判一个无向图是不是至少3-连通
方法是枚举点，去掉一个点之后判余下的图中是否存在割点，如果都存在，则是至少3-连通，否则就不是

//Problem: 3713  User: Uriel
//Memory: 2164K  Time: 3922MS
//Language: C++  Result: Accepted
//enumration & cut vertex
//2010.08.07
#include<stdio.h>
#include<stdlib.h>
#include<algorithm>
using namespace std;
#define N 510

int n,m,cnt,root;
int ancestor[N];
int mark[N];
int deep[N];
int b[N][N];
int adj[N][N];
bool flag;

int DFS(int u,int fa,int t){
    mark[u]=1;
    deep[u]=ancestor[u]=t;
    int cnt=0;
    int i,v;
    for(i=b[u][0];i>=1;i--){
        v=b[u][i];
        if(mark[v]==2)continue;
        if(mark[v]==0){
            DFS(v,u,t+1);
            cnt++;
            ancestor[u]=min(ancestor[u],ancestor[v]);
            if(u==root && cnt>1 || u!=root && ancestor[v]>=deep[u])flag=1;
        }
        else if(mark[v]==1 && v!=fa){
            ancestor[u]=min(ancestor[u],deep[v]);
        }
    }
    return 0;
}

int main(){
    int x,y,i,j;
    while(scanf("%d %d",&n,&m),n|m){
        memset(adj,0,sizeof(adj));
        for(i=0;i<n;i++)b[i][0]=0;
        for(i=0;i<m;i++){
            scanf("%d %d",&x,&y);
            adj[x][y]=adj[y][x]=1;
        }
        for(i=0;i<n;i++){
            for(j=0;j<n;j++){
                if(adj[i][j])b[i][++b[i][0]]=j;
            }
        }
        flag=false;
        for(i=0;i<n;i++){
            memset(mark,0,sizeof(mark));
            memset(deep,0,sizeof(deep));
            memset(ancestor,0,sizeof(ancestor));
            mark[i]=2;
            if(i==0)root=1;
            else
                root=0;
            DFS(root,-1,1);
            for(j=0;j<n;j++){
                if(mark[j]==0){
                    flag=true;
                    break;
                }
            }
            if(flag)break;
        }
        if(flag)puts("NO");
        else
            puts("YES");
     }
     return 0;
}

Uriel 2010-08-08 20:54 发表评论

POJ 1144 Network---求无向图割点数

Uriel — Sun, 08 Aug 2010 11:46:00 GMT

赤裸裸的求割点数的题目，上模板~~
关于求割点的方法，黑书上讲得比较详细了

//Problem: 1144  User: Uriel
//Memory: 228K  Time: 16MS
//Language: C++  Result: Accepted
//Cut Vertex
//2010.08.06
#include<stdio.h>
#include<stdlib.h>
#include<algorithm>
using namespace std;
#define N 110

int n;
int ancestor[N];
int mark[N];
int deep[N];
int cnt[N];
int adj[N][N];

int DFS(int i,int father,int dth,int b[][N],int cut[]){
    int j,k,sons=0,count=0;
    mark[i]=1;
    deep[i]=ancestor[i]=dth;
    for(k=b[i][0];k>=1;k--){
        j=b[i][k];
        if(j!=father && mark[j]==1)ancestor[i]=min(ancestor[i],deep[j]);
        if(mark[j]==0){
            count+=DFS(j,i,dth+1,b,cut);
            sons++;
            ancestor[i]=min(ancestor[i],ancestor[j]);
            if((father==-1&&sons>1)||(father!=-1&&ancestor[j]>=deep[i]))cut[i]=1;
            //if(ancestor[j]>deep[i]) brige[i][j]=1;
        }
    }
    mark[i]=2;
    return count+cut[i];
}

int cutpoints(int a[][N],int n,int cut[]){
    int i,j,b[N][N],count=0;
    for(i=0;i<n;i++)cut[i]=b[i][0]=0;
    memset(mark,0,sizeof(mark));
    for(i=0;i<n;i++)for(j=0;j<n;j++)if(a[i][j])b[i][++b[i][0]]=j;
    for(i=0;i<n;i++)if(mark[i]==0)count+=DFS(0,-1,0,b,cut);
    return count;
}

int main(){
    int x,y;
    char c;
     while(1){
         scanf("%d",&n);
         if(!n)return 0;
         memset(adj,0,sizeof(adj));
         while(1){
             scanf("%d",&x);
             if(x==0)break;
             x--;
             while(1){
                 scanf("%d",&y);
                 y--;
                 adj[x][y]=1;adj[y][x]=1;
                 c=getchar();
                 if(c=='\n')break;
             }
         }
         printf("%d\n",cutpoints(adj,n,cnt));
     }
     return 0;
}

Uriel 2010-08-08 19:46 发表评论

POJ 1523 SPF---DFS

Uriel — Sun, 08 Aug 2010 11:42:00 GMT

求割点没什么问题，但是求去掉这个割点图被分为几个连通分支没什么思路，看了解题报告恍然大悟，DFS就行了啊~~
枚举点，去掉看暴力DFS看被分成几块，判了是否是割点的同时也知道了连通分支数

//Problem: 1523  User: Uriel
//Memory: 4136K  Time: 16MS
//Language: C++  Result: Accepted
//Version: DFS
//2010.08.07
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#define N 1010
int num;
int map[N][N],arr[N];
bool vis[N];

int Loc(int k){
    int i=0;
    for(i=0;i<num;i++){
        if(arr[i]==k)break;
    }
    if(i==num){
        arr[num++]=k;
        return num-1;
    }
    return i;
}

void DFS(int x){
    vis[x]=true;
    for(int j=0;j<num;j++){
        if(!vis[j] && map[x][j])DFS(j);
    }
}

int main(){
    int g=1,x,y,a,b,l,s;
    while(scanf("%d",&x),x){
        scanf("%d",&y);
        num=0;
        a=Loc(x);b=Loc(y);
        map[a][b]=map[b][a]=1;
        while(scanf("%d",&x),x){
            scanf("%d",&y);
            a=Loc(x);b=Loc(y);
            map[a][b]=map[b][a]=1;
        }
        printf("Network #%d\n",g++);
        l=num-1;
        for(int i=0;i<num;i++){
            memset(vis,false,sizeof(vis));
            s=0;
            vis[i]=true;
            for(int j=0;j<num;j++)
                if(!vis[j]){
                    DFS(j);
                    s++;
                }
            if(s>1){
                printf("  SPF node %d leaves %d subnets\n",arr[i],s);
            }
            else
                l--;
        }
        if(l==-1)printf("  No SPF nodes\n");
        printf("\n");
        memset(map,0,sizeof(map));
    }
    return 0;
}

Uriel 2010-08-08 19:42 发表评论

POJ 3715 Blue and Red---二分图匹配

Uriel — Sun, 08 Aug 2010 11:34:00 GMT

摘要: 想到二分图，但是不知道那种贪心思想对不对以及如何优化。。参考了解题报告感谢大牛的解题报告：http://hi.baidu.com/liuzhe/blog/item/793e0c24efa6602cd407428a.html基本上照着写的，然后加了读入输出优化，效果惊人。。0MS。。暂时Rank1了。。还是头一次刷到这么前。。 //Problem: 3715 Use... 阅读全文

Uriel 2010-08-08 19:34 发表评论

POJ 2724 Purifying Machine---二分图匹配

Uriel — Mon, 02 Aug 2010 14:42:00 GMT

建图啊建图。。
思路完全膜拜AC大牛的代码http://hi.baidu.com/aekdycoin/blog/item/2dac891ea09ef9f2e1fe0b67.html....

代码如下：

//Problem: 2724  User: Uriel
//Memory: 384K  Time: 110MS
//Language: G++  Result: Accepted
//floyd+Bipartite Graph
//2010.08.02

#include<stdio.h>
#include<stdlib.h>
#include<string.h>

int n,m;
char a[1050][11];
int d[1050];
int nx,num[1050];
int f[11]={1,2,4,8,16,32,64,128,256,512,1024};
int link[1050],in[1050];

bool map(int a , int b){
    int t1=num[a],t2=num[b];
    int m=t1^t2;
    if(m && ((m&(m-1))==0))return true;
    return false;
}

bool fun(int k){
    int i;
    for(i=0;i<nx;i++){
        if(i==k)continue;
        if(map(i,k) && !in[i]){
            in[i]=1;
            if(link[i]==-1 || fun(link[i])){
                link[i]=k;
                return true;
            }
        }
    }
    return false;
}

int main(){
    while(scanf("%d%d",&n,&m)){
        if(n==0 && m==0)break;
        nx=0;
        int i,j,t;
        memset(d,0,sizeof(d));
        for(i=0;i<m;i++){
            scanf("%s",(a[i]));
            int flag=-1;
            t=0;
            for(j=0;j<n;j++){
                if(a[i][j]=='*'){
                    flag=j;
                    continue;
                }
                t+=f[j]*(a[i][j]-'0');
            }
            if(!d[t]){
                d[t]=1;
                num[nx++]=t;
            }
            if(flag!=-1 && !d[t+f[flag]]){
                d[t+f[flag]]=1;
                num[nx++]=t+f[flag];
            }
        }
        int res=0;
        memset(link,-1,sizeof(link));
        for(i=0;i<nx;i++){
            memset(in,0,sizeof(in));
            if(fun(i))res++;
        }
        printf("%d\n",nx-res/2);
    }
    return 0;
}

Uriel 2010-08-02 22:42 发表评论

POJ 1984 Navigation Nightmare---并查集

Uriel — Mon, 02 Aug 2010 11:12:00 GMT

题意：给出点数和边数
            已知每一条边的起点终点编号和方向
            几个Query，每个Query求两个编号的点在第 I 步之后的曼哈顿距离

思路：先将Query按操作步数的大小排序，从小到大
            遍历每个操作，执行并集函数，然后如果这步有Query则对Query的两个点执行查集函数，如果在一个集合当中，他们的曼哈顿距离即为他们的父亲结点的曼哈顿距离，否则输出-1

            并集函数写法：如果第x条边(第x 步操作产生的)两个端点不在同一集合，则合并，注意合并之后作为儿子的结点到父亲结点的曼哈顿距离要有变化，要分边的方向不同考虑，详见代码

            查集函数注意x，y 相对距离的累加

这题G++始终过不了，开始以为是数组越界或是浮点误差，但是都排除了，还怀疑是abs()的问题，手写了ABS也不行。。有路过的大牛请不吝赐教~

另外发现改手写ABS()之后从172MS降到125MS，很神奇~~

代码如下：

//Problem: 1984  User: Uriel
//Memory: 1352K  Time: 125MS
//Language: C++  Result: Accepted
//Union-Find Set
//2010.08.02

#include<stdio.h>
#include<stdlib.h>
#include<algorithm>
using namespace std;
#define MAXN 40010
#define ABS(x) ((x)>=0?(x):(-x))

struct Edge{
    int a,b,w;
    char dir;
}e[MAXN];

struct query{
    int idx,a,b;
}q[10010];

struct node{
    int f,a,b;
}v[MAXN];

int n,m,k;
char d[4]={'N','E','S','W'};
int step[4][2]={{0,1},{1,0},{0,-1},{-1,0}};

int ex_dir(char c){
    for(int i=0;i<4;i++)
        if(d[i]==c)return i;
}

bool cmp(query a,query b){
    return a.idx<b.idx;
}

int findset(int idx,int &x,int &y){
    int t=idx;
    x=y=0;
    while(v[t].f!=-1){
        x+=v[t].a;
        y+=v[t].b;
        t=v[t].f;
    }
    return t;
}

void unionset(int x){
    int x1,x2,y1,y2;
    int v1=e[x].a,v2=e[x].b;
    int f1=findset(v1,x1,y1),f2=findset(v2,x2,y2);
    if(f1==f2)return;
    int dd=ex_dir(e[x].dir);
    v[f2].f=f1;
    v[f2].a=step[dd][0]*e[x].w+x1-x2;
    v[f2].b=step[dd][1]*e[x].w+y1-y2;
}

int main(){
    int i,j,f1,f2,x1,x2,y1,y2;
    scanf("%d %d",&n,&m);
    for(i=0;i<m;i++){
        scanf("%d %d %d %c",&e[i].a,&e[i].b,&e[i].w,&e[i].dir);
    }
    scanf("%d",&k);
    for(i=0;i<k;i++){
        scanf("%d %d %d",&q[i].a,&q[i].b,&q[i].idx);
    }
    sort(q,q+k,cmp);
    for(i=0;i<=n;i++)v[i].f=-1;
    for(i=0,j=0;i<m;i++){
        unionset(i);
        while(j<k && q[j].idx<=i+1){
            f1=findset(q[j].a,x1,y1);
            f2=findset(q[j].b,x2,y2);
            if(f1==f2)printf("%d\n",ABS(x1-x2)+ABS(y1-y2));
            else
                printf("-1\n");
            j++;
        }
    }
    return 0;
}

Uriel 2010-08-02 19:12 发表评论

POJ 3216 Repairing Company---多源最短路+二分图匹配

Uriel — Mon, 02 Aug 2010 10:53:00 GMT

继续纠结各种建图...

题意：Q个blocks，M个tasks
            Q*Q的矩阵，给出每个边的情况
            对于每个task，已知pi(位于哪个block)，ti(最迟开始时间)，di(每个人完成它所需时间)

思路：建图还是很纠结，看了大牛们的思路...
            先用floyd求两点之间的最短路，然后建图。
            建图方法是如果task i 的开始时间+task i 需要的时间+从task i 所在地到task j 所在地所需时间<=task j的开始时间，则i，j连一条边
            然后就是求最小路径覆盖

代码如下：

//Problem: 3216  User: Uriel
//Memory: 552K  Time: 63MS
//Language: G++  Result: Accepted
//floyd+Bipartite Graph
//2010.08.02

#include<stdio.h>
#include<stdlib.h>
#include<algorithm>
using namespace std;
#define INF 1000000000

struct work{
    int p,t,d;
}b[210];

int Q,M;
int d[210][210],use[210],link[210],adj[210][210],n[210];

void floyd(){
    for(int k=1;k<=Q;++k)
        for(int i=1;i<=Q;++i)
            for(int j=1;j<=Q;++j)
                d[i][j]=min(d[i][j],d[i][k]+d[k][j]);
}

bool ok(int i,int j){
    if(b[i].t+b[i].d+d[b[i].p][b[j].p]<=b[j].t)return true;
    return false;
}

void Build_Graph(){
    int i,j,k;
    memset(adj,0,sizeof(adj));
    for(i=1;i<=M;i++){
        for(j=1;j<=M;j++){
            if(i==j)continue;
            if(ok(i,j))adj[i][++n[i]]=j;
        }
    }
}

int can(int t){
    int i,j;
    for(i=1;i<=n[t];++i){
        j=adj[t][i];
        if(!use[j]){
            use[j]=1;
            if(link[j]==-1 || can(link[j])){
                link[j]=t;
                return 1;
            }
        }
    }
    return 0;
}

int maxMatch(){
    int i,num=0;
    memset(link,-1,sizeof(link));
    for(i=1;i<=M;++i){
        memset(use,0,sizeof(use));
        if(can(i))num++;
    }
    return num;
}

int main(){
    int i,j,k;
    while(scanf("%d %d",&Q,&M),Q|M){
        for(i=1;i<=Q;++i){
            for(j=1;j<=Q;++j){
                scanf("%d",&d[i][j]);
                    if(d[i][j]==-1)d[i][j]=INF;
            }
        }
        floyd();
        for(i=1;i<=M;++i){
            scanf("%d %d %d",&b[i].p,&b[i].t,&b[i].d);
            n[i]=0;
        }
        Build_Graph();
        printf("%d\n",M-maxMatch());
    }
    return 0;
}

Uriel 2010-08-02 18:53 发表评论

POJ 1962 Corporative Network---并查集

Uriel — Sun, 01 Aug 2010 13:58:00 GMT

        很久没碰并查集，正好看见分类上还有一道就做了...
        菜啊，并查集稍微变化了一下就被搞死了...

题意：有两种操作：
            E x 询问节点x到他的根节点的长度
            I x y 将节点y成为x的父节点两点间长度为abs(x-y)%1000
            一直在想记录路径那里怎么改一下，搞来搞去还是不对，参考了CY的Blog之后恍然大悟... - -

            用结构体存每个节点，结构体里有两个变量，一个记录该节点的父节点，一个记录该节点到当前其父节点的距离
            每次I x y 操作之后将x的父节点设为y，x, y之间的距离即为abs(x-y)%1000
            每次E x 操作之后更新x的父节点，x到其父节点这一条线上的结点的父节点和与父节点的距离也依次更新

代码如下：

//Problem: 1962  User: Uriel
//Memory: 516K  Time: 610MS
//Language: G++  Result: Accepted
//Union-Find Set
//2010.08.01
#include<stdio.h>
#include<stdlib.h>

struct node{
    int f;
    int sum;
}p[20010];

int findset(int x){
    int i,ans=0;
    for(i=x;p[i].f!=i;i=p[i].f)ans+=p[i].sum;
    int s=ans;
    while(x!=i){
        int t=p[x].f,q=p[x].sum;
        p[x].f=i;
        p[x].sum=s;
        s-=q;
        x=t;
    }
    return ans;
}

int main(){
    int i,j,k,x,y,n,m;
    char str[10];
    scanf("%d",&m);
    while(m--){
        scanf("%d",&n);
        for(i=1;i<=n;i++){
            p[i].f=i;
            p[i].sum=0;
        }
        while(scanf("%s",str),str[0]!='O'){
            if(str[0]=='E'){
                scanf("%d",&x);
                printf("%d\n",findset(x));
            }
            else if(str[0]=='I'){
                scanf("%d %d",&x,&y);
                p[x].f=y;
                p[x].sum=abs(x-y)%1000;
            }
        }
    }
    return 0;
}

Uriel 2010-08-01 21:58 发表评论

POJ 图论---2

Uriel — Sat, 31 Jul 2010 14:34:00 GMT

刚开始学图论不久，这个是自己做过的一点图论水题，不一定全，有什么错误或者大家有其他该类的题的话欢迎提供~

Bipartite Graph：
POJ 1274
POJ 1422
POJ 1466
POJ 1469
POJ 1719
POJ 2060
POJ 2195 KM
POJ 2226
POJ 2239
POJ 2446
POJ 2516 KM
POJ 2536
POJ 2771
POJ 3041
POJ 3686 KM

Strongly Connected Components：
POJ 1236
POJ 1904
POJ 2186
POJ 2553
POJ 2762
POJ 3180

Uriel 2010-07-31 22:34 发表评论

POJ 图论---1

Uriel — Sat, 31 Jul 2010 14:26:00 GMT

刚开始学图论不久，这个是自己做过的一点图论水题，不一定全，有什么错误或者大家有其他该类的题的话欢迎提供~

Dijkstra && Floyd && Bellman_Ford && SPFA：
POJ 1062
POJ 1125
POJ 1130
POJ 1178
POJ 1237
POJ 1502
POJ 1511
POJ 1529
POJ 1797
POJ 1847
POJ 1860
POJ 1932
POJ 1975
POJ 2139
POJ 2240
POJ 2263
POJ 2387
POJ 2496
POJ 2502
POJ 2570
POJ 2594
POJ 3013
POJ 3259
POJ 3275
POJ 3311
POJ 3615
POJ 3653
POJ 3660
POJ 3767

Prim && Kruskal：
POJ 1251
POJ 1258
POJ 1287
POJ 1751
POJ 1789
POJ 1861
POJ 2031
POJ 2075
POJ 2253
POJ 2349
POJ 2377
POJ 2395
POJ 2421
POJ 2485
POJ 2560
POJ 3268
POJ 3723

System of Difference Constraints：
POJ 1201
POJ 1275
POJ 1364
POJ 2983
POJ 3159
POJ 3169
POJ 3687

Topological Sort：
POJ 1094
POJ 1128
POJ 1634
POJ 2367

Uriel 2010-07-31 22:26 发表评论

POJ 3180---The Cow Prom 强连通分量

Uriel — Sat, 24 Jul 2010 14:29:00 GMT

摘要: 题意依旧纠结，建图依旧纠结。。求含结点数>=2的连通分量数。。 PS：难得有图论题搜不到DY大牛的解题报告~ &nb... 阅读全文

Uriel 2010-07-24 22:29 发表评论

POJ 1904---King's Quest 强连通分量

Uriel — Sat, 24 Jul 2010 02:56:00 GMT

摘要: 最近在做图论，并且打算暑假军训之前(再次抱怨学校大二升大三暑假集训的制度...)除了每周5场组队赛之外大部分时间用来看图论，虽然也不知道来不来得及。。。强连通的问题看了好几天，因为组队赛占用很多时间，所以相关题也没做多少，依然觉得建图很困... 阅读全文

Uriel 2010-07-24 10:56 发表评论

POJ题目分类---二分法(Version 2010.07.18)

Uriel — Sun, 18 Jul 2010 10:17:00 GMT

评述的话主要是自己以后找起题来方便。。大牛们忽略就好。。

1064 不用__int64 ,int 就行，只能用C++交。。与浮点有关慎用G++啊。。

1186 Hash+二分

1434 几何相关。。

1505 二分+判可行

1650 二分枚举。。找最接近一个小数的分数。。半抄半写。。

1727 跟几何相关。。

1977 矩阵链乘，转移矩阵(A+I),每个baker状态为一行向量!

2018 二分+dp做的，还是O(n)的数学方法好

2153 排序加二分查找。。3922Ms水过。。

2366 hash不知道怎么做。。二分过了。。

2456 直接二分结果+判可行就行

3070 求斐波拉契数列某个数的后四位，题目已经给了转移矩阵，用二分求的。。搞了好一会儿0Ms，很开心。。结果得知不用二分也就16Ms

3104 二分+判断，注意时间k-1的转化

3122 二分面积

3233 矩阵链乘之和|A A|，右上角数即为所求，强大啊。。Orz !

|0 I |

3273 二分+判断

3320 快排+二分

3497 直接二分结果

3685 参考了别人的代码。。

3734 矩阵乘法。2 1 1 0 N次幂（用平方加速），左上角的数即为所求。

1 2 0 1

1 0 2 1

0 1 1 2

3735 矩阵乘法

3744 矩阵乘法。。初始矩阵：A=(0.0,1.0)，设下一个雷在x步以后，则将A乘以[0 1-p]的x次方（一定要用平方加速）。设最后得到了A=(a,b)，输出

[1 p ]
a*(1-p)+b*p即可

Uriel 2010-07-18 18:17 发表评论