zoj 3228 Searching the String【ac自动机】

2896454 2012-05-18 16:35:19 Accepted 3228 C++ 2100 72084 zz_1215
2896441 2012-05-18 16:30:35 Compilation Error 3228 C++ 0 0 zz_1215
2896425 2012-05-18 16:26:32 Compilation Error 3228 C++ 0 0 zz_1215
2896395 2012-05-18 16:11:15 Wrong Answer 3228 C++ 1980 78916 zz_1215
2896390 2012-05-18 16:05:16 Time Limit Exceeded 3228 C++ 7001 75364 zz_1215
2896381 2012-05-18 15:57:51 Time Limit Exceeded 3228 C++ 7001 0 zz_1215      

               一开始用了个vector保存id,以为能ac,没想到数据这么变态,肯定有好多重复的,所以超时了,注释掉find()函数之后再交,1980ms,说明其他的操作不会超时,加一个link数组解决这个问题,在建立字典树的时候顺便把相同的连接起来,没想到居然会CE,我了个擦的我可不是拿记事本写的啊!!在codeblock上和devc上都编译通过了啊,再交一次,一看返回,居然有个函数的名字也是link。。。我就晕了,把link改成linking。。ac
#include<iostream>
#include<vector>
#include<algorithm>
#include<cstdio>
#include<queue>
#include<stack>
#include<string>
#include<map>
#include<set>
#include<cmath>
#include<cassert>
#include<cstring>
#include<iomanip>
using namespace std;

#ifdef _WIN32
#define i64 __int64
#define out64 "%I64d\n"
#define in64 "%I64d"
#else
#define i64 long long
#define out64 "%lld\n"
#define in64 "%lld"
#endif

#define FOR(i,a,b)      for( int i = (a) ; i <= (b) ; i ++)
#define FF(i,a)         for( int i = 0 ; i < (a) ; i ++)
#define FFD(i,a)        for( int i = (a)-1 ; i >= 0 ; i --)
#define S64(a)          scanf(in64,&a)
#define SS(a)           scanf("%d",&a)
#define LL(a)           ((a)<<1)
#define RR(a)           (((a)<<1)+1)
#define pb              push_back
#define CL(Q)           while(!Q.empty())Q.pop()
#define MM(name,what)   memset(name,what,sizeof(name))
#define read            freopen("in.txt","r",stdin)
#define write           freopen("out.txt","w",stdout)

const int inf = 0x3f3f3f3f;
const i64 inf64 = 0x3f3f3f3f3f3f3f3fLL;
const double oo = 10e9;
const double eps = 1e-10;
const double pi = acos(-1.0);
const int maxz = 600011;
const int maxn = 100011;
const int head = 0;

struct zz
{
    int to[26];
    int fail;
    int id;
}zx[maxz];

char sc[maxn];
string rs,s[maxn];
int flag[maxn];
int n,use;
int at[maxn];
int have[maxn][3];
int linking[maxn];
bool ok[maxn];

inline int get()
{
    use++;
    zx[use].fail = 0;
    zx[use].id = 0;
    MM(zx[use].to,-1);
    return use;
}

void insert(int id,string & p)
{
    int now = head;
    int c;
    for(int i=0;i<p.length();i++)
    {
        c = p[i]-'a';
        if(zx[now].to[c]!=-1)
        {
            now = zx[now].to[c];
        }
        else
        {
            zx[now].to[c] = get();
            now =  zx[now].to[c];
        }
    }
    if(zx[now].id)
    {
        ok[id] = true;
        linking[id] = zx[now].id;
    }
    else
    {
        zx[now].id = id;
    }
    return ;
}

void ac()
{
    queue<int>q;
    CL(q);
    q.push(head);
    int now,to,temp;
    while(!q.empty())
    {
        now = q.front();
        q.pop();
        for(int i=0;i<26;i++)
        {
            if(zx[now].to[i]!=-1)
            {
                to = zx[now].to[i];
                q.push(to);
                temp = now;
                while(temp)
                {
                    temp = zx[temp].fail;
                    if(zx[temp].to[i]!=-1)
                    {
                        temp = zx[temp].to[i];
                        break;
                    }
                }
                zx[to].fail = temp;
            }
        }
    }
    return ;
}

