传送门

三个限制都可以数位 $dp$ , $dfs$ 是维护当前位,之前各位总和模 $7$ 意义下的值,之前填的数模 $7$ 意义下的值,是否贴着限制

主要现在求的是各个合法数的平方的和,比较恶心

开一个结构体维护一下三个东西,合法数的个数,合法数的和,合法数的平方和

前两个好维护,平方和发现其实可以拆开,比如下一层 $dfs$ 传上来的和为 $x$ ,当前位填的数为 $k$ ,考虑当前所有合法数的平方和看成每一个合法数的平方的和

设下一层某一个合法数传上来的值为 $y$,那么当前层的值即为

$10^p \cdot k + y$ ,平方后即为 $10^{2p} \cdot k^2 + y^2 + 2ky \cdot 10^p$

发现所有的 $y^2$ 加起来以后其实就是下一层的合法数的平方和

对于 $2ky \cdot 10^p$ ,把 $2k \cdot 10^p$ 提出来,那么也可以把所有 $y$ 加起来即为下一层的合法数的和

对于 $10^{2p} \cdot k^2$ ,因为加了下一层合法数的个数次,那么贡献即为下一层合法数数量乘 $10^{2p} \cdot k^2$

所以我们完全可以用下一层的数据维护出当前位的数据

具体实现可以看代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
using namespace std;
typedef long long ll;
inline ll read()
{
ll x=,f=; char ch=getchar();
while(ch<''||ch>'') { if(ch=='-') f=-; ch=getchar(); }
while(ch>=''&&ch<='') { x=(x<<)+(x<<)+(ch^); ch=getchar(); }
return x*f;
}
const int mo=1e9+;
inline int fk(int x) { return x>=mo ? x-mo : x; }
ll T,L,R;
struct dat {
ll tot,sum,sum2;
dat (ll _tot=,ll _sum=,ll _sum2=) { tot=_tot,sum=_sum,sum2=_sum2; }
inline dat operator + (const dat &tmp) const {
return dat(fk(tot+tmp.tot),fk(sum+tmp.sum),fk(sum2+tmp.sum2));
}
inline dat operator * (const ll &tmp) const {
return dat(tot, fk(tot*tmp%mo + sum) , fk( fk((tot*tmp%mo)*tmp%mo + (2ll*tmp%mo)*sum%mo) + sum2 ) );
}
}f[][][][];
ll Power[];
vector <int> V;
dat dfs(int p,int sum,int now,bool lim)
{
if(!p) return (sum&&now) ? dat(,,) : dat(,,);
dat &t=f[p][sum][now][lim];
if(t.tot!=-) return t;
t.tot=;
for(int k=;k<=;k++)
{
if(lim&&k>V[p-]) break;
if(k==) continue;
dat add=dfs(p-,(sum+k)%,(now*+k)%,lim&(k==V[p-]));
t=t + ( add * (Power[p-]*k%mo) );
}
// cout<<p<<" "<<sum<<" "<<now<<" "<<lim<<" "<<t.tot<<" "<<t.sum<<" "<<t.sum2<<endl;
return t;
}
inline void Clr()
{
for(int i=;i<=;i++)
for(int j=;j<;j++)
for(int k=;k<;k++)
f[i][j][k][]=f[i][j][k][]=dat(-,,);
}
int main()
{
T=read();
Power[]=; for(int i=;i<=;i++) Power[i]=Power[i-]*%mo;
while(T--)
{
L=read()-,R=read();
ll now=R; V.clear(); Clr();
while(now) V.push_back(now%),now/=;
ll ans=dfs(V.size(),,,).sum2;
now=L; V.clear(); Clr();
while(now) V.push_back(now%),now/=;
printf("%d\n",fk(ans-dfs(V.size(),,,).sum2+mo));
}
return ;
}

