社区讨论

我爱TLE

P2152[SDOI2009] SuperGCD参与者 3已保存回复 6

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
6 条
当前快照
1 份
快照标识符
@mk29dsrt
此快照首次捕获于
2026/01/06 15:19
上个月
此快照最后确认于
2026/01/09 20:20
上个月
查看原帖
CPP
#include<bits/stdc++.h>
#define int short
using namespace std;
int a[114514],b[114514],c[114514];
const int x=__gcd(1,2);
/*
相等情况:如果 a = b,则 gcd(a, b) = a。
任一为零:如果 a = 0 或 b = 0,则 gcd(a, b) = max(a, b)。
均为偶数:如果 a 和 b 都是偶数,则 gcd(a, b) = 2 * gcd(a/2, b/2)。
一偶一奇:如果 a 是偶数,b 是奇数,则 gcd(a, b) = gcd(a/2, b)。
均为奇数:如果 a 和 b 都是奇数,则 gcd(a, b) = gcd(a-b, b)。 
这种算法通过位运算和递归的方式,能够显著提高计算效率。
*/
inline bool compare()
{
	if(a[0]==b[0])
	{
		for(int i=1;i<=a[0];i++)
		{
			if(a[i]!=b[i]) return 0;
		}
		return 1;
	}
	else return 0;
}
inline void dbt(int s[])//Divide by two
{
	int tw=0;
	for(int i=s[0];i>=1;i--)
	{
		int now=tw*10+s[i];
		s[i]=now>>1;
		tw=now&1;
	}
	while(s[0]>0&&s[s[0]]==0)
	{
		s[0]--;
	}
}
inline void jian()
{
	int wei=0;
	memset(c,0,sizeof(c));
	c[0]=a[0]; 
	for(int i=1;i<=c[0];i++)
	{
		c[i]+=a[i]-(i<=b[0]?b[i]:0);
		if(c[i]<0) 
		{
			c[i]+=10;
			c[i+1]--;
		}
	}
	while(c[0]>0&&c[c[0]]==0)
	{
		c[0]--;
	}
	memset(a,0,sizeof(a));
	a[0]=c[0];
	for(int i=1;i<=a[0];i++)
	{
		a[i]=c[i];
	}
}
inline bool a_ge_b()
{
	if(a[0]>b[0]) return 1;
	if(a[0]<b[0]) return 0;
	for(int i=a[0];i>=1;i--)
	{
		if(a[i]>b[i]) return 1;
		if(a[i]<b[i]) return 0;
	}
	return 1;
}
inline void swap_ab()
{
	int t[114514];
	memcpy(t,a,sizeof(a));
	memcpy(a,b,sizeof(b));
	memcpy(b,t,sizeof(b));
}
inline void mul2k(int s[],int k)
{
	while(k--)
	{
		int o=0;
		for(int i=1;i<=s[0];i++)
		{
			int now=s[i]*2+o;
			s[i]=now%10;
			o=now/10;
		}
		if(o>0)
		{
			s[++s[0]]=o;
		}
	}
}
inline string exchange(int ss[])
{
	string s="";
	if(ss[0]==0) return "0";
	for(int i=ss[0];i>=1;i--)
	{
		s+=(ss[i]+'0');
	}
	return s;
}
inline bool ou(int s[])
{
	if(s[0]==0) return 1;
	return !(s[1]&1);
}
inline string supergcd()
{
	int k=0; 
	while(1)
	{
		if(compare()) break;
		if(a[0]==0) {swap_ab();break;}
		if(b[0]==0) break;
		if(ou(a)&&ou(b))
		{
			k++;
			dbt(a);
			dbt(b);
		}
		else if(ou(a)&&!ou(b))
		{
			dbt(a);
		}
		else if(!ou(a)&&ou(b))
		{
			dbt(b);
		}
		else
		{
			if(!a_ge_b()) swap_ab();
			jian();
		}
	}
	mul2k(a,k);
	return exchange(a);
}
signed main()
{
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	string aa,bb;
	cin>>aa>>bb;
	a[0]=aa.size();
	for(int i=1;i<=a[0];i++)
	{
		a[i]=aa[aa.size()-i]-'0';
	}
	b[0]=bb.size();
	for(int i=1;i<=b[0];i++)
	{
		b[i]=bb[bb.size()-i]-'0';
	}
	cout<<supergcd();
	return 0;
}

屎山代码求卡常

回复

6 条回复,欢迎继续交流。

正在加载回复...