C++博客-yaoyaozii-随笔分类-POJ

poj 1061 线性同余

Acaini.yaoyaozii — Wed, 31 Mar 2010 14:54:00 GMT

//推出同余式后直接上算法导论结论

#include <iostream>
using namespace std;
__int64 extgcd(__int64 a, __int64 b, __int64 &x, __int64 &y)
{
    if(b==0)
    {
        x=1,y=0;return a;
    }
    __int64 r=extgcd(b,a%b,x,y);
    __int64 t=x; x=y; y=t-a/b*y;
    return r;
}
int main()
{
    __int64 x,y,m,n,l,k,t;
    __int64 a,b,c,d;
    while(scanf("%I64d%I64d%I64d%I64d%I64d",&x,&y,&n,&m,&l)!=EOF)
    {

        a=m-n;
        b=x-y;
        if(m-n<0)
        {
            a=-a;
            b=(-b+l)%l;
        }
        else b=(b+l)%l;
        d=extgcd(a,l,k,t);
        if(b%d) { printf("Impossible\n"); continue;}
        k=(k*(b/d)%l+l)%(l/d);

        printf("%I64d\n", k);
    }
    return 0;
}

Acaini.yaoyaozii 2010-03-31 22:54 发表评论

poj 2115 线性同余

Acaini.yaoyaozii — Wed, 31 Mar 2010 14:48:00 GMT

A+C*X=B(%2^K)
C*X=B-A(%2^K)
令a=c,b=B-A,n=2^K;
利用以下结论（具体证明见《算法导论）：
推论1：方程ax=b(mod n)对于未知量x有解，当且仅当gcd(a,n) | b。
推论2：方程ax=b(mod n)或者对模n有d个不同的解，其中d=gcd(a,n)，或者无解。
定理1：设d=gcd(a,n)，假定对整数x和y满足d=ax+by(比如用扩展Euclid算法求出的一组解)。如果d | b，则方程ax=b(mod n)有一个解x0满足x0=x*(b/d) mod n 。特别的设e=x0+n，方程ax=b(mod n)的最小整数解x1=e mod (n/d)，最大整数解x2=x1+(d-1)*(n/d)。
定理2：假设方程ax=b(mod n)有解，且x0是方程的任意一个解，则该方程对模n恰有d个不同的解(d=gcd(a,n))，分别为：xi=x0+i*(n/d) mod n 。

a*x=b(%n) => a*x+n*y=b
d=ext_gcd(a,n,x0,y0)
最小整数解x1=(x0*(b/d)%n+n)%(n/d);

#include <iostream>
using namespace std;
__int64 exgcd(__int64 a, __int64 b, __int64 &x, __int64 &y)
{
    if(b==0)
    {
        x=1;y=0;return a;
    }
    __int64 r=exgcd(b, a%b, x, y);
    __int64 t=x;x=y;y=t-a/b*y;
    return r;
}
int main()
{
    __int64 A,B,C,K;
    __int64 a,b,n,d,x,y,e;
    while(scanf("%I64d%I64d%I64d%I64d",&A,&B,&C,&K)!=EOF)
    {
        if(A==0 && B==0 && C==0 && K==0) break;
        a=C;
        n=((__int64)1<<K); //小心溢出
        b=B-A;
        d=exgcd(a, n, x, y);
        if(b%d) {printf("FOREVER\n"); continue;}
        e=x*(b/d)%n+n;
        printf("%I64d\n",e%(n/d));
    }
    return 0;
};

Acaini.yaoyaozii 2010-03-31 22:48 发表评论

poj 3368 rmq(st)

Acaini.yaoyaozii — Sat, 27 Mar 2010 06:50:00 GMT

//离散化+st ?

#include <iostream>
#include <math.h>
using namespace std;

int a[100005];
int index[100005];
int f[100005];
int start[100005];
int mf[100005][100];
int n;

void rmq_init()
{
    int i,j,k=floor(log(double(n))/log(2.0));
    for(i=1; i<=n; i++) mf[i][0]=f[i];
    for(j=1; j<=k; j++)
        for(i=1; i+(1<<(j-1))<=n; i++)
            mf[i][j]=max(mf[i][j-1],mf[i+(1<<(j-1))][j-1]);
}
int rmq(int l,int r)
{
    int il,ir, m,k;
    il=index[a[l]],ir=index[a[r]];
    if(il==ir)
        return r-l+1;
    else
    if(il+1==ir)
        return max(start[ir]-l, r-start[ir]+1);
    else
    {
        m=max( start[il+1]-l, r-start[ir]+1);
        il++,ir--;
        k=floor(log(double(ir-il+1))/log(2.0));
        m=max( m, max(mf[il][k], mf[ir-(1<<k)+1][k]));
        return m;
    }
}
int main()
{
    int i,N,Q,l,r;
    while(scanf("%d", &N)!=EOF && N)
    {
        scanf("%d", &Q);
        memset(f, 0, sizeof(f));
        for(a[0]=0,i=1, n=0; i<=N; i++)
        {
            scanf("%d", a+i);
            if(a[i]!=a[i-1])
                index[ a[i] ]=++n, start[n]=i;
            f[n]++;
        }
        rmq_init();
        for(i=0; i<Q; i++)
        {
            scanf("%d%d", &l, &r);
            printf("%d\n", rmq(l,r));
        }
    }
    return 0;
}

Acaini.yaoyaozii 2010-03-27 14:50 发表评论

poj 3264 RMQ

Acaini.yaoyaozii — Fri, 26 Mar 2010 18:07:00 GMT

跑了8616K 2016MS，相当慢，不知道大牛怎么跑3、400MS左右

#include <iostream>
#include <math.h>
using namespace std;
int mx[50001][21];
int mi[50001][21];
int d[50001];
int n;
void rmq_init()
{
    int i,j;
    for(i=1;i<=n;i++) { mx[i][0]=d[i]; mi[i][0]=d[i]; }
    int m=floor(log((double)n)/log(2.0));
    for(i=1;i<=m;i++)
        for(j=0;j+(1<<(i-1))<=n;j++)
        {
            mx[j][i]=max(mx[j][i-1],mx[j+(1<<(i-1))][i-1]);
            mi[j][i]=min(mi[j][i-1],mi[j+(1<<(i-1))][i-1]);
        }
}

int rmq(int l,int r)
{
    int m=floor(log((double)(r-l+1))/log(2.0));
    int a=max(mx[l][m],mx[r-(1<<m)+1][m]);
    int b=min(mi[l][m],mi[r-(1<<m)+1][m]);
    return a-b;
}

int main()
{
    int q,l,r,i;
    while(scanf("%d%d", &n, &q)!=EOF)
    {
        for(i=1; i<=n; i++)
            scanf("%d", &d[i]);
        rmq_init();
        while(q--)
        {
            scanf("%d%d", &l, &r);
            printf("%d\n", rmq(l,r));
        }
    }
    return 0;
}

Acaini.yaoyaozii 2010-03-27 02:07 发表评论

poj 2352 树状数组

Acaini.yaoyaozii — Fri, 26 Mar 2010 14:07:00 GMT

//x从0开始，树状数组要求从一开始，故x++
//getsum(x)，求xi//level(sum)++,level为sum的star数++
#include <iostream>
using namespace std;
#define M 32001
int level[32005];//原数组
int c[32005];//树状数组
//k&(k^(k-1))
int getsum(int k)
{
    int sum;
    for(sum=0; k>0; k-=((-k)&k) ) sum+=c[k];
    return sum;
}
void modify(int k, int detal) //k:position ,detal:increase value
{
    for(; k<=M; k+=((-k)&k) ) c[k]+=detal;
}
int main()
{
    int x,y,n,i;
    while(scanf("%d", &n)!=EOF)
    {
        memset(c, 0, sizeof(c));
        memset(level, 0, sizeof(level));
        for(i=0; i<n; i++)
        {
            scanf("%d%d", &x, &y);
            level[ getsum(++x) ] ++;
            modify( x, 1);
        }
        for(i=0; i<n; i++)
        {
            printf("%d\n", level[i]);
        }
    }
    return 0;
}

Acaini.yaoyaozii 2010-03-26 22:07 发表评论

poj 1469 二分匹配

Acaini.yaoyaozii — Thu, 25 Mar 2010 08:29:00 GMT

#include <iostream>
#include <memory.h>
using namespace std;
#define P 101
#define M 301
bool map[P][M];
int n,m;
int chk[M];
int match[M];

