codeforces #274 C. Riding in a Lift dp+前缀和优化
codeforces #274 ?C. Riding in a Lift ? dp+前綴和優(yōu)化
Imagine that you are in a building that has exactly?n?floors. You can move between the floors in a lift. Let's number the floors from bottom to top with integers from?1?to?n. Now you're on the floor number?a. You are very bored, so you want to take the lift. Floor number?b?has a secret lab, the entry is forbidden. However, you already are in the mood and decide to make?k?consecutive trips in the lift.
Let us suppose that at the moment you are on the floor number?x?(initially, you were on floor?a). For another trip between floors you choose some floor with number?y?(y?≠?x) and the lift travels to this floor. As you cannot visit floor?b?with the secret lab, you decided that the distance from the current floor?x?to the chosen?y?must be strictly less than the distance from the current floor?x?to floor?b?with the secret lab. Formally, it means that the following inequation must fulfill:?|x?-?y|?<?|x?-?b|. After the lift successfully transports you to floor?y, you write down number?y?in your notepad.
Your task is to find the number of distinct number sequences that you could have written in the notebook as the result of?k?trips in the lift. As the sought number of trips can be rather large, find the remainder after dividing the number by?1000000007?(109?+?7).
InputThe first line of the input contains four space-separated integers?n,?a,?b,?k?(2?≤?n?≤?5000,?1?≤?k?≤?5000,?1?≤?a,?b?≤?n,?a?≠?b).
OutputPrint a single integer — the remainder after dividing the sought number of sequences by?1000000007?(109?+?7).
Sample test(s) input 5 2 4 1 output 2 input 5 2 4 2 output 2 input 5 3 4 1 output 0 NoteTwo sequences?p1,?p2,?...,?pk?and?q1,?q2,?...,?qk?are?distinct, if there is such integer?j?(1?≤?j?≤?k), that?pj?≠?qj.
Notes to the samples:
dp[i][j]表示第i步到達第j層的種數(shù)。
dp[i][j]=∑ dp[i-1][l] ? ?(abs(l-j)<dp(l-b)) ?即從第l層轉移到第j層。
直接轉移是n^3。由于abs(l-j)<dp(l-b)限制了l的范圍,且范圍是連續(xù)的,因此可以用前綴和維護。
#include<bits/stdc++.h> #define REP(i,a,b) for(int i=a;i<=b;i++) #define MS0(a) memset(a,0,sizeof(a))using namespace std;typedef long long ll; const int INF=(1<<29); const int maxn=5010;const ll p=1000000007; ll n,a,b,k; ll dp[maxn][maxn]; ll sum[maxn];int main() {freopen("in.txt","r",stdin);while(cin>>n>>a>>b>>k){MS0(dp);for(int i=1;i<=n;i++) dp[1][i]=(abs(a-i)<abs(a-b))&&(i!=a);sum[0]=0;for(int i=1;i<=n;i++) sum[i]=(sum[i-1]+dp[1][i])%p;for(int i=2;i<=k;i++){for(int j=1;j<=n;j++){if(j<b){dp[i][j]=(dp[i][j]+(sum[(b+j-1)/2]-sum[0]-(sum[j]-sum[j-1])+3*p)%p)%p;}else if(j>b){dp[i][j]=(dp[i][j]+(sum[n]-sum[(b+j)/2]-(sum[j]-sum[j-1])+3*p)%p)%p;}}for(int j=1;j<=n;j++){sum[j]=(sum[j-1]+dp[i][j]%p)%p;}}ll ans=0;for(int i=1;i<=n;i++){ans=(ans%p+dp[k][i]%p)%p;}cout<<ans<<endl;}return 0; } View Code?
轉載于:https://www.cnblogs.com/--560/p/4755401.html
總結
以上是生活随笔為你收集整理的codeforces #274 C. Riding in a Lift dp+前缀和优化的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: oracle学习笔记一:用户管理(2)创
- 下一篇: 转【FullPage.js 应用参数参考