CodeForces 623B Array GCD （数论+DP）*好题

题目链接：http://codeforces.com/problemset/problem/623/B

题目大意：

给定一个整数序列，
给定两种操作，一种是删除一整段区间，
一种是对一个指定的数修改成x+1,或者x-1,
其中要求是操作一不能删除全部的序列且只能删除一段，
操作二对每个数只能进行一次（不是操作只能进行一次，是对每个数~）

题目分析： 

这道题我是看了别人的思路才码出来的，
实现层面来讲这道题不是很难，重点在想清楚。
首先根据题意，要求是不能删除全部的序列，
那么简单的容斥下不难发现，第一个数和最后一个数至少有一个是存在的，
那么我们不难暴力枚举其整个序列的公共因子gcd。
假设我们已经确定了gcd为fac，
对于每个数我们都指定其执行第二次操作时需要的代价，
如果本身就可以整除则代价为0，
如果需要浮动1才可以整除则代价为b，
如果以上两种条件都行不通则直接赋值无穷大，意味着这一位操作不可取。

对于删除整个区间的操作，我们设定三种状态，dp(i,0)表示没有开始进行删除操作，
dp(i,1)表示正在删除中，dp(i,2)表示删除已经结束时候的状态。
递推方程不难想到，详见代码。
时间复杂度：O(n*C)。

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
#define pii pair<int,int>
#define fi first
#define se second
#define mk(x,y) make_pair(x,y)
const int mod=1e9+7;
const int maxn=1e6+5;
const int ub=1e6;
const double inf=1e-4;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
题目大意：
给定一个整数序列，
给定两种操作，一种是删除一整段区间，
一种是对一个指定的数修改成x+1,或者x-1,
其中要求是操作一不能删除全部的序列且只能删除一段，
操作二对每个数只能进行一次（不是操作只能进行一次，是对每个数~）

题目分析：
这道题我是看了别人的思路才码出来的，
实现层面来讲这道题不是很难，重点在想清楚。
首先根据题意，要求是不能删除全部的序列，
那么简单的容斥下不难发现，第一个数和最后一个数至少有一个是存在的，
那么我们不难暴力枚举其整个序列的公共因子gcd。
假设我们已经确定了gcd为fac，
对于每个数我们都指定其执行第二次操作时需要的代价，
如果本身就可以整除则代价为0，
如果需要浮动1才可以整除则代价为b，
如果以上两种条件都行不通则直接赋值无穷大，意味着这一位操作不可取。

对于删除整个区间的操作，我们设定三种状态，dp(i,0)表示没有开始进行删除操作，
dp(i,1)表示正在删除中，dp(i,2)表示删除已经结束时候的状态。
递推方程不难想到，详见代码。
时间复杂度：O(n*C)。
*/
int n,a,b,x;
int seq[maxn];
ll dp[maxn][3],val[maxn],INF;
ll solve(int fac){
    INF=1e18;
    ll ans=INF;
    rep(i,1,n+1)
        if(seq[i]%fac==0) val[i]=0;
        else if((seq[i]+1)%fac==0 || (seq[i]-1)%fac==0) val[i]=b;
        else val[i]=INF;///
    ///rep(i,0,maxn) rep(j,0,3) dp[i][j]=INF;
    dp[0][0]=dp[0][1]=dp[0][2]=0;///初始化
    rep(i,1,n+1){
        dp[i][0]=dp[i-1][0]+val[i];
        dp[i][0]=min(dp[i][0],INF);
        dp[i][1]=min(dp[i-1][1],dp[i-1][0])+a;
        dp[i][1]=min(dp[i][1],INF);
        dp[i][2]=min(dp[i-1][2],dp[i-1][1])+val[i];
        dp[i][2]=min(dp[i][2],INF);
    }
    return min(min(dp[n][0],dp[n][1]),dp[n][2]);
}
///四种情况，从1到n-1,从2到n，b操作用完，和b操作没用。
int f[maxn/10],cnt=0;
void divid(int x){
    for(int i=2;i*i<=x;i++) if(x%i==0){
        while(x%i==0) x/=i;
        f[cnt++]=i;
    }
    if(x>1) f[cnt++]=x;
}
int main(){
    cin>>n>>a>>b;
    rep(i,1,n+1) cin>>seq[i];
    divid(seq[1]+1),divid(seq[1]-1),divid(seq[1]);
    divid(seq[n]+1),divid(seq[n]-1),divid(seq[n]);
    sort(f,f+cnt);cnt=unique(f,f+cnt)-f;
    ll ans=1e18;
    rep(i,0,cnt) ans=min(ans,solve(f[i]));
    cout<<ans<<endl;
    return 0;
}