优化到极限 好点

题目大意:有一个N*M的网格,左下角的坐标是(0, 0) 右上角的坐标是(m, n), 即一共有(N+1)*(M+1)个整点, 在任意两个整点连一条线段,若没有其他整点在这条线段上,且该线段长度在[L...H]之间,那么这两个整点叫做好点,求有多少对好点。答案mod B。1 <= M, N <= 100000鄙视

首先要解决的一个问题,微笑若给出两个整点,要判断它们是不是好点。这很简单,只要看它们的横坐标的差和纵坐标的差是不是互质的(为什么,自己想想)。

暴力的话,枚举每两两个整点,时间复杂度 O (n ^ 4)。

通过观察,当纵坐标的差值是一定的,那么我们先枚举每个纵坐标差h,再统计横坐标的(枚举上面的每每两个点),那么乘上(n - h + 1), 时间复杂度 O(n ^ 3)。

既然纵坐标能够优化,那么横坐标也可以,枚举每个纵坐标的差h和横坐标的差w,若合法,则加上(n - h + 1) *(m - w + 1),时间复杂度 O(n ^ 2)。

通过一系列比较简单的优化,时间已经降下来不少,但还是远远不够尴尬

我们可以发现,当h确定时,只要w和h是互质的,就是合法的(可以用所有状态减去不合法的,也就是于h不互质的),这可以用容斥原理来算出来,使用一下等差数列(这个要自己画画图才能思路清晰),就差不多了。时间复杂度O(n * (2 ^ 8(最多只会有8个不同的质因数)))

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;

typedef long long ll;
#define Q 100000

ll f[Q + 7];
ll n, m, L, H, mod;

void Init() {
    cin >> n >> m >> L >> H >> mod;
    
    memset(f, 0, sizeof(f));
    
    //算出所有数的最小质因数 
    for (int i = 2; i <= Q; i ++)
        if (f[i] == 0)
            for (int j = i; j <= Q; j += i) 
                if (f[j] == 0) f[j] = i;
}

// 因为有负数的情况,取模要注意 
inline ll MOD(ll t) {
    return ((t % mod) + mod) % mod;
}

ll Solve(ll len) {
    if (len == 0LL) return 0LL;
    // 线段长度为1的。 
    ll res = MOD((n + 1) * m + (m + 1) * n);
    
    ll w = 0LL;
    // 枚举高度 
    for (ll h = n; h > 0; h --) {
        // 宽度是不下降的 
        while (w < m && h * h + (w + 1) * (w + 1) <= len) w ++;
        if (w == 0LL) continue;
        
        int p = 0;
        ll t = h, que[10];
        
        // 分解质因数 
        for (; t > 1; p ++) {
            que[p] = f[t];
            for (; t % que[p] == 0; t /= que[p]);
        }
        
        // 所有状态 
        ll sum = MOD(w * (m - w + 1) + w * (w - 1) / 2);
        
        // 统计不合法状态 
        for (int set = 1; set < (1 << p); set ++) {
            ll num = 1LL, type = 1LL;
            for (int i = 0; i < p; i ++)
                if ((set & (1 << i)) > 0) {
                    num *= que[i];
                    type *= -1LL;
                }
            
            sum += (w / num) * (m - w + 1L) * type;
            sum = MOD(sum);
            
            t = w / num;
            sum += ((t - 1) * t / 2LL * num + (w - t * num) * t) * type;
            sum = MOD(sum);
        }
        
        res += (sum << 1) * (n - h + 1);
        res = MOD(res);
    }
    
    return res;
}

int main() {
    freopen("GoodPoint.in", "r", stdin);
    freopen("GoodPoint.out", "w", stdout);
    
    Init();
    
    // 分成两个问题解决 
    cout << MOD(Solve(H * H) - Solve(L * L - 1)) << endl;;
    
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 4、下载使用后,可先查看README.md文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值