1. 程式人生 > >JZOJ5864. 【NOIP2018模擬9.11】很多序列

JZOJ5864. 【NOIP2018模擬9.11】很多序列

這裡寫圖片描述
這裡寫圖片描述

題解

觀察資料範圍,就發現x1很小,而其他很大,於是就可以根據這裡來想。
將其他數換成ax1+b的形式。
可以知道只有出現了連續x1個數在序列裡面,
那麼最大的數就是在這連續x1的數前面一個數。
fi表示,用x2,x3,...,xn來組成的數,在%x1之後為i,最小在上面時候可以出現在數列中。
這樣fix1+i就是%x1=i第一次出現的時候。
連邊就是t連向t+xi%x1,邊權為xi/x1。如果這個數大於x1,就取模,然後邊權加一。
答案就是(m

axfi1)x1+i

code

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string.h>
#include <cmath>
#include <math.h>
#define ll long long
#define P putchar
#define G getchar
using namespace std;
char ch;
void read(ll &n)
{
    n=0
; ch=G(); while((ch<'0' || ch>'9') && ch!='-')ch=G(); ll w=1; if(ch=='-')w=-1,ch=G(); while('0'<=ch && ch<='9')n=(n<<3)+(n<<1)+ch-'0',ch=G(); n*=w; } ll n,a[13],ans,f[1000003],v[13],w; int q[2000003],head,tail,b[13],x,t; bool bz[1000003]; int main() { freopen("sequence.in"
,"r",stdin); freopen("sequence.out","w",stdout); read(n); for(int i=1;i<=n;i++)read(a[i]),b[i]=a[i]%a[1],v[i]=a[i]/a[1]; if(n==2) { printf("%lld",(a[1]-1)*(a[2]-1)-1); return 0; } memset(f,127,sizeof(f)); memset(bz,1,sizeof(bz)); for(head=q[tail=1]=f[0]=bz[0]=0;head<tail;) { x=q[++head]; for(int i=2;i<=n;i++) { t=x+b[i];w=f[x]+v[i]; if(t>=a[1])w++,t=t-a[1]; if(f[t]>w) { f[t]=w; if(bz[t])bz[q[++tail]=t]=0; } } bz[x]=1; } for(int i=0;i<a[1];i++) if(f[i]>=ans)ans=f[i],x=i; printf("%lld\n",(ans-1)*a[1]+x); return 0; }