int dfs(int p)
{
    int i,t;
    for(i=1; i<=m; i++)
        if(map[p][i] && !chk[i])
        {
            chk[i]=1;
            t=match[i];
            match[i]=p;
            if(t==0 || dfs(t)) return 1;
            match[i]=t;
        }
        return 0;
}

int mat()
{
    int i, res=0;
    for(i=1; i<=n; i++)
    {
        memset(chk, 0, sizeof(chk));
        if(dfs(i)) res++;
    }
    return res;
}
int main()
{
    int i,j,count,mi,kase;
    scanf("%d",&kase);
    while(kase--)
    {
        scanf("%d%d",&n, &m);
        memset(map, 0, sizeof(map));
        memset(match, 0, sizeof(match));
        for(i=1; i<=n; i++)
        {
            scanf("%d",&count);
            for(j=1; j<=count; j++)
            {
                scanf("%d", &mi);
                map[i][mi]=1; //map[mi][i]=true;
            }
        }
        if(mat()==n) printf("YES\n");
        else printf("NO\n");
    }
    return 0;
}

Acaini.yaoyaozii 2010-03-25 16:29 发表评论

poj 2186 有向图强连通分支

Acaini.yaoyaozii — Thu, 25 Mar 2010 08:28:00 GMT

#include <iostream>
using namespace std;
struct edge
{
    int v,next;
};
edge grap[60005];
edge grapt[61005];
int cnt;
int n,m;
bool visit[10005];
int finish[10005], id;
int grp[10005], g;
int outdegree[10005];
void add(int u, int v)
{
    grap[cnt].v=v;
    grap[cnt].next=grap[u].next;
    grap[u].next=cnt;
    grapt[cnt].v=u;
    grapt[cnt].next=grapt[v].next;
    grapt[v].next=cnt;
    cnt++;
}
void dfs(int u)
{
    int p;
    visit[u]=1;
    p=grap[u].next;
    while(p)
    {
        if(!visit[grap[p].v])
            dfs(grap[p].v);
        p=grap[p].next;
    }
    finish[++id]=u;
}
void dfst(int u)
{
    int p;
    visit[u]=1;
    p=grapt[u].next;
    grp[u]=g;
    while(p)
    {
        if(!visit[grapt[p].v])
            dfst(grapt[p].v);
        p=grapt[p].next;
    }
}
void scc()
{
    int i;
    g=0,id=0;
    memset(visit, 0, sizeof(visit));
    for(i=1; i<=n; i++)
        if(!visit[i])
            dfs(i);
    memset(visit, 0, sizeof(visit));
    for(i=n; i>0 ; --i)
    {
        if(!visit[finish[i]])
        {
            g++;
            dfst(finish[i]);
        }
    }
}
void solve()
{
    int i,j, p, tmp;
    memset(outdegree, 0, sizeof(outdegree));
    for(i=1; i<=n; i++)
    {
        p=grap[i].next;
        while(p)
        {
            if(grp[i]!=grp[ grap[p].v])
                ++outdegree[grp[i]];
            p=grap[p].next;
        }
    }
    for(i=1, tmp=0; i<=g; i++)
    {
        //printf("%d %d\n", i, outdegree[i]);
        if(outdegree[i] == 0)
        {
            ++tmp, j=i;
        }
    }
    //printf("%d %d\n", tmp, j);
    if(tmp!=1) // >1 || 0
    {
        printf("0\n");
    }
    else
    {
        tmp=0;
        for(i=1; i<=n; i++)
            if(grp[i]== j )
                tmp++;
        printf("%d\n", tmp);
    }
}
int main()
{
    int i, j, u, v;
    while(scanf("%d%d",&n, &m)!=EOF)
    {
        cnt=n+1;
        memset(grap, 0, sizeof(grap));
        memset(grapt, 0, sizeof(grapt));
        for(i=0; i<m; i++)
        {
            scanf("%d%d", &u, &v);
            add(u, v);
        }
        scc();
        solve();
    }
    return 0;
}

Acaini.yaoyaozii 2010-03-25 16:28 发表评论

poj 2239 二分图最大匹配

Acaini.yaoyaozii — Thu, 25 Mar 2010 08:26:00 GMT

#include <iostream>
using namespace std;
#define N 301
#define M 96
bool map[N][M];
int match[M];
bool chk[M];
int n, m=85;
int dfs(int p)
{
    int i,t;
    for(i=1; i<=m; i++)
    {
        if(!chk[i] && map[p][i] )
        {
            chk[i]=1;
            t=match[i];
            match[i]=p;
            if(t == -1 || dfs(t) ) return 1;
            match[i]=t;
        }
    }
    return 0;
}
int maxMah()
{
    int i, ans=0;
    for(i=1; i<=n; i++)
    {
        memset(chk, 0, sizeof(chk));
        if(dfs(i)) ans++;
    }
    return ans;
}

int main()
{
    int i,j,p,q,t;

    while(scanf("%d", &n)!=EOF)
    {
        //printf("%d\n",n);
        memset(map,0, sizeof(map));
        memset(match, -1, sizeof(match));
        for(i=1; i<=n; i++)
        {
            //printf("%d\n",i);
            scanf("%d", &t);
            //printf("%d\n", t);
            for(j=0; j<t; j++)
            {
                scanf("%d%d",&p, &q);
                map[i][(p-1)*12+q]=1;
            }
        }
        printf("%d\n",maxMah());
    }
    return 0;
}

Acaini.yaoyaozii 2010-03-25 16:26 发表评论

poj 3748 位运算

Acaini.yaoyaozii — Thu, 25 Mar 2010 07:50:00 GMT

#define setbit(x,y) x|=(1<<(y)) //将X的第Y位置1
#define clrbit(x,y) x&=~(1<<(y)) //将X的第Y位清0
...个人菜B做法：

#include <iostream>
using namespace std;
unsigned r,rp,x,y;
unsigned a[32]={0xfffffffe,0xfffffffd,0xfffffffb,0xfffffff7,
                0xffffffef,0xffffffdf,0xffffffbf,0xffffff7f,
                0xfffffeff,0xfffffdff,0xfffffbff,0xfffff7ff,
                0xffffefff,0xffffdfff,0xffffbfff,0xffff7fff,
                0xfffeffff,0xfffdffff,0xfffbffff,0xfff7ffff,
                0xffefffff,0xffdfffff,0xffbfffff,0xff7fffff,
                0xfeffffff,0xfdffffff,0xfbffffff,0xf7ffffff,
                0xefffffff,0xdfffffff,0xbfffffff,0x7fffffff,};
int main()
{
    while(scanf("%x,%d,%d", &r, &x, &y)!=EOF)
    {
        printf("%x\n",(((r&a[x])|(1<<y))|(1<<(y-1)))&a[y-2]);
    }
    return 0;
}

Acaini.yaoyaozii 2010-03-25 15:50 发表评论

poj 2724 二分匹配

Acaini.yaoyaozii — Wed, 24 Mar 2010 17:23:00 GMT

#include <iostream>
using namespace std;
#define N 1100
#define M 1100
bool map[N][M];
int match[M];
bool chk[M];
int n, m;
bool num[1100];
int hash[1100];
int dfs(int p)
{
    int i,t;
    for(i=1; i<=m; i++)
    {
        if(!chk[i] && map[p][i] )
        {
            chk[i]=1;
            t=match[i];
            match[i]=p;
            if(t == -1 || dfs(t) ) return 1;
            match[i]=t;
        }
    }
    return 0;
}
int maxMah()
{
    int i, ans=0;
    for(i=1; i<=n; i++)
    {
        memset(chk, 0, sizeof(chk));
        if(dfs(i)) ans++;
    }
    return ans;
}

