社区讨论
我爱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 条回复,欢迎继续交流。
正在加载回复...