编程马拉松17C

最新推荐文章于 2025-06-18 20:28:59 发布

lovewangtaotao

最新推荐文章于 2025-06-18 20:28:59 发布

阅读量3k

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.csdn.net/lovewangtaotao/article/details/52369826

本文探讨了一个涉及模数运算的复杂路径计数问题，通过预处理质数及其符号来优化计算，并采用动态规划方法解决给定图中特定路径的数量问题。文中详细介绍了算法实现细节，包括初始化过程、路径矩阵的更新以及最终结果的计算。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include <algorithm>
#include <assert.h>
#include <complex>
#include <ctype.h>
#include <functional>
#include <iostream>
#include <limits.h>
#include <locale.h>
#include <map>
#include <math.h>
#include <queue>
#include <set>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <string>
#include <time.h>
#include <vector>
//#include <unordered_set>
//#include <unordered_map>

#pragma warning(disable:4996)
using namespace std;

#define mp make_pair
typedef long long ll;
typedef unsigned long long ull;
typedef double db;
typedef long double ldb;
typedef pair <int, int> pii;
typedef pair <ll, ll> pll;
typedef pair <ll, int> pli;
typedef pair <ldb, ldb> pdd;

int IT_MAX = 1 << 21;
const int MOD = 1000000007;
const int INF = 1034567891;
const ll LL_INF = 1234567890123456789ll;
const db PI = 3.141592653589793238;
const ldb ERR = 1E-12;
#define mem(x ,y ) memset(x , y , sizeof(x))
struct road{
    int x,y,w ;
}edg[50010];
ll w[105][105] ;
ll d[105][105][105] ;
ll sum[105] ;
const int mod = 1e9+7 ;
int n , m ;
void add(int s,int e,int w,int t){
    for(int z=1;z<=sqrt(w);z++)if(w % z == 0){
        int tmp1 = z , tmp2 = w / z ;
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                ll sav1 = t * (d[tmp1][i][s] * d[tmp1][e][j] ) ;
                ll sav2 = t * (d[tmp2][i][s] * d[tmp2][e][j] ) ;
                d[tmp1][i][j] += sav1 ;
                if(tmp1 != tmp2) d[tmp2][i][j] += sav2 ;
                sum[tmp1] += sav1 ;
                if(tmp1 != tmp2) sum[tmp2] += sav2 ;
            }
        }
    }
}
int u[105] ;
int v[105] ;
void init(){
    for(int i=1;i<=100;i++) v[i] = 0 , u[i] = 1 ;
    u[1] = 1 ;
    for(int i=2;i<105;i++)if(!v[i]){
        for(int j=i;j<105;j+=i){
            v[j] = 1 ;
            if(j % (i*i) == 0) u[j] = 0 ;
            else u[j] *= -1 ;
        }
    }
}
int main() {
    init() ;
    cin >> n >> m ;
    for(int i=0;i<m;i++){
        scanf("%d%d%d" , &edg[i].x , &edg[i].y , &edg[i].w) ;
    }
    mem(d , 0) ;
    mem(sum , 0) ;
    ll ans = 0 ;
    for(int x=1;x<=100;x++){
        for(int i=1;i<=n;i++) d[x][i][i] = 1 ;
        for(int t=0;t<m;t++){
            if(edg[t].w % x == 0) d[x][edg[t].x][edg[t].y] ++ ;
        }
        for(int k=1;k<=n;k++){
            for(int i=1;i<=n;i++){
                for(int j=1;j<=n;j++)if(i!=j&&k!=j&&i!=k){
                    d[x][i][j] += d[x][i][k] * d[x][k][j] ;
                    d[x][i][j] %= mod ; 
                }
            }
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++)if(i!=j){
                sum[x] += d[x][i][j] ;
                sum[x] %= mod ; 
            }
        }
        ans += sum[x] * u[x] ;
    }
    cout << ((ans % mod)+mod)%mod << endl ;
    int Q ; cin >> Q ;
    int ta , tb ;
    for(int q=0;q<Q;q++){
        cin >> ta >> tb ; ta -- ;
        int ts = edg[ta].x ; int te = edg[ta].y ;
        add(ts , te , edg[ta].w , -1) ;
        add(ts , te , tb , 1) ;
        edg[ta].w = tb ;
        ans = 0 ;
        for(int i=1;i<=100;i++){
            ans += sum[i] * u[i] ;
        }
        cout << ((ans % mod)+mod)%mod << endl ;
    }
}

用题解的方法还是没过，现在写在这，以后可以看代码，

我用的题解方法。不取mod不超时。取mod就超时。日了狗。以后再看啊