int getnum(char *s)
{
    int i,ret=0, l=strlen(s);
    for(i=0; s[i]; i++)
    {
        ret+= (s[i]=='1'?(1<<(l-i-1)):0);
    }
    return ret;
}
int main()
{
    int x,y, i, j, k, c;//,ans=0;
    char tmp[15];
    while(scanf("%d%d", &x, &y)!=EOF && x!=0 && y!=0 )
    {
        memset(map, 0, sizeof(map));
        memset(num,0, sizeof(num));
        memset(hash,0, sizeof(hash));
        memset(match, -1, sizeof(match));
        for(i=0; i<y; i++)
        {
            scanf("%s",tmp);
            for(j=0; tmp[j]; j++)
                if(tmp[j]=='*')
                    break;
            if(tmp[j])
            {
                tmp[j]='0';
                num[ getnum(tmp) ] =true ;
                tmp[j]='1';
                num[ getnum(tmp) ] =true ;
            }
            else num[ getnum(tmp) ] =true ;
        }
        for(i=0, k=0; i<(1<<x); i++)
            if(num[i])
            {
                hash[++k]=i;
            }

        for(i=1; i<=k; i++)
            for(j=1; j<=k; j++)
            {

                c=hash[i]^hash[j];//判是否相差一位
                if(c && ((c&(c-1))==0) )//判是否相差一位
                    map[i][j]=map[j][i]=1;
            }
        n=m=k;
        printf("%d\n",k-maxMah()/2); //(k-maxMah)+maxMah/2=k-maxMah()/2
    }
    return 0;
}

Acaini.yaoyaozii 2010-03-25 01:23 发表评论

pku poj 2954[pick+gcd+叉积]

Acaini.yaoyaozii — Tue, 06 Oct 2009 10:28:00 GMT

#include<iostream>
using namespace std;
int x[4];
int y[4];
int gcd(int a,int b)
{
    if(b==0) return a;
    else return gcd(b,a%b);
}
int main()
{
    int area;
    int b,i;
    while(scanf("%d%d%d%d%d%d",x+1,y+1,x+2,y+2,x+3,y+3)!=EOF)
    {
        if(x[1]==0 && y[1]==0
            && x[2]==0 && y[2]==0
            && x[3]==0 && y[3]==0 ) break;
        area=abs( (x[2]-x[1])*(y[3]-y[1])-(x[3]-x[1])*(y[2]-y[1]) );
        b=gcd(abs(x[2]-x[1]),abs(y[2]-y[1]))+gcd(abs(x[3]-x[2]),abs(y[3]-y[2]))+gcd(abs(x[1]-x[3]),abs(y[1]-y[3]));
        i=(area+2-b)>>1;
        printf("%d\n",i);

    }
    return 0;
}

Acaini.yaoyaozii 2009-10-06 18:28 发表评论

poj 1265[计算几何]

Acaini.yaoyaozii — Tue, 06 Oct 2009 09:58:00 GMT

学习了pick定理
给定顶点座标均是整点（或正方形格点）的简单多边形，皮克定理说明了其面积A和内部格点数目i、边上格点数目b的关系：A = i + b/2 - 1。

#include<iostream>
#include<math.h>
using namespace std;

#define MAXN 500

struct point{int x,y;};
point p[MAXN];
int gcd(int a,int b)
{
    if(b==0) return a;
    else return gcd(b,a%b);
}
//计算多边形面积,顶点按顺时针或逆时针给出
double area_polygon(int n,point* p){
    double s1=0,s2=0;
    int i;
    for (i=0;i<n;i++)
        s1+=p[(i+1)%n].y*p[i].x,s2+=p[(i+1)%n].y*p[(i+2)%n].x;
    return fabs(s1-s2)/2;
}

int main()
{
    int tc,m,I,E,dx,dy,i,j,k;
    scanf("%d",&tc);
    for(k=1;k<=tc;k++)
    {
        scanf("%d",&m);
        p[0].x=0;
        p[0].y=0;
        E=0;
        for(i=1;i<=m;i++)
        {
            scanf("%d%d",&dx,&dy);
            p[i].x=p[i-1].x+dx;
            p[i].y=p[i-1].y+dy;
            E+=gcd(abs(dx),abs(dy));
        }
        double area=area_polygon(m+1,p);
        I=(2*area+2-E)/2;
        printf("Scenario #%d:\n%d %d %.1lf\n\n",k,I,E,area);
    }
    return 0;
}

Acaini.yaoyaozii 2009-10-06 17:58 发表评论

pku 1654[计算几何]

Acaini.yaoyaozii — Tue, 06 Oct 2009 05:12:00 GMT

500MS弄到79MS，求多边形面积，有可能是凹，用叉积就可以了，不能存所有点，会MLE
一次读入比单个读入快

#include<iostream>
#include<math.h>
using namespace std;
struct point {__int64 x,y;};
__int64 dx[]={0,-1,0,1,-1,0,1,-1,0,1};
__int64 dy[]={0,-1,-1,-1,0,0,0,1,1,1};

//__int64 xmult(point p0,point p1,point p2){
//    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
//}
char path[1000010];
int main()
{
    int t;
    point p0,p1,p2;
    __int64 ans;
    scanf("%d",&t);
    while(t--)
    {
        p0.x=0;
        p0.y=0;
        p2=p1=p0;
        ans=0;
        scanf(" %s",path);
        int i=0;
        for(i=0;path[i]!='5';i++)
        {
            p1=p2;
            p2.x+=dx[path[i]-'0'];
            p2.y+=dy[path[i]-'0'];
            ans+=(p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
        }
        if(ans<0) ans*=-1;
        printf("%I64d",(ans>>1));
        if(((ans>>1)<<1)==ans) printf("\n");
        else printf(".5\n");
    }
    return 0;
}

Acaini.yaoyaozii 2009-10-06 13:12 发表评论

pku 1066[计算几何]

Acaini.yaoyaozii — Mon, 05 Oct 2009 07:05:00 GMT

//枚举4条边界的点a，然后求线段out(trea,a)与其他wall相交的最小数
//代码中的max改为min有助理解

#include<iostream>
using namespace std;
#define MAXN 30
#define eps 1e-8
#define zero(x) (((x)>0?(x):-(x))

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

line wall[MAXN+1];
int n;
point trea;
//计算cross product (P1-P0)x(P2-P0)
double xmult(point p1,point p2,point p0){
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}

//判两点在线段异侧,点在线段上返回0
int opposite_side(point p1,point p2,line l){
    return xmult(l.a,p1,l.b)*xmult(l.a,p2,l.b)<-eps;
}
//判两线段相交,不包括端点和部分重合
int intersect_ex(line u,line v){
    return opposite_side(u.a,u.b,v)&&opposite_side(v.a,v.b,u);
}
int main()
{
    int i;
    scanf("%d",&n);
    for(i=0;i<n;i++)
        scanf("%lf%lf%lf%lf",&(wall[i].a.x),&(wall[i].a.y),&(wall[i].b.x),&(wall[i].b.y));
    scanf("%lf%lf",&(trea.x),&(trea.y));
    line out;
    out.a=trea;
    int max=n+1;
    double x,y;
    for(x=0,out.b.y=100;x<=100;x+=0.01)
    {
        out.b.x=x;
        int t=0;
        for(i=0;i<n;i++)
            if(intersect_ex(out,wall[i]))
                t++;
        if(t<max) max=t;
    }
    for(x=0,out.b.y=0;x<=100;x+=0.01)
    {
        out.b.x=x;
        int t=0;
        for(i=0;i<n;i++)
            if(intersect_ex(out,wall[i]))
                t++;
        if(t<max) max=t;
    }
    for(y=0,out.b.x=0;y<=100;y+=0.01)
    {
        out.b.y=y;
        int t=0;
        for(i=0;i<n;i++)
            if(intersect_ex(out,wall[i]))
                t++;
        if(t<max) max=t;
    }
    for(y=0,out.b.x=100;y<=100;y+=0.01)
    {
        out.b.y=y;
        int t=0;
        for(i=0;i<n;i++)
            if(intersect_ex(out,wall[i]))
                t++;
        if(t<max) max=t;
    }
    printf("Number of doors = %d\n",max+1);
    return 0;
}

Acaini.yaoyaozii 2009-10-05 15:05 发表评论

pku 1696[计算几何]

Acaini.yaoyaozii — Mon, 05 Oct 2009 04:12:00 GMT

#include<stdio.h>
#include<math.h>
#include<memory.h>
//using namespace std;
#define eps 1e-8
#define MAXN 50
#define PI 3.1415926535898
struct point{ double x,y;};

point planet[MAXN+1];
bool visit[MAXN+1];
int ans[MAXN+1];

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

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

int main()
{
    int t,indice,i,j,k,n;
    double max;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(i=1;i<=n;i++)
            scanf("%d%lf%lf",&indice,&(planet[i].x),&(planet[i].y));

        memset(visit,0,sizeof(visit));

        int next=1;

        for(i=2;i<=n;i++)
            if(planet[next].y>planet[i].y)
                next=i;
            else if(planet[next].y == planet[i].y)
                if(planet[next].x>planet[i].x)
                    next=i;
        point p1,p0;
        p1.x=0;
        p1.y=planet[next].y;
        p0.x=planet[next].x;
        p0.y=planet[next].y;
        k=0;
        double thta,max;
        for(j=1;j<=n;j++)
        {
            visit[next]=true;
            ans[k++]=next;
            max=2*PI;
            for(i=1;i<=n ; i++)
            {
                if(visit[i]) continue;
                double xmul=xmult(p0,p1,planet[i]);
                double dmul=dmult(p0,p1,planet[i]);
                if( xmul<=0)
                {
                    if(dmul==0 )
                        thta=3*PI/2.0;
                    else
                        if(dmul<0 ) thta=PI+atan(xmul/dmul);
                        else  thta=2*PI+atan(xmul/dmul);

                }
                if( thta<max ||
                    (fabs(thta-max)<eps && fabs(distance(p0,planet[i])-distance(p0,planet[next]))<eps)
                    )
                { next=i;max=thta;}
            }
            p1.x=p0.x;
            p1.y=p0.y;
            p0.x=planet[next].x;
            p0.y=planet[next].y;
        }
        printf("%d",k);
        for(i=0;i<k;i++)
            printf(" %d",ans[i]);
        printf("\n");
    }
    return 0;
}

