https://codeforces.com/problemset/problem/798/C
题意:给定一个数列,可以进行以下的操作,选相邻两个数进行x-y,x+y的操作。问最少多少次使得整个序列gcd>1;
思路:
模拟可以发现最多两次就可以出gcd==2的两个数。
而且如果原来就gcd!=1,是不会通过1次变成gcd>1的;
假设t|x,t|y--->t|ax+by;因此如果1次后gcd可以>1,原来就可以>1。因此不会发生这样事。
我们继续探究为什么两次就能gcd==2
两个数A.B.无非四种情况:
奇数,奇数--------------->操作后变成 偶数,偶数 +1
奇数,偶数--------------->操作后变成 奇数,奇数 ----偶偶 +2
偶数,奇数--------------->操作后变成 奇数,奇数 ---偶偶 +2
偶数,偶数--------------->操作后变成 偶数,偶数--- +0
所以,对于三个数以及以上的,考虑当前是奇数变偶还是不变等。进行dp.
dp[i][0/1]:前i个数且当前数字修改为偶数/奇数时,前i个数字的gcd>1的最小代价
0 0--->0 (不会有偶数变成奇数的状态)
0 1--->2
1 0--->2
1 1--->1
#include<iostream>
#include<vector>
#include<queue>
#include<cstring>
#include<cmath>
#include<map>
#include<set>
#include<cstdio>
#include<algorithm>
#define debug(a) cout<<#a<<"="<<a<<endl;
using namespace std;
const int maxn=1e5+1000;
typedef long long LL;
const LL inf=0x3f3f3f3f;
inline LL read(){LL x=0,f=1;char ch=getchar(); while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}while (isdigit(ch)){x=x*10+ch-48;ch=getchar();}
return x*f;}
LL a[maxn];
LL dp[maxn][2];
///dp[i][0/1]:前i个数且当前数字修改为偶数/奇数时,前i个数字的gcd>1的最小代价
///0 0--->0 (不会有偶数变成奇数的状态)
///0 1--->2
///1 0--->2
///1 1--->1
int main(void)
{
cin.tie(0);std::ios::sync_with_stdio(false);
LL n;cin>>n;
LL gcd=0;
for(LL i=1;i<=n;i++){
cin>>a[i];
gcd=__gcd(gcd,a[i]);
}
if(gcd>1){
cout<<"YES"<<"\n";
cout<<"0"<<"\n";
}
else if(gcd==1){
dp[1][!(a[1]%2)]=inf;///第一个数偶数,奇数就无法变;第一个数是奇数,偶数就不会变
for(LL i=2;i<=n;i++){
if(a[i]%2==0){
dp[i][0]=min(dp[i-1][0],dp[i-1][1]+2);
dp[i][1]=inf;///当前偶数不会变成奇数的
}
else if(a[i]&1){
///把当前的奇数改成偶数
///和前面的偶数改,和前面变成的奇数改
dp[i][0]=min(dp[i-1][0]+2,dp[i-1][1]+1);
dp[i][1]=min(inf,dp[i-1][0]);///当前仍然奇数不动,继承前面改好的偶数状态
}
}
///总会yes的
if(dp[n][0]!=inf){
cout<<"YES"<<"\n";
cout<<dp[n][0]<<"\n";
}
else cout<<"NO"<<"\n";
}
return 0;
}
当然了,dp可能略有麻烦。我们发现可以贪心,奇数奇数最优+1,优先把奇数一块都变了。然后再把剩下的变了。不会更差
#include<iostream>
#include<vector>
#include<queue>
#include<cstring>
#include<cmath>
#include<map>
#include<set>
#include<cstdio>
#include<algorithm>
#define debug(a) cout<<#a<<"="<<a<<endl;
using namespace std;
const int maxn=1e5+1000;
typedef long long LL;
inline LL read(){LL x=0,f=1;char ch=getchar(); while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}while (isdigit(ch)){x=x*10+ch-48;ch=getchar();}
return x*f;}
LL a[maxn];
int main(void)
{
cin.tie(0);std::ios::sync_with_stdio(false);
LL n;cin>>n;
LL gcd=0;
for(LL i=1;i<=n;i++){
cin>>a[i];gcd=__gcd(gcd,a[i]);
}
if(gcd>1){
cout<<"YES"<<"\n";
cout<<"0"<<"\n";
}
else if(gcd==1){
for(LL i=1;i<=n;i++){
if(a[i]&1) a[i]=1;
else if(a[i]%2==0) a[i]=0;
}
LL ans=0;
for(LL i=1;i<n;i++){
if(a[i]==1&&a[i+1]==1) ans+=1,a[i]=0,a[i+1]=0;
}
for(LL i=1;i<n;i++){
if(a[i]==1&&a[i+1]==0) a[i]=0,a[i+1]=0,ans+=2;
else if(a[i]==0&&a[i+1]==1) a[i]=0,a[i+1]=0,ans+=2;
}
cout<<"YES"<<"\n";
cout<<ans<<"\n";
}
return 0;
}