HDU - 1556 - Color the ball ( 差分 )

关注

原创最新推荐文章于 2026-03-21 01:39:59 发布 649 阅读 AI 写同款· GEO 优化 ›

碳酸钙岩石酸蚀酸溶COMSOL建模：非均质对酸蚀通道的影响多孔介质中的流动与化学反应耦合是石油工程和地下资源开发中的核心问题，酸液与岩石矿物的溶蚀过程直接关系到酸化增产效果。在碳酸盐岩储层中，盐酸与方解石的反应速率极快，酸液会优先沿高渗透区域推进，形成指状酸蚀通道，即wormhole。传统均质模型难以描述这种局部溶蚀行为，而数值模拟工具COMSOL通过耦合达西定律、稀物质传递与孔隙度演化方程，能够定量刻画酸蚀路径与渗透率的变化规律。工程实践中，模拟结果可用于预测酸液突破体积、优化注入排量，并为基质酸化方案设计提供依据。本文围绕碳酸钙岩石酸蚀酸溶COMSOL模型，探阅读详情

题目链接：点击进入

题目

在这里插入图片描述

思路

每次 [ a , b ] 颜色次数+1，这不典型区间修改，最后求每个气球颜色修改次数，区间操作即对差分数组操作，对差分数组求前缀和就是原数组，因此此题差分。

代码

// Problem: Color the ball
// Contest: Virtual Judge - HDU
// URL: https://vjudge.net/problem/HDU-1556
// Memory Limit: 32 MB
// Time Limit: 3000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

//#pragma GCC optimize(3)//O3
//#pragma GCC optimize(2)//O2
#include<iostream>
#include<string>
#include<map>
#include<set>
//#include<unordered_map>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstring>
#include<stack>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<fstream>
#define X first
#define Y second
#define best 131 
#define INF 0x3f3f3f3f3f3f3f3f
#define pii pair<int,int>
#define lowbit(x) x & -x
#define inf 0x3f3f3f3f
//#define int long long
//#define double long double
//#define rep(i,x,y) for(register int i = x; i <= y;++i)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const double pai=acos(-1.0);
const int maxn=1e6+10;
const int mod=1e9+7;
const double eps=1e-9;
const int N=5e3+10;
/*--------------------------------------------*/
inline int read()
{
    int k = 0, f = 1 ;
    char c = getchar() ;
    while(!isdigit(c)){if(c == '-') f = -1 ;c = getchar() ;}
    while(isdigit(c)) k = (k << 1) + (k << 3) + c - 48 ,c = getchar() ;
    return k * f ;
}
/*--------------------------------------------*/

int n,c[maxn],dis[maxn];

int main() 
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);cout.tie(0);
	while(cin>>n)
	{
		if(n==0) break;
		for(int i=1;i<=n;i++)
		{
			int l,r;
			cin>>l>>r;
			dis[l]++;
			dis[r+1]--;
		}
		for(int i=1;i<=n;i++)
		{
			dis[i]=dis[i]+dis[i-1];
			if(i==n)
				cout<<dis[i]<<endl;
			else
				cout<<dis[i]<<' ';
		}
		for(int i=1;i<=n+1;i++) dis[i]=0;
	}
    return 0;
}