Acaini.yaoyaozii 2009-10-05 12:12 发表评论

pku 2653[计算几何]

Acaini.yaoyaozii — Sun, 04 Oct 2009 13:39:00 GMT

// 同样是暴力，不同的暴力，一个TLE，一个5**MS
//判线段相交

#include<iostream>
#include <math.h>
using namespace std;
#define MAXN 100000
#define eps 1e-8
#define zero(x) (((x)>0?(x):-(x))

struct point{double x,y;};
struct line{ point a,b;};
line sti[MAXN+1];
bool under[MAXN+1];

//计算cross product (P1-P0)x(P2-P0)
double xmult(point p1,point p2,point p0){
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}

int dots_inline(point p1,point p2,point p3){
    return zero(xmult(p1,p2,p3));
}
int same_side(point p1,point p2,line l){
    return xmult(l.a,p1,l.b)*xmult(l.a,p2,l.b)>eps;
}
int dot_online_in(point p,line l){
    return zero(xmult(p,l.a,l.b))&&(l.a.x-p.x)*(l.b.x-p.x)<eps&&(l.a.y-p.y)*(l.b.y-p.y)<eps;
}

int intersect_in(line u,line v){
    if (!dots_inline(u.a,u.b,v.a)||!dots_inline(u.a,u.b,v.b))
        return !same_side(u.a,u.b,v)&&!same_side(v.a,v.b,u);
    return dot_online_in(u.a,v)||dot_online_in(u.b,v)||dot_online_in(v.a,u)||dot_online_in(v.b,u);
}

int main()
{
    int n,i,j;
    while(scanf("%d",&n)!=EOF && n)
    {
        memset(under,0,sizeof(under));
        for(i=1;i<=n;i++)
            scanf("%lf%lf%lf%lf",&(sti[i].a.x),&(sti[i].a.y),&(sti[i].b.x),&(sti[i].b.y));
        for(i=1;i<n;i++)
            for(j=i+1;j<=n;j++)
            {
                if(intersect_in(sti[i],sti[j]))
                {
                    under[i]=true;
                    break;
                }
            }
        printf("Top sticks:");
        for(i=1;i<n;i++)
            if(!under[i])
                printf(" %d,",i);
        printf(" %d.\n",n);
    }
    return 0;
}

Acaini.yaoyaozii 2009-10-04 21:39 发表评论

pku 1296[计算几何]

Acaini.yaoyaozii — Sun, 04 Oct 2009 10:15:00 GMT

//学习判定两直线关系：平行重合相交
//求交点

#include<iostream>
#include <math.h>
using namespace std;
#define eps 1e-8
#define zero(x) (((x)>0?(x):-(x))
struct point{double x,y;};
struct line{point a,b;};

double xmult(line u,line v){
    return (u.b.x-u.a.x)*(v.b.y-v.a.y)-(v.b.x-v.a.x)*(u.b.y-u.a.y);
}
//计算两直线交点,注意事先判断直线是否平行!
point intersection(line u,line v){
    point ret=u.a;
    double t=((u.a.x-v.a.x)*(v.a.y-v.b.y)-(u.a.y-v.a.y)*(v.a.x-v.b.x))
            /((u.a.x-u.b.x)*(v.a.y-v.b.y)-(u.a.y-u.b.y)*(v.a.x-v.b.x));
    ret.x+=(u.b.x-u.a.x)*t;
    ret.y+=(u.b.y-u.a.y)*t;
    return ret;
}

int main()
{
    line u,v;
    int n;
    scanf("%d",&n);
    printf("INTERSECTING LINES OUTPUT\n");
    while(n--)
    {
        scanf("%lf%lf%lf%lf",&(u.a.x),&(u.a.y),&(u.b.x),&(u.b.y));
        scanf("%lf%lf%lf%lf",&(v.a.x),&(v.a.y),&(v.b.x),&(v.b.y));
        double a1,b1,c1,a2,b2,c2;
        //xishu
        a1=u.a.y-u.b.y;
        b1=-(u.a.x-u.b.x);
        c1=u.a.x*u.b.y-u.a.y*u.b.x;

        a2=v.a.y-v.b.y;
        b2=-(v.a.x-v.b.x);
        c2=v.a.x*v.b.y-v.a.y*v.b.x;
        //if( fabs(xmult(u,v))
        if(fabs(a1*b2-a2*b1)<eps)
        {
            if( fabs(a1*c2-a2*c1)<eps && fabs(b1*c2-b2*c1)<eps )
                printf("LINE\n");
            else printf("NONE\n");
        }
        else
        {
            point res;
            //res=intersection(u,v);
            res.x=(b1*c2-b2*c1)/(a1*b2-a2*b1);
            res.y=(a2*c1-a1*c2)/(a1*b2-a2*b1);
            printf("POINT %.2lf %.2lf\n",res.x,res.y);
        }
    }
    printf("END OF OUTPUT\n");
    return 0;
}

Acaini.yaoyaozii 2009-10-04 18:15 发表评论

pku 2398[计算几何]

Acaini.yaoyaozii — Sun, 04 Oct 2009 05:12:00 GMT

//和2318差不多，只是输入的board是无序的，输出是玩具数为t(>0)的箱子数(>0)

#include<iostream>
#include<stdlib.h>
using namespace std;
#define MAXN 5000+1
#define eps 1e-8
#define _sign(x) ((x)>eps?1:((x)<-eps?2:0))

struct point{int x,y;};
struct line{point a,b;};

int ans[MAXN];
line board[MAXN];

