CodeForces - 182D Common Divisors( kmp + 最大公约数 )

题目链接:点击进入
题目

在这里插入图片描述
在这里插入图片描述

题意

若串b可以由x个串a连接组成,那么就认为a是b的除数,例如abab的除数为ab与abab
求两个字符串公共除数的数量

思路

看样例大概可以得知,除数实际上就是最小循环节 t 以及最小循环节的倍数 xt ( 长度上要求 len % xt == 0 ),x 转换一下就是,周期 T = len / t 的所有因子数。( 这个结论是保证在串是完全由最小循环节构成的,否则,串的最小循环节就是它自己 )
先判断两个串的最小循环节是否相等,不等就不存在相同除数。
若是相等,就对两个串的周期T,看有多少个相同因子,而求两个数有多少个相同因子,转换一下就是求两个数最大公约数的因子数。

代码
// Problem: Common Divisors
// Contest: Virtual Judge - CodeForces
// URL: https://vjudge.net/problem/CodeForces-182D
// Memory Limit: 262 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

//#pragma GCC optimize(3)//O3
//#pragma GCC optimize(2)//O2
#include<iostream>
#include<string>
#include<map>
#include<set>
//#include<unordered_map>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstring>
#include<stack>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<fstream>
#define X first
#define Y second
#define base 233 
#define INF 0x3f3f3f3f3f3f3f3f
#define pii pair<int,int>
#define lowbit(x) x & -x
#define inf 0x3f3f3f3f
//#define int long long
//#define double long double
//#define rep(i,x,y) for(register int i = x; i <= y;++i)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const double pai=acos(-1.0);
const int maxn=1e6+10;
const int mod=1e9+7;
const double eps=1e-9;
const int N=5e3+10;
/*--------------------------------------------*/
inline int read()
{
    int k = 0, f = 1 ;
    char c = getchar() ;
    while(!isdigit(c)){if(c == '-') f = -1 ;c = getchar() ;}
    while(isdigit(c)) k = (k << 1) + (k << 3) + c - 48 ,c = getchar() ;
    return k * f ;
}
/*--------------------------------------------*/

string s1,s2;
int next1[maxn],next2[maxn],ans;

void getnext(string s,int len,int Next[])
{
	int k=-1,j=0;
	Next[0]=-1;
	while(j<len)
	{
		if(k==-1||s[j]==s[k])
		{
			k++;
			j++;
			Next[j]=k;
		}
		else k=Next[k];
	}
}

int look(int x)
{
	int sum=0;
	for(int i=1;i*i<=x;i++)
	{
		if(x%i==0)
		{
			sum++;
			if(i*i!=x)
				sum++;
		}
	}
	return sum;
}

int gcd(int a,int b)
{
	if(!b) return a;
	else return gcd(b,a%b);
}

int main() 
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);cout.tie(0);
	cin>>s1;
	cin>>s2;
	int l1=s1.size();
	int l2=s2.size();
	getnext(s1,l1,next1);
	getnext(s2,l2,next2);
	int a1=l1-next1[l1];
	int a2=l2-next2[l2];
	string ss1,ss2;
	if(l1%a1) a1=l1;
	if(l2%a2) a2=l2;
	for(int i=0;i<a1;i++)
		ss1+=s1[i];
	for(int i=0;i<a2;i++)
		ss2+=s2[i];
	if(ss1==ss2)
		ans=look(gcd(l1/a1,l2/a2));
	else ans=0;
	cout<<ans<<endl;
    return 0;
} 
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值