hdu_4507_吉哥系列故事——恨7不成妻(鬼畜数位DP)
2016-05-29 23:35
411 查看
题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=4507
题意:中文,不解释,注意的是求的是合法数字的平方和,即(a+b+c+……)^2
题解:数位DP,要维护三个值,cnt为合法的个数,sum为这cnt个合法个数的和,sqsum为这cnt个数的平方和dp[i][j][k]为第i位前i位的数字和%7的余数为j,前i位的真实值%7为k的状态,求平方和 要展开 如:(a+b)^2=a^2+2*a*b+b^2
#include <cstdio>
#include <cstring>
#define F(i,a,b) for(int i=a;i<=b;++i)
typedef long long LL;
const LL mod=1e9+7;
struct node{ LL cnt,sum,sqsum;}dp[20][10][10];
int bit[20],len;LL p[20],l,r;
node dfs(int pos,int mo,int now,bool inf){
if(pos==-1)return (node){(mo!=0&&now!=0&&!inf),0,0};
if(!inf && dp[pos][mo][now].cnt!=-1)return dp[pos][mo][now];
int end=inf?bit[pos]:9;
node ans,tmp;ans.cnt=ans.sqsum=ans.sum=0;
F(i,0,end){
if(i==7)continue;
tmp=dfs(pos-1,(mo+i)%7,(now*10+i)%7,inf&&i==end),ans.cnt+=tmp.cnt,ans.cnt%=mod;
ans.sum+=(tmp.sum+((p[pos]*i)%mod)*tmp.cnt%mod)%mod,ans.sum%=mod;
ans.sqsum+=(tmp.sqsum+((2*p[pos]*i)%mod)*tmp.sum)%mod,ans.sqsum%=mod;
ans.sqsum+=((tmp.cnt*p[pos])%mod*p[pos]%mod*i*i%mod),ans.sqsum%=mod;
}
if(!inf)dp[pos][mo][now]=ans;
return ans;
}
LL fuck(LL n){for(len=0;n;n/=10)bit[len++]=n%10;return dfs(len-1,0,0,1).sqsum;}
int main(){
int T;p[0]=1;
F(i,1,19)p[i]=(p[i-1]*10)%mod;
F(i,0,19)F(j,0,9)F(k,0,9)dp[i][j][k].cnt=-1;
scanf("%d",&T);
while(T--){
scanf("%I64d%I64d",&l,&r);
LL ans=fuck(r+1);
ans-=fuck(l);
ans=(ans%mod+mod)%mod;
printf("%I64d\n",ans);
}
return 0;
}
题意:中文,不解释,注意的是求的是合法数字的平方和,即(a+b+c+……)^2
题解:数位DP,要维护三个值,cnt为合法的个数,sum为这cnt个合法个数的和,sqsum为这cnt个数的平方和dp[i][j][k]为第i位前i位的数字和%7的余数为j,前i位的真实值%7为k的状态,求平方和 要展开 如:(a+b)^2=a^2+2*a*b+b^2
#include <cstdio>
#include <cstring>
#define F(i,a,b) for(int i=a;i<=b;++i)
typedef long long LL;
const LL mod=1e9+7;
struct node{ LL cnt,sum,sqsum;}dp[20][10][10];
int bit[20],len;LL p[20],l,r;
node dfs(int pos,int mo,int now,bool inf){
if(pos==-1)return (node){(mo!=0&&now!=0&&!inf),0,0};
if(!inf && dp[pos][mo][now].cnt!=-1)return dp[pos][mo][now];
int end=inf?bit[pos]:9;
node ans,tmp;ans.cnt=ans.sqsum=ans.sum=0;
F(i,0,end){
if(i==7)continue;
tmp=dfs(pos-1,(mo+i)%7,(now*10+i)%7,inf&&i==end),ans.cnt+=tmp.cnt,ans.cnt%=mod;
ans.sum+=(tmp.sum+((p[pos]*i)%mod)*tmp.cnt%mod)%mod,ans.sum%=mod;
ans.sqsum+=(tmp.sqsum+((2*p[pos]*i)%mod)*tmp.sum)%mod,ans.sqsum%=mod;
ans.sqsum+=((tmp.cnt*p[pos])%mod*p[pos]%mod*i*i%mod),ans.sqsum%=mod;
}
if(!inf)dp[pos][mo][now]=ans;
return ans;
}
LL fuck(LL n){for(len=0;n;n/=10)bit[len++]=n%10;return dfs(len-1,0,0,1).sqsum;}
int main(){
int T;p[0]=1;
F(i,1,19)p[i]=(p[i-1]*10)%mod;
F(i,0,19)F(j,0,9)F(k,0,9)dp[i][j][k].cnt=-1;
scanf("%d",&T);
while(T--){
scanf("%I64d%I64d",&l,&r);
LL ans=fuck(r+1);
ans-=fuck(l);
ans=(ans%mod+mod)%mod;
printf("%I64d\n",ans);
}
return 0;
}
相关文章推荐
- 详解Android应用中屏幕尺寸的获取及dp和px值的转换
- 基于Android中dp和px之间进行转换的实现代码
- Android中dip、dp、sp、pt和px的区别详解
- LFC1.0.0 版本发布
- Android dpi,dip,dp的概念以及屏幕适配
- Android px、dp、sp之间相互转换
- HP data protector软件学习1--基本角色与基本工作流程
- HP data protector软件学习2--软件组成与界面介绍
- android中像素单位dp、px、pt、sp的比较
- Android对px和dip进行尺寸转换的方法
- Android根据分辨率进行单位转换-(dp,sp转像素px)
- android 尺寸 dp,sp,px,dip,pt详解
- DP问题各种模型的状态转移方程
- POJ-1695-Magazine Delivery-dp
- nyoj-1216-整理图书-dp
- TYVJ1193 括号序列解题报告
- 对DP的一点感想
- TYVJ上一些DP的解题报告
- soj1005. Roll Playing Games
- 01背包问题