HUD 4507 吉哥系列故事——恨7不成妻的更多相关文章

  1. HDU - 4507 - 吉哥系列故事——恨7不成妻(数位DP,数学)

    链接: https://vjudge.net/problem/HDU-4507 题意: 单身! 依然单身! 吉哥依然单身! DS级码农吉哥依然单身! 所以,他生平最恨情人节,不管是214还是77,他都 ...

  2. HDU 4507 吉哥系列故事――恨7不成妻(数位DP+结构体)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4507 题目大意:如果一个整数符合下面3个条件之一,那么我们就说这个整数和7有关 1.整数中某一位是7: ...

  3. HDU 4507 吉哥系列故事——恨7不成妻

    需要推下平方和的式子..维护个数,和,平方和. #include<iostream> #include<cstdio> #include<cstring> #inc ...

  4. HDU 4507 吉哥系列故事——恨7不成妻 (数位DP)

    题意: 如果一个整数符合下面3个条件之一,那么我们就说这个整数和7有关: 1.整数中某一位是7: 2.整数的每一位加起来的和是7的整数倍: 3.这个整数是7的整数倍: 给定一个区间[L,R],问在此区 ...

  5. 吉哥系列故事——恨7不成妻(数位DP)

    吉哥系列故事——恨7不成妻 http://acm.hdu.edu.cn/showproblem.php?pid=4507 Time Limit: 1000/500 MS (Java/Others)   ...

  6. [HDU4507]吉哥系列故事——恨7不成妻

    [HDU4507]吉哥系列故事--恨7不成妻 试题描述 单身!依然单身!吉哥依然单身!DS级码农吉哥依然单身!所以,他生平最恨情人节,不管是214还是77,他都讨厌!吉哥观察了214和77这两个数,发 ...

  7. 吉哥系列故事——恨7不成妻(数位dp)

    吉哥系列故事--恨7不成妻 传送门 Problem Description 单身! 依然单身! 吉哥依然单身! DS级码农吉哥依然单身! 所以,他生平最恨情人节,不管是214还是77,他都讨厌! 吉哥 ...

  8. Day9 - J - 吉哥系列故事——恨7不成妻 HDU - 4507

    单身! 依然单身! 吉哥依然单身! DS级码农吉哥依然单身! 所以,他生平最恨情人节,不管是214还是77,他都讨厌! 吉哥观察了214和77这两个数,发现: 2+1+4=7 7+7=7*2 77=7 ...

  9. 吉哥系列故事――恨7不成妻 HDU - 4507

    题目: 单身! 依然单身! 吉哥依然单身! DS级码农吉哥依然单身! 所以,他生平最恨情人节,不管是214还是77,他都讨厌! 吉哥观察了214和77这两个数,发现: 2+1+4=7 7+7=7*2 ...

随机推荐

  1. tracert命令与tracert (IP地址)-d有什么区别?

    他们的意义基本相同,都是路由追踪,返回从源到目标的路由情况:但tracert -d不解析各路由器的名称,只返回路由器的IP地址.而tracert 不仅返回各路由器的IP地址,而且返回其名称.简单来说, ...

  2. 安装conda后取消命令行前出现的base,取消每次启动自动激活conda的基础环境, 使用ubuntu 自带的Python环境

    方法一: 退出base环境回到系统自带的环境 conda deactivate 方法二 1,通过将auto_activate_base参数设置为false实现: conda config --set ...

  3. JVM----双亲委派模型

    加载类的开放性 我们在了解双亲委派模型之前,不得不先了解一下什么是类加载器.虚拟机设计团队之初是希望类加载过程“通过一个类的全限定名来获取描述该类的二进制字节流”这个动作能放到虚拟机外部实现,以便于让 ...

  4. kotlin array

    (1)Kotlin语言使用Array表示数组. (2)[] 可以用于访问数组的元素, [] 被进行了操作符的重载,调用的是 Array 类的 setter 和 getter 方法 2.创建数组 (1) ...

  5. Fragment入门代码

    让一个activity和加载多个布局文件 package com.example.fragment; import android.app.Activity; import android.app.F ...

  6. 利用phpStudy 探针 提权网站服务器

    声明: 本教程仅仅是演示管理员安全意识不强,存在弱口令情况.网站被非法入侵的演示,请勿用于恶意用途! 今天看到论坛有人发布了一个通过这phpStudy 探针 关键字搜索检索提权网址服务器,这个挺简单的 ...

  7. RabbitMQ学习之:(三)第一个RMQ的程序 (转贴+我的评论)

    RabbitMQ for Windows: Building Your First Application Posted by Derek Greer on March 7, 2012 This is ...

  8. 小D课堂 - 新版本微服务springcloud+Docker教程_4-01 常用的服务间调用方式讲解

    笔记 第四章 服务消费者ribbon和feign实战和注册中心高可用 1.常用的服务间调用方式讲解     简介:讲解常用的服务间的调用方式 RPC:             远程过程调用,像调用本地 ...

  9. ping指定IP的指定端口号

    ping是简单的测试网络连接情况的小工具,对于一般用户很实用,但是ping有个缺点就是,不能指定端口 tcping命令是针对tcp监控的,也可以看到ping值,即使源地址禁ping也可以通过tcpin ...

  10. ES6拓展符修改对象

    // ES6 拓展符合并两个对象let ab = { ...a, ...b }; // 等同于 let ab = Object.assign({}, a, b); // 修改对象部分属性.用户自定义的 ...