double xmult(point p0,point p1,point p2){//叉积，囧
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
int cmp(const void *a,const void *b)
{
    line *aa=(line*)a,*bb=(line*)b;
    if(aa->a.x == bb->a.x)
        return aa->b.x - bb->b.y;
    else return aa->a.x - bb->a.x;
}
int main()
{
    int n,m,x1,y1,x2,y2,u1,u2,i;
    point toys;
    while(scanf("%d",&n)!=EOF && n )
    {
        memset(ans,0,sizeof(ans));
        scanf("%d%d%d%d%d",&m,&x1,&y1,&x2,&y2);
        for(i=0;i<n;i++)
        {
            scanf("%d%d",&u1,&u2);
            board[i].a.x=u1;
            board[i].a.y=y1;
            board[i].b.x=u2;
            board[i].b.y=y2;
        }
        qsort(board,n,sizeof(board[0]),cmp);
        board[n].a.x=x2;
        board[n].a.y=y1;
        board[n].b.x=x2;
        board[n].b.y=y2;
        for(i=0;i<m;i++)
        {
            scanf("%d%d",&(toys.x),&(toys.y));
            int low=0,hig=n,mid;
            while(low+1<hig)
            {
                mid=(low+hig)>>1;
                double res=xmult(board[mid].a,board[mid].b,toys);
                if(res>0) low=mid;
                else hig=mid;
            }
            if(xmult(board[low].a,board[low].b,toys)<0) ans[low]++;
            else  ans[low+1]++;
        }
        printf("Box\n");
        int t=1,c;
        for(t=1;t<=m;t++)
        {
            for(i=0,c=0;i<=n;i++)
            {
                if(ans[i]==t)
                {
                    c++;
                    ans[i]=0;
                }
            }
            if(c>0)
                printf("%d: %d\n",t,c);
            m=m-c*t;
        }

    }
    return 0;
}

Acaini.yaoyaozii 2009-10-04 13:12 发表评论

【转】POJ 计算几何入门题目推荐

Acaini.yaoyaozii — Sun, 04 Oct 2009 04:18:00 GMT

//今天开始按着这个来做题训练计算几何了

计算几何题的特点与做题要领：
1.大部分不会很难，少部分题目思路很巧妙
2.做计算几何题目，模板很重要，模板必须高度可靠。
3.要注意代码的组织，因为计算几何的题目很容易上两百行代码，里面大部分是模板。如果代码一片混乱，那么会严重影响做题正确率。
4.注意精度控制。
5.能用整数的地方尽量用整数，要想到扩大数据的方法（扩大一倍，或扩大sqrt2）。因为整数不用考虑浮点误差，而且运算比浮点快。

一。点，线，面，形基本关系，点积叉积的理解

POJ 2318 TOYS（推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=2318
POJ 2398 Toy Storage（推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=2398
一个矩形，有被若干直线分成N个格子，给出一个点的坐标，问你该点位于哪个点中。
知识点：其实就是点在凸四边形内的判断，若利用叉积的性质，可以二分求解。

POJ 3304 Segments
http://acm.pku.edu.cn/JudgeOnline/problem?id=3304
知识点：线段与直线相交，注意枚举时重合点的处理

POJ 1269 Intersecting Lines
http://acm.pku.edu.cn/JudgeOnline/problem?id=1269
知识点：直线相交判断，求相交交点

POJ 1556 The Doors （推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=1556
知识点：简单图论＋简单计算几何，先求线段相交，然后再用Dij求最短路。

POJ 2653 Pick-up sticks
http://acm.pku.edu.cn/JudgeOnline/problem?id=2653
知识点：还是线段相交判断

POJ 1066 Treasure Hunt
http://acm.pku.edu.cn/JudgeOnline/problem?id=1066
知识点：线段相交判断，不过必须先理解“走最少的门”是怎么一回事。

POJ 1410 Intersection
http://acm.pku.edu.cn/JudgeOnline/problem?id=1410
知识点：线段与矩形相交。正确理解题意中相交的定义。
详见：http://hi.baidu.com/novosbirsk/blog/item/68c682c67e8d1f1d9d163df0.html

POJ 1696 Space Ant （推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=1696
德黑兰赛区的好题目。需要理解点积叉积的性质

POJ 3347 Kadj Squares
http://acm.pku.edu.cn/JudgeOnline/problem?id=3347
本人的方法极度猥琐。复杂的线段相交问题。这个题目是计算几何的扩大数据运算的典型应用，扩大根号2倍之后就避免了小数。

POJ 2826 An Easy Problem?! （推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=2826
问：两条直线组成一个图形，能容纳多少雨水。很不简单的Easy Problem，要考虑所有情况。你不看discuss看看能否AC。（本人基本不能）提示一下，水是从天空垂直落下的。

POJ 1039 Pipe
http://acm.pku.edu.cn/JudgeOnline/problem?id=1039
又是线段与直线相交的判断，再加上枚举的思想即可。

POJ 3449 Geometric Shapes
http://acm.pku.edu.cn/JudgeOnline/problem?id=3449
判断几何体是否相交，不过输入输出很恶心。
此外，还有一个知识点，就是给出一个正方形（边不与轴平行）的两个对角线上的顶点，需要你求出另外两个点。必须掌握其方法。

POJ 1584 A Round Peg in a Ground Hole
http://acm.pku.edu.cn/JudgeOnline/problem?id=1584
知识点：点到直线距离，圆与多边形相交，多边形是否为凸

POJ 2074 Line of Sight （推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=2074
与视线问题的解法，关键是求过两点的直线方程，以及直线与线段的交点。数据有一个trick，要小心。

二。凸包问题

POJ 1113 Wall
http://acm.pku.edu.cn/JudgeOnline/problem?id=1113
知识点：赤裸裸的凸包问题，凸包周长加上圆周。

POJ 2007 Scrambled Polygon
http://acm.pku.edu.cn/JudgeOnline/problem?id=2007
知识点：凸包，按极角序输出方案

POJ 1873 The Fortified Forest （推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=1873
World Final的水题，先求凸包，然后再搜索。由于规模不大，可以使用位运算枚举。
详见：http://hi.baidu.com/novosbirsk/blog/item/333abd54c7f22c52574e0067.html

POJ 1228 Grandpa's Estate （推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=1228
求凸包顶点数目，很多人求凸包的模板是会多出点的，虽然求面积时能得到正确答案，但是在这个题目就会出问题。此外，还要正确理解凸包的性质。

POJ 3348 Cows
http://acm.pku.edu.cn/JudgeOnline/problem?id=3348
凸包面积计算

三。面积问题，公式问题

POJ 1654 Area
http://acm.pku.edu.cn/JudgeOnline/problem?id=1654
知识点：利用有向面积（叉积）计算多边形面积

POJ 1265 Area
http://acm.pku.edu.cn/JudgeOnline/problem?id=1265
POJ 2954 Triangle
http://acm.pku.edu.cn/JudgeOnline/problem?id=2954
Pick公式的应用，多边形与整点的关系。（存在一个GCD的关系）

四。半平面交

半平面交的主要应用是判断多边形是否存在核，还可以解决一些与线性方程组可行区域相关的问题（就是高中时的那些）。

POJ 3335 Rotating Scoreboard
http://acm.pku.edu.cn/JudgeOnline/problem?id=3335
POJ 3130 How I Mathematician Wonder What You Are!
http://acm.pku.edu.cn/JudgeOnline/problem?id=3130
POJ 1474 Video Surveillance
http://acm.pku.edu.cn/JudgeOnline/problem?id=1474
知识点：半平面交求多边形的核，存在性判断

POJ 1279 Art Gallery
http://acm.pku.edu.cn/JudgeOnline/problem?id=1279
半平面交求多边形的核，求核的面积

POJ 3525 Most Distant Point from the Sea （推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=3525
给出一个多边形，求里面的一个点，其距离离多边形的边界最远，也就是多边形中最大半径圆。
可以使用半平面交+二分法解。二分这个距离，边向内逼近，直到达到精度。

POJ 3384 Feng Shui （推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=3384
半平面交实际应用，用两个圆覆盖一个多边形，问最多能覆盖多边形的面积。
解法：用半平面交将多边形的每条边一起向“内”推进R，得到新的多边形，然后求多边形的最远两点。

POJ 1755 Triathlon （推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=1755
半平面交判断不等式是否有解。注意不等式在转化时正负号的选择，这直接影响到半平面交的方向。

POJ 2540 Hotter Colder
http://acm.pku.edu.cn/JudgeOnline/problem?id=2540
半平面交求线性规划可行区域的面积。

POJ 2451 Uyuw's Concert
http://acm.pku.edu.cn/JudgeOnline/problem?id=2451
Zzy专为他那篇nlogn算法解决半平面交问题的论文而出的题目。

五。计算几何背景，实际上解题的关键是其他问题（数据结构、组合数学，或者是枚举思想）
若干道经典的离散化＋扫描线的题目，ACM选手必做题目

POJ 1151 Atlantis （推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=1151
POJ 1389 Area of Simple Polygons
http://acm.pku.edu.cn/JudgeOnline/problem?id=1389
矩形离散化，线段树处理，矩形面积求交

POJ 1177 Picture （推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=1177
矩形离散化，线段树处理，矩形交的周长，这个题目的数据比较强。线段树必须高效。

POJ 3565 Ants （推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=3565
计算几何中的调整思想，有点像排序。要用到线段相交的判断。
详见：http://hi.baidu.com/novosbirsk/blog/item/fb668cf0f362bec47931aae2.html

POJ 3695 Rectangles
http://acm.pku.edu.cn/JudgeOnline/problem?id=3695
又是矩形交的面积，但是由于是多次查询，而且矩形不多，使用组合数学中的容斥原理解决之最适合。线段树是通法，但是除了线段树，还有其他可行的方法。

POJ 2002 Squares
http://acm.pku.edu.cn/JudgeOnline/problem?id=2002
枚举思想，求平面上若干个点最多能组成多少个正方形，点的Hash

POJ 1434 Fill the Cisterns!（推荐）
http://acm.pku.edu.cn/JudgeOnline/problem?id=1434
一开始发昏了，准备弄个线段树。其实只是个简单的二分。

六。随机算法
POJ 2420 A Star not a Tree?
http://acm.pku.edu.cn/JudgeOnline/problem?id=2420
多边形的费马点。所谓费马点，就是多边形中一个点P，该点到其他点的距离之和最短。四边形以上的多边形没有公式求费马点，因此可以使用随机化变步长贪心法。
详见：http://hi.baidu.com/novosbirsk/blog/item/75983f138499f825dd54019b.html

七。解析几何
这种题目本人不擅长，所以做得不多，模板很重要。当然，熟练运用叉积、点积的性质还是很有用的。
POJ 1375 Intervals
http://acm.pku.edu.cn/JudgeOnline/problem?id=1375
知识点：过圆外一点求与圆的切线

POJ 1329 Circle Through Three Points
http://acm.pku.edu.cn/JudgeOnline/problem?id=1329
求三角形外接圆

POJ 2354 Titanic
http://acm.pku.edu.cn/JudgeOnline/problem?id=2354
求球面上两个点的距离，而且给的是地理经纬坐标。

POJ 1106 Transmitters
http://acm.pku.edu.cn/JudgeOnline/problem?id=1106
角度排序，知道斜率求角度，使用atan函数。

POJ 1673 EXOCENTER OF A TRIANGLE
http://acm.pku.edu.cn/JudgeOnline/problem?id=1673
可以转化为三角形的垂心问题。

八。旋转卡壳

POJ 2187 Beauty Contest
http://acm.pku.edu.cn/JudgeOnline/problem?id=2187
凸包求最远点对。可以暴力枚举，也可以使用旋转卡壳。

POJ 3608 Bridge Across Islands（难）
http://acm.pku.edu.cn/JudgeOnline/problem?id=3608
两个凸包的最近距离。本人的卡壳始终WA。郁闷。

九。其他问题
POJ 1981 Circle and Points
http://acm.pku.edu.cn/JudgeOnline/problem?id=1981
求单位圆最多能覆盖平面上多少个点

Acaini.yaoyaozii 2009-10-04 12:18 发表评论

pku 2318[计算几何]

Acaini.yaoyaozii — Sun, 04 Oct 2009 04:14:00 GMT

//计算几何第一题
//学习了叉积和利用叉积判断左右位置关系
//叉积+二分
//延续了这两天比赛的PE黑手

#include<iostream>
using namespace std;
#define MAXN 5000+1
#define eps 1e-8
#define _sign(x) ((x)>eps?1:((x)<-eps?2:0))

struct point{int x,y;};
struct line{point a,b;};

int ans[MAXN];
point toys[MAXN];
line board[MAXN];

double xmult(point p0,point p1,point p2){//叉积，囧
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
int main()
{
    int n,m,x1,y1,x2,y2,u1,u2,i;
    bool flag=false;
    while(scanf("%d",&n)!=EOF && n )
    {
        if(flag) printf("\n");
        memset(ans,0,sizeof(ans));
        scanf("%d%d%d%d%d",&m,&x1,&y1,&x2,&y2);
        for(i=0;i<n;i++)
        {
            scanf("%d%d",&u1,&u2);
            board[i].a.x=u1;
            board[i].a.y=y1;
            board[i].b.x=u2;
            board[i].b.y=y2;
        }
        board[n].a.x=x2;
        board[n].a.y=y1;
        board[n].b.x=x2;
        board[n].b.y=y2;
        for(i=0;i<m;i++)
        {
            scanf("%d%d",&(toys[i].x),&(toys[i].y));
            int low=0,hig=n,mid;
            while(low+1<hig)
            {
                mid=(low+hig)>>1;
                double res=xmult(board[mid].a,board[mid].b,toys[i]);
                //if(res>0)
                if(res>0) low=mid;
                else hig=mid;
            }
            if(xmult(board[low].a,board[low].b,toys[i])<0) ans[low]++;
            else  ans[low+1]++;
        }
        for(i=0;i<=n;i++)
        {
            printf("%d: %d\n",i,ans[i]);
        }
        flag=true;
    }
    return 0;
}

Acaini.yaoyaozii 2009-10-04 12:14 发表评论

pku 2406

Acaini.yaoyaozii — Mon, 03 Aug 2009 06:25:00 GMT

传说可以用kmp的精髓next[]来做，但是怎么都领悟不了，所以就凭感觉直接模拟了
因为s=a^n
取前缀，分别从长度1到strlen，然后与后面剩下的比较，如果ok就说明这个长度是a的最短的长度，然后sterlen(s)/strlen(a)

#include<iostream>
using namespace std;
#define M 1000000

char t[M+1],p[M+1];
int lent,lenp;

bool kmp(char *t,char *p)
{
    int i,j;
    for(i=lenp,j=0;i<lent;i++)
    {
        if(t[i]!=p[j]) return false;
        if(t[i]==p[j]) j++;
        if(j==lenp) j=0;
    }
    return true;
}

int main()
{
    int i;
    while(scanf(" %s",&t)!=EOF && t[0]!='.')
    {
        lent=strlen(t);
        for(i=0;i<lent;i++)
        {
            p[i]=t[i];
            p[i+1]='\0';
            lenp=i+1;
            if(lent%lenp==0 && kmp(t,p))
            {
                printf("%d\n",lent/lenp);
                break;
            }
        }
    }
    return 0;
}

Acaini.yaoyaozii 2009-08-03 14:25 发表评论

pku 2031

Acaini.yaoyaozii — Wed, 29 Jul 2009 02:05:00 GMT

#include<iostream>
#include<cmath>
using namespace std;
struct BALL
{
    double x,y,z,r;
};
struct Edge
{
    int x,y;
    double w;
};
BALL ball[200];
Edge edge[10000];
int parent[200];
int n;

int find(int c)
{
    if(parent[c]<0) return c;
    else return find(parent[c]);
}

bool uni(int x,int y)
{
    int a,b,t;
    a=find(x);
    b=find(y);
    if(a!=b)
    {
        t=parent[a]+parent[b];
        if(parent[a]<parent[b])
        {
            parent[b]=a;
            parent[a]=t;
        }
        else
        {
            parent[a]=b;
            parent[b]=t;
        }
        return true;
    }
    return false;
}
double Kruskal(int en)
{
    int i;
    double mst=0;
    memset(parent,-1,sizeof(parent));
    for(i=0;i<en;i++)
        if(uni(edge[i].x,edge[i].y))
            mst+=edge[i].w;
    return mst;
}
int cmp(const void *a, const void *b)
{
    double ta=(*(Edge*)a).w,tb=(*(Edge*)b).w;
    if(ta<tb) return -1;
    else if(ta==tb) return 0;
    else return 1;
}

int main()
{
    int i,j,k;
    double dx,dy,dz,ts;
    while(scanf("%d",&n)!=EOF && n)
    {


        for(i=0;i<n;i++)
            scanf("%lf%lf%lf%lf",&(ball[i].x),&(ball[i].y),&(ball[i].z),&(ball[i].r));
        for(i=0,k=0;i<n;i++)
            for(j=i+1;j<n;j++)
            {
                edge[k].x=i;
                edge[k].y=j;
                dx=ball[i].x-ball[j].x;
                dy=ball[i].y-ball[j].y;
                dz=ball[i].z-ball[j].z;
                ts=pow(dx*dx+dy*dy+dz*dz,0.5);
                if(ts>ball[i].r+ball[j].r)
                    edge[k++].w=ts-ball[i].r-ball[j].r;
                else edge[k++].w=0;
            }
        qsort(edge,k,sizeof(edge[0]),cmp);
        printf("%.3lf\n",Kruskal(k));
    }

}

Acaini.yaoyaozii 2009-07-29 10:05 发表评论

pku 1019 poj

Acaini.yaoyaozii — Tue, 28 Jul 2009 01:43:00 GMT

#include<iostream>
#include<cmath>
using namespace std;
#define M 40000
__int64 sum[M+1];
__int64 len[M+1];
void init()
{
    int i;
    len[0]=0;
    sum[0]=0;
    for(i=1;i<=M;i++)
    {
        len[i]=len[i-1]+(int)log10(double(i))+1;
        sum[i]=sum[i-1]+len[i];
    }

}
int search(__int64 n)
{
    __int64 k,w,li,i,j;
    k=1;
    while(sum[k]<n) k++;//di k zu
    w=n-sum[k-1];//
    i=1;
    while(w-int(log10(double(i))+1)>0)
    {
        w-=int(log10(double(i))+1);
        i++;
    }
    li=(int)log10((double)i)+1;
    for(j=1;j<=li-w;j++)
        i/=10;
    return i%10;
}


int main()
{
    init();
    int t;
    __int64 n;

    scanf("%d",&t);

    for(;t--;)
    {
        scanf("%I64d",&n);
        printf("%d\n",search(n));
    }
    return 0;
}

Acaini.yaoyaozii 2009-07-28 09:43 发表评论

pku poj 3126

Acaini.yaoyaozii — Thu, 23 Jul 2009 08:27:00 GMT

筛个素数表+BFS
貌似数据里没有impossible的情况

#include<iostream>
using namespace std;
#define M 10000
bool prime[M];
bool visit[M];
int s,e;
int step;
void crearprimetable()//筛素数表
{
       int i,j;
       //prime[0]=prime[1]=1;
       for(i=2;i<M;i++)      //判断i是不是素数
           if(prime[i]==0)        //是
               for(j=i+i;j<M;j+=i) //把是i倍数的数j记上1,不是素数.
                   prime[j]=1;

}

void bfs()
{
    int i,j,k,num,t;
    int que[10000];
    int front,rear,trear;
    step=0;
    if(s==e) return ;
    front=rear=0;
    que[rear++]=s;
    while(front<rear)
    {
        step++;
        trear=rear;
        while(front<trear)
        {
            num=que[front++];
            for(i=1;i<10;i++)//换个位
            {
                t=num/10*10+i;
                if(!prime[t] && !visit[t])
                    if(t==e) return ;
                    else  {que[rear++]=t;visit[t]=true;}
            }
            for(i=0;i<10;i++)//换十位
            {
                t=num/100*100+i*10+num%10;
                if(!prime[t] && !visit[t])
                    if(t==e) return ;
                    else  {que[rear++]=t;visit[t]=true;}
            }
            for(i=0;i<10;i++)//换百位
            {
                t=num/1000*1000+i*100+num%100;
                if(!prime[t] && !visit[t])
                    if(t==e) return ;
                    else  {que[rear++]=t;visit[t]=true;}
            }
            for(i=1;i<10;i++)//换百位
            {
                t=i*1000+num%1000;
                if(!prime[t] && !visit[t])
                    if(t==e) return ;
                    else  {que[rear++]=t;visit[t]=true;}
            }
        }
    }

}
int main()
{
    int n;
    crearprimetable();
    scanf("%d",&n);
    while(n--)
    {
        scanf("%d%d",&s,&e);
        memset(visit,0,sizeof(visit));
        bfs();
        printf("%d\n",step);
    }
    return 0;
}

Acaini.yaoyaozii 2009-07-23 16:27 发表评论

pku poj 3087

Acaini.yaoyaozii — Thu, 23 Jul 2009 07:12:00 GMT

//模拟洗牌，不知道为什么分类是bfs

#include<iostream>
using namespace std;
char s1[101];
char s2[101];
char s12[201];//每次洗牌后的串
char s[201]; //目标串
char first[201];//第一次洗牌后的串
int n,len,len2;
int step;
void shuf()
{
    int i,j;
    step=1;
    while(step)
    {
        for(i=0,j=0;i<len;i++)
        {
            s12[j++]=s2[i];
            s12[j++]=s1[i];
        }
        len2=j;
        if(step==1) strcpy(first,s12);
        if(strcmp(s,s12)==0) return ;
        if(step!=1 && strcmp(s12,first)==0) { step=-1;return ;}
        for(j=0;j<len;j++)
            s1[j]=s12[j];
        for(;j<len2;j++)
            s2[j-len]=s12[j];
        step++;
    }
}
int main()
{
    int i,j,k;
    scanf("%d",&n);
    for(k=1;k<=n;k++)
    {
        memset(s,0,sizeof(s));
        memset(first,0,sizeof(first));
        memset(s12,0,sizeof(s12));
        scanf(" %d %s %s %s",&len,&s1,&s2,&s);
        shuf();
        if(step==-1) printf("%d -1\n",k);
        else printf("%d %d\n",k,step);
    }
    return 0;
}

Acaini.yaoyaozii 2009-07-23 15:12 发表评论

pku poj 1426

Acaini.yaoyaozii — Thu, 23 Jul 2009 05:57:00 GMT

//菜鸟做法，BFS,n==99，内存空间无敌的大...

#include<iostream>
using namespace std;
//usig64_18,446,744,073,709,551,615_20wei
__int64 que[800000];
__int64 n;
__int64 t;
__int64 bfs()
{

    int front=0,rear=0;
    que[rear++]=1;
    while(front<rear)
    {
        t=que[front++];
        t*=10;
        if(t%n==0) return t;
        que[rear++]=t;
        t+=1;
        if(t%n==0) return t;
        que[rear++]=t;
    }
}
int main()
{
    while(scanf("%I64d",&n)!=EOF && n)
    {

        printf("%I64d\n",bfs());
    }
    return 0;
}

Acaini.yaoyaozii 2009-07-23 13:57 发表评论

pku poj 2251

Acaini.yaoyaozii — Thu, 23 Jul 2009 03:39:00 GMT

这题想起魔方了,改改A了2225

#include<iostream>
using namespace std;
char dung[35][35][35];
bool visit[35][35][35];
int dx[6]={-1,0,0,0,0,1};
int dy[6]={0,-1,0,0,1,0};
int dz[6]={0,0,-1,1,0,0};
int que[82000];//3*30*30*30,一开始开得太小了2700,27000,30000,40000,50000,640000全TM的WA||RE
int L,R,C;
int sx,sy,sz,ex,ey,ez;
int step;
void bfs()
{
    int i,x,y,z,tx,ty,tz;
    int rear,front,trear;
    bool reach=false;
    memset(que,0,sizeof(que));
    rear=front=0;
    que[rear++]=sx;
    que[rear++]=sy;
    que[rear++]=sz;
    visit[sx][sy][sz]=true;
    while(front<rear )//&& !visit[ex][ey][ez])
    {
        if(visit[ex][ey][ez]) {reach=true;break;}//找到E，停
        trear=rear;
        step++;
        while(front<trear )//&& !visit[ex][ey][ez])
        {
            //if(visit[ex][ey][ez]) {reach=true;break;}
            x=que[front++];
            y=que[front++];
            z=que[front++];
            for(i=0;i<6;i++)
            {
                tx=x+dx[i];
                ty=y+dy[i];
                tz=z+dz[i];
                if(tx>=0 && tx<L && ty>=0 && ty<R && tz>=0 && tz<C &&  !visit[tx][ty][tz] && dung[tx][ty][tz]!='#')
                {
                    visit[tx][ty][tz]=true;
                    que[rear++]=tx;
                    que[rear++]=ty;
                    que[rear++]=tz;
                }
            }//for
        }//while
    }//while
    if(!reach) step=-1;//找不到E,trapped
}
int main()
{
    int i,j,k;
    while(scanf("%d%d%d",&L,&R,&C)!=EOF && L && R && C)
    {
        memset(visit,0,sizeof(visit));
        for(i=0;i<L;i++)
            for(j=0;j<R;j++)
                scanf(" %s",dung[i][j]);
        for(i=0;i<L;i++)
            for(j=0;j<R;j++)
                for(k=0;k<C;k++)
                {
                    if(dung[i][j][k]=='S') {sx=i;sy=j;sz=k;}
                    else if(dung[i][j][k]=='E') {ex=i;ey=j;ez=k;}
                }
        step=0;
        bfs();
        if(step==-1) printf("Trapped!\n");
        else printf("Escaped in %d minute(s).\n",step);

    }
    return 0;
}

Acaini.yaoyaozii 2009-07-23 11:39 发表评论

pku poj 3009

Acaini.yaoyaozii — Wed, 22 Jul 2009 15:46:00 GMT

//很有趣的一题，让我想起口袋怪兽
WA了老半天，原来是两个if()放反了，居然不判RE判WA，o(╯□╰)o

#include<iostream>
using namespace std;
char map[21][21];
int dx[4]={0,-1,0,1};
int dy[4]={-1,0,1,0};
int h,w,sx,sy;
int step;
int res;
void dfs(int x,int y)
{
    int dir,tx,ty;
    if(step>=10) return ;
    for(dir=0;dir<4;dir++)
    {
        tx=x;ty=y;
        while(1)
        {
            tx=tx+dx[dir]; ty=ty+dy[dir];
            if(tx<0 || tx>=h || ty<0 || ty>=w   )break;
            if(map[tx][ty]=='3')
            {
                step++;
                //printf("--%d--\n",step);
                if(res>step) res=step;
                step--;
                return ;
            }
            else
                if(map[tx][ty]=='1')
                {
                    tx=tx-dx[dir];ty=ty-dy[dir];
                    if(tx!=x || ty!=y)
                    {
                        map[tx+dx[dir]][ty+dy[dir]]='0';
                        step++;
                        dfs(tx,ty);
                        step--;
                        map[tx+dx[dir]][ty+dy[dir]]='1';
                    }
                    break;
                }
        }//while
    }//for
}//dfs

int main()
{
    int i,j;

    while(scanf("%d%d",&w,&h)!=EOF && w)
    {
        step=0;
        for(i=0;i<h;i++)
            for(j=0;j<w;j++)
            scanf(" %c",&map[i][j]);
        for(i=0;i<h;i++)
            for(j=0;j<w;j++)
                if(map[i][j]=='2') {sx=i;sy=j;break;}
        res=9999999;
        dfs(sx,sy);
        if(res>10) printf("-1\n");
        else printf("%d\n",res);
    }
    return 0;
}

Acaini.yaoyaozii 2009-07-22 23:46 发表评论

pku 2151

Acaini.yaoyaozii — Sun, 19 Jul 2009 09:13:00 GMT

//解释转的~~~~~~
『题目大意』
一次比赛中，共M道题，T个队，p[i][j]表示队i解出题j的概率；问每队至少解出一题且

冠军队至少解出N道题的概率。

『算法』
设a[i][j][k]表示第i队在前j道题中共解出k道题的概率，易得a[i][j][k]有如下递推
关系(另需考虑边界条件)：

a[i][j][k] = a[i][j-1][k-1] * p[i][j] + a[i][j-1][k] * (1-p[i][j])

设s[i][j]表示a[i][M][0] + a[i][M][1] + ... + a[i][M][j]

问题的解可以转化为：每队均至少做一题的概率（用P1表示）减去每队做题数均在1到N-1

之间的概率（用P2表示）。

P1 = (s[1][M] - s[1][0])*(s[2][M]-s[2][0])*...*(s[T][M]-s[T][0])
P2 = (s[1][N-1] - s[1][0])*(s[2][N-1]-s[2][0])*...*(s[T][N-1]-s[T][0])

『算法复杂度』
O(T*M^2)

『说明』
感谢UESTC的zhucheng在poj的提示！

#include<iostream>
using namespace std;
#define MM 30
#define MT 1000

double p[MT+1][MM+1];
double d[MT+1][MM+1][MM+1];
double MTO[MT+1]; //每队至少做出一题的概率
double LTN[MT+1];//少于N道，亦即1N-1

int main()
{
    int i,j,k;
    int M,T,N;
    double tmp1,tmp2;
    while(scanf("%d%d%d",&M,&T,&N)!=EOF && M)
    {
        memset(MTO,0,sizeof(MTO));
        memset(LTN,0,sizeof(LTN));
        for(i=1;i<=T;i++)
            for(j=1;j<=M;j++)
                scanf("%lf",&p[i][j]);
        for(i=1;i<=T;i++)
        {
            d[i][0][0]=1;
            for(j=1;j<=M;j++)
            {
                d[i][j][0] = d[i][j-1][0]*(1-p[i][j]);
                for(k=1;k<=M;k++)
                    d[i][j][k]=p[i][j]*d[i][j-1][k-1]+(1-p[i][j])*d[i][j-1][k];
            }
        }
        tmp1=tmp2=1.0;
        for(i=1;i<=T;tmp1*=MTO[i],i++)
            for(k=1;k<=M;k++)
                MTO[i]+=d[i][M][k];
        for(i=1;i<=T;tmp2*=LTN[i],i++)
            for(k=1;k<N;k++)
                LTN[i]+=d[i][M][k];
        printf("%.3lf\n",tmp1-tmp2);
    }

    return 0;
}

附上discuss上一组数据：
10 20 10
0.1 0.9 0.8 1 0.9 0.8 1 0.9 0.8 0.2
1 0.9 0.8 1 0.9 0.8 1 0.9 0.8 0.8
0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.7
0.2 0.23 0.56 0.2 0.23 0.56 0.2 0.23 0.56 0.88
0.56 0.2 0.23 0.56 0.88 0.56 0.2 0.23 0.56 0.88
0.37 0.99 0.12 0.82 0.47 0.37 0.99 0.12 0.82 0.47
0.82 0.47 0.37 0.99 0.12 0.82 0.47 0.37 0.99 0.12
0.37 0.99 0.12 0.82 0.472 0.373 0.99 0.12 0.82 0.47
0.472 0.373 0.99 0.12 0.82 0.472 0.373 0.99 0.12 0.82
0.1 0.9 0.8 1 0.9 0.8 1 0.9 0.8 0.2
1 0.9 0.8 1 0.9 0.8 1 0.9 0.8 0.8
0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.7
0.2 0.23 0.56 0.2 0.23 0.56 0.2 0.23 0.56 0.88
0.56 0.2 0.23 0.56 0.88 0.56 0.2 0.23 0.56 0.88
0.37 0.99 0.12 0.82 0.47 0.37 0.99 0.12 0.82 0.47
0.82 0.47 0.37 0.99 0.12 0.82 0.47 0.37 0.99 0.12
0.37 0.99 0.12 0.82 0.472 0.373 0.99 0.12 0.82 0.47
0.472 0.373 0.99 0.12 0.82 0.472 0.373 0.99 0.12 0.82
0.56 0.88 0.56 0.2 0.23 0.373 0.99 0.12 0.82 0.472
0.472 0.373 0.99 0.12 0.82 0.82 0.472 0.373 0.99 0.33

结果：0.740

Acaini.yaoyaozii 2009-07-19 17:13 发表评论

pku 1094

Acaini.yaoyaozii — Mon, 13 Jul 2009 06:00:00 GMT

#include<iostream>
using namespace std;
#define M 26
bool grap[M+1][M+1];
int squ[M+1];
char equ[M*M][4];
int ind[M+1];
int indcop[M+1];
int cnt,n,m;

int topsort()
{
    int i,j,p,num=0,
        flag=0,
        now=0;
    for(i=0;i<n;i++)
        indcop[i]=ind[i];

    for(i=0;i<n;i++)
    {
        num=0;
        for(j=0;j<n;j++)
            if(indcop[j]==0)
            {
                num++;
                p=j;
            }
            if(num==0)          //有回路
                return 0;
            if(num>1)           //不确定
                flag=1;
                //return -1; //不确定的情况下还要判断有没回路，o(╯□╰)o
            indcop[p]=-1;           //置为负值
            squ[now++]=p;
            for(j=0;j<n;j++)     //删除相关的边
                if(grap[p][j])
                    indcop[j]--;
    }
    if(flag)       return -1; //不确定
    return 1;
}

int main()
{
    int i,j,ret;
    char l,r;
    bool done;
    while(scanf("%d%d",&n,&m)==2 && n )
    {
        for(i=0;i<m;i++)
            scanf(" %s",&equ[i]);
        done=false;
        memset(ind,0,sizeof(ind));
        memset(grap,0,sizeof(grap));
        for(i=0;i<m;i++)
        {
            l=equ[i][0];r=equ[i][2];
            grap[l-'A'][r-'A']=true;
            ind[r-'A']++;
            ret=topsort();
            if( ret==0 ) { printf("Inconsistency found after %d relations.\n",i+1); done=true; break;}
            if( ret==1)
            {
                printf("Sorted sequence determined after %d relations: ",i+1);
                for(j=0;j<n;j++)
                    printf("%c",'A'+squ[j]);
                printf(".\n");
                done=true;
                break;
            }
        }
        if(!done) printf("Sorted sequence cannot be determined.\n");
    }
    return 0;
}

Acaini.yaoyaozii 2009-07-13 14:00 发表评论