P4173 残缺的字符串(FFT)

本文详细解析了Luogu4173题目,通过定义匹配函数和完全匹配函数,利用快速傅里叶变换(FFT)进行字符串匹配,解决了复杂的大规模字符串比较问题。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

[Luogu4173]

题解

\(1.\)定义匹配函数

\(2.\)定义完全匹配函数

\(3.\)快速计算每一位的完全匹配函数值

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#define debug(...) fprintf(stderr,__VA_ARGS__)
#define Debug(x) cout<<#x<<"="<<x<<endl
using namespace std;
typedef long long LL;
const int INF=1e9+7;
inline LL read(){
    register LL x=0,f=1;register char c=getchar();
    while(c<48||c>57){if(c=='-')f=-1;c=getchar();}
    while(c>=48&&c<=57)x=(x<<3)+(x<<1)+(c&15),c=getchar();
    return f*x;
}
inline char readc(){
    register char c=getchar();
    while(c==' '||c=='\n'||c=='\t') c=getchar();
    return c;
}

const int MAXN=2e6+5;
const double Pi=acos(-1);

struct cmpx{
    double x,y;
    inline cmpx(){}
    inline cmpx(double _x,double _y){x=_x,y=_y;}
    inline friend cmpx operator + (cmpx a,cmpx b){return cmpx(a.x+b.x,a.y+b.y);}
    inline friend cmpx operator - (cmpx a,cmpx b){return cmpx(a.x-b.x,a.y-b.y);}
    inline friend cmpx operator * (cmpx a,cmpx b){return cmpx(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);}
    inline friend cmpx operator * (cmpx a,double b){return cmpx(a.x*b,a.y*b);}
}A[MAXN],B[MAXN],C[MAXN];

char s[MAXN],t[MAXN];
int ans[MAXN],a[MAXN],b[MAXN];
int n,m;

namespace F_F_T{
    int rev[MAXN],limit,l;
    inline void init(int n){
        for(limit=1,l=0;limit<=n;limit<<=1) l++;
        for(int i=0;i<limit;i++)
            rev[i]=(rev[i>>1]>>1)|((i&1)<<(l-1));
    }
    inline void FFT(cmpx *A,int type){
        for(int i=0;i<limit;i++)
            if(i<rev[i]) swap(A[i],A[rev[i]]);
        for(int len=1;len<limit;len<<=1){
            cmpx Wn=(cmpx){cos(Pi/len),type*sin(Pi/len)};
            for(int i=0;i<limit;i+=(len<<1)){
                cmpx w=(cmpx){1,0};
                for(int j=0;j<len;j++,w=w*Wn){
                    cmpx x=A[i+j],y=w*A[i+len+j];
                    A[i+j]=x+y;
                    A[i+len+j]=x-y;
                }
            }
        }
        if(type==-1){
            for(int i=0;i<limit;i++) A[i].x/=limit;
        }
    }
}using namespace F_F_T;

int main(){
    n=read(),m=read();
    init(n+m-2);
    for(int i=n-1;i>=0;i--){
        char c=readc();
        a[i]=(c!='*')?c-'a'+1:0;
    }
    for(int i=0;i<m;i++){
        char c=readc();
        b[i]=(c!='*')?c-'a'+1:0;
    }

    for(int i=0;i<n;i++) A[i]=cmpx(a[i]*a[i]*a[i],0);
    for(int i=0;i<m;i++) B[i]=cmpx(b[i],0);
    FFT(A,1);FFT(B,1);
    for(int i=0;i<limit;i++) C[i]=A[i]*B[i];

    for(int i=0;i<limit;i++) A[i]=B[i]=cmpx(0,0);
    for(int i=0;i<n;i++) A[i]=cmpx(a[i],0);
    for(int i=0;i<m;i++) B[i]=cmpx(b[i]*b[i]*b[i],0);
    FFT(A,1);FFT(B,1);
    for(int i=0;i<limit;i++) C[i]=C[i]+A[i]*B[i];

    for(int i=0;i<limit;i++) A[i]=B[i]=cmpx(0,0);
    for(int i=0;i<n;i++) A[i]=cmpx(a[i]*a[i],0);
    for(int i=0;i<m;i++) B[i]=cmpx(b[i]*b[i],0);
    FFT(A,1);FFT(B,1);
    for(int i=0;i<limit;i++) C[i]=C[i]-A[i]*B[i]*2.0;

    FFT(C,-1);

    for(int i=n-1;i<m;i++){
        if(fabs(C[i].x)<0.5) ans[++ans[0]]=i-n+2;
        //printf("%.3lf\n",C[i].x);
    }
    printf("%d\n",ans[0]);
    for(int i=1;i<=ans[0];i++)
        printf("%d ",ans[i]);
}

转载于:https://www.cnblogs.com/lizehon/p/10504826.html

内容概要:本文探讨了在MATLAB/SimuLink环境中进行三相STATCOM(静态同步补偿器)无功补偿的技术方法及其仿真过程。首先介绍了STATCOM作为无功功率补偿装置的工作原理,即通过调节交流电压的幅值和相位来实现对无功功率的有效管理。接着详细描述了在MATLAB/SimuLink平台下构建三相STATCOM仿真模型的具体步骤,包括创建新模型、添加电源和负载、搭建主电路、加入控制模块以及完成整个电路的连接。然后阐述了如何通过对STATCOM输出电压和电流的精确调控达到无功补偿的目的,并展示了具体的仿真结果分析方法,如读取仿真数据、提取关键参数、绘制无功功率变化曲线等。最后指出,这种技术可以显著提升电力系统的稳定性与电能质量,展望了STATCOM在未来的发展潜力。 适合人群:电气工程专业学生、从事电力系统相关工作的技术人员、希望深入了解无功补偿技术的研究人员。 使用场景及目标:适用于想要掌握MATLAB/SimuLink软件操作技能的人群,特别是那些专注于电力电子领域的从业者;旨在帮助他们学会建立复杂的电力系统仿真模型,以便更好地理解STATCOM的工作机制,进而优化实际项目中的无功补偿方案。 其他说明:文中提供的实例代码可以帮助读者直观地了解如何从零开始构建一个完整的三相STATCOM仿真环境,并通过图形化的方式展示无功补偿的效果,便于进一步的学习与研究。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值