C. Mike and gcd problem(思维+数论+dp/贪心好题)

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;
}