void find(int pos,int now)
{
    if(now)
    {
        int len,id = zx[now].id;
        if(id)
        {
            have[id][0]++;
            len = s[id].length();
            if( pos- len + 1 > at[id] + len - 1)
            {
                have[id][1]++;
                at[id] = pos - s[id].length() + 1;
            }
        }
        find(pos,zx[now].fail);
    }
    return ;
}

void ak(string & p)
{
    int now = head;
    int temp;
    int c;
    for(int i=0;i<p.length();i++)
    {
        c = p[i] - 'a';
        if(zx[now].to[c] != -1)
        {
            now = zx[now].to[c];
        }
        else
        {
            temp = now;
            while(now)
            {
                now = zx[now].fail;
                if(zx[now].to[c]!= -1)
                {
                    now = zx[now].to[c];
                    break;
                }
            }
            zx[temp].to[c] = now;
        }
        find(i,now);
    }
    return ;
}

void start()
{
    ac();

    for(int i=0;i<=n;i++)
    {
        at[i] = -inf;
        have[i][0] = have[i][1] = 0;
    }

    ak(rs);

    return ;
}

int main()
{
    int tt = 1;
    char sd[11];
    while(scanf("%s",sc)!=EOF)
    {
        MM(ok,false);
        rs = sc;
        SS(n);
        use = -1;
        get();
        for(int i=1;i<=n;i++)
        {
            SS(flag[i]);
            scanf("%s",sd);
            s[i] = sd;
            insert(i,s[i]);
        }
        start();
        printf("Case %d\n",tt++);
        for(int i=1;i<=n;i++)
        {
            if(!ok[i])
            {
                printf("%d\n",have[i][flag[i]]);
            }
            else
            {
                printf("%d\n",have[linking[i]][flag[i]]);
            }
        }
        printf("\n");
    }
    return 0;
}


以下是ZOJ1626的C++ AC代码,使用了旋转卡壳算法: ```c++ #include <iostream> #include <cstdio> #include <algorithm> #include <cmath> #include <cstring> #define MAXN 100010 #define eps 1e-8 #define INF 1e20 using namespace std; struct point { double x,y; friend point operator -(point a,point b) { point res; res.x=a.x-b.x; res.y=a.y-b.y; return res; } friend bool operator <(point a,point b) { if(fabs(a.x-b.x)<eps) return a.y<b.y; return a.x<b.x; } friend double operator *(point a,point b) { return a.x*b.y-a.y*b.x; } friend double dis(point a,point b) { return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)); } }a[MAXN],b[MAXN],st[MAXN]; int n; double ans=INF; int cmp(point a,point b) { double tmp=(a-b)*(a[1]-b); if(fabs(tmp)<eps) return dis(a,a[1])-dis(b,a[1])<0; return tmp>0; } int main() { while(~scanf("%d",&n) && n) { for(int i=1;i<=n;i++) scanf("%lf%lf",&a[i].x,&a[i].y); sort(a+1,a+n+1); int tot=0; for(int i=1;i<=n;i++) { while(tot>=2 && (st[tot]-st[tot-1])*(a[i]-st[tot])<0) tot--; st[++tot]=a[i]; } int k=tot; for(int i=n-1;i>=1;i--) { while(tot>k && (st[tot]-st[tot-1])*(a[i]-st[tot])<0) tot--; st[++tot]=a[i]; } tot--; for(int i=1;i<=tot;i++) b[i]=st[i]; int tmp=1; for(int i=2;i<=tot;i++) if(b[i].y<b[tmp].y) tmp=i; swap(b[1],b[tmp]); sort(b+2,b+tot+1,cmp); st[1]=b[1]; st[2]=b[2]; k=2; for(int i=3;i<=tot;i++) { while(k>1 && (st[k]-st[k-1])*(b[i]-st[k])<=0) k--; st[++k]=b[i]; } double ans=0; if(k==2) ans=dis(st[1],st[2]); else { st[k+1]=st[1]; for(int i=1;i<=k;i++) for(int j=1;j<=k;j++) ans=max(ans,dis(st[i],st[j])); } printf("%.2lf\n",ans/2); } return 0; } ``` 其中,结构体 `point` 表示二维平面上的一个点,包含了点的坐标和一些基本操作。函数 `cmp` 是旋转卡壳算法中的比较函数,按照点到起点的极角从小到大排序。在主函数中,先使用 Graham 扫描法求出点集的凸包,然后按照旋转卡壳的步骤,求出凸包上的最远点对距离作为最小直径。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值